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Résumé : 


L'étude statique ou dynamique des poutres et des plaques reposant sur la surface ou à l’intérieur 
des milieux élastiques reste un sujet d’actualité vu l’utilisation accrue de ce genre de structures 
dans les différents domaines de construction, notamment celui du génie civil, génie mécanique, 
secteur naval, secteur aéronautique etc. Ce genre d’étude à réaliser doit être accomplie de façon 
authentique afin d’aboutir, d’une manière rationnelle et précise, au maximum du comportement 
réel statique et dynamique des structures en interaction avec des milieux élastiques. Dans ce 
contexte, ce travail de thèse de doctorat est focalisé sur une approche semi-analytique, basant 
sur le couplage de deux études essentielles relatives à la structure et au milieu élastique, 
permettant l’analyse dynamique des plaques rectangulaires reposant sur la surface d’un semi- 
infini élastique de propriétés inertielles (modèle de Lamb). 

L’approche utilisée est basée sur une méthode de calcul connue dans la théorie d’élasticité sous 
le nom ‘la méthode mixte”, et connue aussi dans certaines littératures sous le nom ‘la méthode 
de Zhemochkin’. Le principe de cette méthode pour notre étude consiste à discrétiser le système 
étudié (plaque rectangulaire reposant sur la surface du semi-infini élastique de propriétés 
inertielles) en un nombre fini d’éléments rectangulaires identiques. Le contact continu entre la 
plaque et la surface du semi-infini élastique est donc remplacé par un contact partiel assuré par 
des liaisons se trouvant au niveau des centres des éléments de discrétisation. 


L’application de l’approche proposée, basée sur la méthode mixte, aboutit finalement à un 
système d’équations canoniques linéaire contenant plusieurs paramètres relatifs à la structure 
étudiée. Les deux paramètres essentiels du système d’équations canoniques représentent les 
déplacements verticaux de la plaque et les déplacements de la surface du semi-infini élastique 
sur laquelle la plaque se repose. La méthode énergétique de Ritz est utilisée pour déterminer les 
déflexions de la plaque étudiée, alors que les déplacements verticaux de la surface du semi- 
infini élastique sont déterminés par l’étude de la fonction de Green pour le modèle de Lamb. 
L’étude de la fonction de Green est la tâche la plus compliquée et la plus difficile 
mathématiquement de ce travail de doctorat. Sa complexité mathématique réside dans la 
détermination de plusieurs intégrales dont certaines sont associées au problème de singularité 
ce qui représente de plus le grand challenge à surmonter de cette étude. Ce challenge est 
surmonté grâce à une attention particulière investie pour l’évaluation de ces intégrales 
analytiquement ou semi analytiquement en jouant sur certains changements de variable et sur 
certains remplacements des fonctions spéciales par des séries polynomiales. Le second grand 
défi rencontré dans cette étude est la détermination des déflexions de la plaque qui surmonté 
par l’application de la méthode énergétique de Ritz en adoptant le principe de la solution de 
Clebsch pour la résolution de l’équation différentielle aux dérivées partielles des déflexions de 
la plaque étudiée. Enfin, le couplage de ces deux études essentielles avec d’autres études 
auxiliaires relatives aux caractéristiques géométriques et mécanique de la structure étudiée est 
accompli par la méthode de Zhemochkin. Après plusieurs simplifications et transformations 
mathématiques, la méthode de résolution finale est présentée sous une forme matricielle 
permettant la détermination des efforts de liaison dans la zone de contact et puis par 


—————— 
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l’application des différentes formules de la théorie d’élasticité, les autres entités physiques sont 
déterminées à savoir les valeurs des fréquences propres de la plaque, les modes propres de la 
plaque, sa réponse dynamique à différents types d’excitations extérieures etc. 


Afin de s’assurer de la fiabilité, l’authenticité et la précision des résultats de calcul donnés par 
cette approche, une comparaison des valeurs des fréquences propres de la plaque et de sa 
réponse dynamique due à une excitation harmonique avec le modèle du semi-infini élastique 
de propriétés distributives (modèle de Boussinesq) est faite. Aussi les valeurs des fréquences 
propres de la plaque sont vérifiées et authentifiées par l’application du principe de superposition 
de la méthode modale. Finalement, les fréquences d’excitation de la plaque sont aussi localisées 
en utilisant l’analyse (FFT) des spectres obtenus à partir des signaux de la réponse dynamique 
de la plaque. 


Mots-clés : Problèmes de Contact ; Plaque ; Modèle de Lamb ; Fonction de Green ; Méthode 
de Ritz ; Fréquences propres ; Modes propres ; Réponse dynamique. 
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Abstract : 


The static or dynamic study of beams and plates resting on the surface or inside elastic 
foundation is a topical subject given the large use of this kind of structures in the various 
construction fields like civil engineering, mechanical engineering, naval sector, aeronautical 
sector etc. The researched study must be authentic in order to achieve and predict with a rational 
and precise manner as much as possible the real static and dynamic behavior of structures in 
interaction with elastic foundation. 

This doctoral thesis focuses on a semi-analytical approach, based on the coupling of two 
essential studies relating the structure and the elastic foundation, thereby allowing the dynamic 
analysis of a rectangular plate resting on the surface of an elastic half-space with inertial 
properties (Lamb's model). 

The coupling method used is known in elasticity theory as the mixed method, also referred to 
as Zhemochkin's method. The principle of this method for consists of discretizing the studied 
system (1.e. the rectangular plate resting on the surface of the elastic half-space with inertial 
properties) into a number of identical rectangular elements. The continuous contact between 
the plate and the surface of the elastic half-space is replaced by a partial contact provided by 
connections located at the centers of the elements. The application of the proposed approach, 
based on the mixed method, leads to a linear canonical system of equations containing several 
parameters relating to the studied structure. The two essential parameters of the canonical 
system represent the vertical deflections of the plate and the vertical displacements of the elastic 
half-space surface on which the plate rests. The energy method of Ritz's is used to determine 
the deflections of the studied plate, while the vertical displacements of the surface of the elastic 
half-space are determined by the study of Green's function for the Lamb's model. The study of 
Green's function is the most complicated and the most mathematically difficult task of this 
doctoral work. Its mathematical complexity resides in the determination of several integrals, 
some of which are associated with the singularity phenomena. The singularity problem 
represents the greatest challenge to be overcome in this study. 

This challenge is overcome through special attention to the evaluation of these integrals 
analytically or semi-analytically by playing on certain points like changes of variable systems 
and certain approximations of special functions by polynomial series. The second major 
challenge encountered in this study is the determination of the plate deflections. This 1s 
addressed by the application of the energy method of Ritz after adopting Clebsch's solution of 
partial differential equation concerning the deflections of the studied plate. Finally, these two 
essential studies are coupled with other auxiliary studies relating to the geometrical and 
mechanical characteristics of the studied structure. This is accomplished with Zhemochkin's 
method. After several simplifications and mathematical transformations, the final resolution 1s 
presented in a matrix form, thereby allowing the determination of the reactive forces in the 
contact zone. Next, through the application of various formulations of the theory of elasticity, 
other physical magnitudes are determined including the plates’ Eigen-frequencies, Eigen- 
shapes, and dynamic response under various external excitations. 
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To ensure the reliability and the precision of the computational results, a comparison of the 
values of the Eigen-frequencies of the plate and its dynamic response due to a harmonic external 
excitation are accomplished through a comparison with Boussinesq's model of the distributive 
properties of an elastic half-space of. Additionally, the values of the Eigen-frequencies of the 
plate are authenticated by applying the superposition principle of the modal method. Finally, 
the excitation frequencies of the plate are localized using a Fast Fourier Transform analysis of 
the spectra obtained from the signals of the dynamic response of the plate. 


Keywords : Contact Problems ; Plate ; Lamb's model ; Green's function ; Ritz's method ; Eigen- 
frequencies ; Eigen-shapes ; Dynamic response. 
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Introduction générale : 


Les études dynamiques ou statiques de poutres, plaques et autres types de structures reposant 
sur les surfaces des milieux élastiques de différents types appartiennent au domaine des 
problèmes du contact. Ce vaste domaine des problèmes de contact est défini comme étant 
l'interaction d’au moins deux systèmes physiques de caractéristiques mécaniques et 
géométriques différentes. Son importance réside dans son vaste domaine d’application comme 
le cas des plaques et poutres représentant des éléments structuraux de base de différentes 
structures notamment dans le domaine du génie mécanique, génie civil et autres. 


Cependant, l’analyse statique ou dynamique de ce genre de structures n’est pas achevée 
complètement à nos jours, car elle est liée à de très grandes difficultés mathématiques. Ceci est 
expliqué par leur nature compliquée ce que l’on peut exprimer par les facteurs suivants : la 
diversité des hypothèses simplificatrices ; le nombre important de paramètres géométriques et 
mécaniques entrant dans le calcul ; les difficultés de modélisation des différentes parties de la 
structure étudiée ; les difficultés de couplage des différentes études liées à ce genre de structures 
etc. Pour ces raisons essentielles, l’analyse dynamique plus précisément des structures en 
interaction avec milieux élastiques n’est pas accomplie d’une manière satisfaisante et définitive 
jusqu’à présent. Dans ce contexte, la nécessité de développer les approches de calculs déjà 
existantes ou de proposer de nouvelles approches plus précises et plus efficaces est 
indispensable et s’impose fortement. 


Dans ce contexte, ce travail préparé dans le cadre de la thèse de doctorat est accompli par une 
approche semi-analytique permettant le calcul dynamique des plaques rectangulaires reposant 
sur la surface d’un milieu élastique de propriétés inertielles. L'approche proposée se base sur la 
méthode de calcul des structures hyperstatiques connue dans la théorie d’élasticité par la 
méthode mixte et connue aussi dans certaines références bibliographiques par la méthode de 
Zhemochkin. L'approche proposée est aménagée de telle sorte de déterminer en premier lieu 
les forces de réaction dans la zone de contact entre la plaque et la surface du semi-infini 
élastique de propriétés inertielles. Dès que les forces de réaction sont déterminées, l’application 
des différentes formules de la théorie d’élasticité faisant relation entre différentes entités 
physiques de la structure étudiée permettra la détermination des restes des entités physiques à 
savoir : les fréquences propres et les modes propres de la plaque, sa réponse dynamique due à 
de différents types de charges harmoniques extérieures etc. 


Pour alléger un peu la tâche de calcul très compliquée, l'amortissement et les forces de 
frottement dans la zone de contact entre la plaque et la surface du milieu élastique sont négligés 
et que la plaque étudiée est considérée de forme carrée. 


Le travail préparé dans le cadre de cette thèse de doctorat est structuré en : un sommaire, un 
résumé en français, un résumé en anglais et un résumé en arabe, une introduction générale, cinq 
chapitres, une conclusion générale, des références bibliographiques et quelques annexes. 
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Chapitre 1 est consacré à l’étude bibliographique sur la thématique de la thèse couvrant 
beaucoup plus les problèmes de contact. Une synthèse bibliographique des différents travaux 
de recherche effectuée dans ce domaine, ainsi que certains rappels de base concernant la théorie 
des plaques sont aussi apportés dans ce chapitre. Aussi bien ce chapitre est enrichi par une 
description de quelques modèles des milieux élastiques connus dans la littérature et une 
conclusion du chapitre. 


Chapitre 2 est consacré à la détermination des déflexions de la plaque de forme carrée en 
utilisant la solution de Clebsch de l’équation différentielles des déflexions de la plaque et la 
méthode énergétique de Ritz. 


Chapitre 3 est consacré à l’évaluation des intégrales lors de l’étude de la fonction de Green 
définissant les déplacements verticaux de la surface du semi-infini élastique de propriétés 
inertielles (modèle de Lamb). Certaines de ces intégrales sont associées à des problèmes de 
singularité ce qui a nécessité des efforts particuliers pour surmonter ce challenge et obtenir 
l’évaluation appropriée de ces intégrales pour qu’elles puissent être appliquées à la structure 
étudiée. 

Chapitre 4 est consacré à la description en détail de l'approche utilisée pour le calcul dynamique 
des plaques rectangulaires reposant sur la surface du milieu élastique de propriétés inertielles. 


Chapitre 5 est consacré à l’application de l’approche proposée pour l’analyse dynamique d’une 
plaque de forme carrée reposant sur la surface du milieu élastique de propriétés inertielles. La 
modélisation finale de la structure étudiée a abouti à un système matriciel permettant la 
détermination des forces de réaction dans la zone de contact ce qui a permis ensuite de 
déterminer les fréquences propres et les modes propres de la plaque et sa réponse dynamique 
due à plusieurs types de charges harmoniques extérieures. Ce chapitre contient aussi la 
vérification des valeurs des fréquences propres de la plaque en utilisant le principe de 
superposition de la méthode modale. La comparaison des valeurs des fréquences propres de la 
plaque et la comparaison de sa réponse dynamique sont effectuées en considérant la même 
structure reposant sur la surface du semi-infini élastique de propriétés distributives (modèle de 
Boussinesq). 
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CHAPITRE I : Etat de l’art et étude bibliographique de la thématique de la thèse 
L1. Introduction : 


De nombreux sujets de grande importance dans la pratique industrielle ou civile peuvent être 
liés à l’utilisation des plaques et des poutres en interaction avec des milieux élastiques. Les 
chaussées en béton armé des autoroutes et des pistes des aéroports, les tôles des différentes 
structures mécaniques et les coques réservoirs etc. ce sont parmi des centaines d’autres 
applications bien connues. L’analyse statique ou dynamique des plaques et d’autres types de 
structures en interaction avec des milieux élastiques est également basée sur l'étude de la 
fonction de Green qui décrit en détail le type et la nature du milieu élastique considéré. 

Dans ce contexte, ce chapitre est consacré aux trois volets essentiels suivants : (1) synthèse de 
références bibliographiques des travaux de recherche relatifs à ce sujet ; (11) certains rappels de 
base concernant la théorie des plaques et poutres ; (1ii) la description détaillée des modèles des 
milieux élastiques connus dans la littérature. Les modèles les plus utilisés sont : le modèle des 
ressorts (Winkler), le modèle du semi-infini élastique de propriétés distributives (Boussinesq), 
le modèle du semi-infini élastique de propriétés inertielles (Lamb), le modèle du semi-infini 
élastique en couches etc. La description de chaque modèle est accompagnée, non seulement, 
par la définition des propriétés mécaniques du modèle et sa conception pratique, mais aussi 
avec des expressions mathématiques du noyau essentiel définissant le problème et les difficultés 
de solution qui se résultent avec quelques recommandations et orientations de solution. 


L2. Synthèse des références bibliographiques : 


L’étude statique ou dynamique de n'importe quelle structure reposant sur la surface d’un milieu 
élastique est liée à des difficultés mathématiques majeures. Ce challenge est expliqué par 
plusieurs facteurs qui entrent dans l’analyse de telles structures. La diversité des hypothèses 
liées à ce domaine des problèmes de contact est un enjeu important et bien évidemment plus 
que le nombre des hypothèses prises est grand, plus que l’erreur de la solution trouvée du 
problème étudié se réduit, mais avec un effort de calcul plus important. La diversité des 
paramètres géométriques et mécaniques, ainsi que la combinaison de plusieurs différentes 
études entrant dans le calcul de ce genre de structures influe considérablement sur la complexité 
de sa solution. Par exemple, le couplage des études d’au moins deux différents types de 
structures comme le cas de la plaque reposant sur la surface du milieu élastique de propriétés 
inertielles s’impose. C’est pourquoi, il est impératif de faire une lecture approfondie de la 
littérature dans ce vaste domaine des problèmes de contact. 


Une solution proche de la forme exacte est donné dans le travail de X. Feng et al [1] concernant 
la détermination du déplacement de la surface d'un semi-infini élastique de modèle de Lamb 
causé par une charge se trouvant à une profondeur de la surface du milieu élastique. Cette étude 
représente dans son ensemble l’analyse de la fonction de Green pour le modèle de Lamb en 3D 
et la solution proposée est bien vérifiée par la comparaison de ses résultats numériques avec la 
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solution de Johnson. Les formules apportées dans cette étude peuvent constituer une solution 
précieuse qui servira comme référence aux autres solutions numériques pour comparaison. Kh. 
Kh. Il’yasov et al [2] ont étudié la solution du problème de Lamb concernant le cas d’une charge 
extérieure harmonique uniformément répartie agissant sur la surface du semi-infini élastique. 
Ce type de charge tient en considération aussi le paramètre de Poisson du noyau du semi-infini 
élastique faisant partie des paramètres de la structure étudiée en vue de trouver une solution 
avec beaucoup d’exactitude. La solution finale est obtenue sur la base des transformations des 
intégrales et aussi avec l’utilisation de la méthode des éléments finis ce qui a abouti à une 
solution analytique du problème des vibrations en régime permanent sous une forme simple. 


Le travail de V. V. Nishawala et al apporté dans la référence [3] rapporte une généralisation du 
problème de Lamb à un semi-infini élastique en considérant les champs aléatoires (RES) de 
densité de masse, soumis à une charge linéaire normale. Les deux champs RFSs corrélées et non 
corrélées sont considérés, sans aucune restriction de bruit faible. Les automates cellulaires (CA) 
sont utilisés pour simuler la propagation des ondes. CA est définie comme étant une méthode 
de calcul local pour la discrétisation rectangulaire du domaine spatial ce qui équivaut à 
appliquer la méthode des différences finies pour résoudre les équations différentielles de 
l'élasticité classique. Enfin, les auteurs de cette étude ont déterminé dans quelle mesure la 
fractale ou le paramètre de Hurst est un facteur significatif dans la modification de la solution 
du problème de Lamb concernant le cas stochastique planaire par l’évaluation du coefficient de 
variation de la réponse par rapport au coefficient de variation des champs aléatoires REs. La 
réponse de la surface du semi-infini élastique dû à l'impact d'un pénétrateur normal à cette 
surface est considérée dans le travail de I. Argatov et al [4]. Les déplacements des points situant 
sur la surface du semi-infini élastique éloignés du point où le pénétrateur est appliqué sont 
décrits en résolvant le problème de Lamb pour ce problème. Dans le cas où les chocs par un 
pénétrateur sont transférés à la surface du semi-infini élastique par une plaque flexible, un 
modèle asymptotique du premier ordre est appliqué qui tient compte des propriétés dissipatives 
élastique dues au transfert d'énergie élastique vers le semi-infini élastique. Les principaux 
résultats de ce travail sont obtenus sous forme approchée. L’article publié par L. À Khajiyeva 
et al [5] est focalisé essentiellement à dériver le modèle asymptotique pour la propagation des 
ondes de surface du semi-infini élastique incompressible précontraint, soumis à une charge de 
surface bien définie. Dans cette étude de modèle dérivé, la désintégration de la surface est régie 
par une équation elliptique pseudo-statique, tandis que la propagation des ondes est décrite par 
une équation hyperbolique sur la surface et l'effet de la précontrainte, à savoir la contrainte de 
Cauchy principale o2, est étudié aussi. Enfin, un exemple illustratif du problème de Lamb est 
considéré, démontrant l'efficacité de l'approche apportée. 


L’étude des vibrations libres d’une plaque rectangulaire isotrope avec des conditions aux limites 
arbitraires basées sur le principe de la matrice de rigidité dynamique est présentée dans le travail 
M. Nefovska-D et al [6]. Dans ce travail, la méthode de superposition de Gorman est exploitée 
pour obtenir la solution des équations de mouvement dominantes et la matrice de rigidité 
dynamique est dérivée à l'aide de la méthode de projection. Les résultats obtenus dans cette 
étude sont en bonne concordance avec les résultats obtenus en utilisant des solutions exactes 
pour certains cas particuliers, ainsi qu'avec des solutions données par la méthode des éléments 
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finis. Le calcul des plaques flexibles reposant sur la surface des milieux élastiques en utilisant 
les équations généralisées de la méthode des différences finies est donné par M. Sali et al [7]. 
L'algorithme utilisé permet de prendre en compte, non seulement, les discontinuités finies de la 
fonction recherchée, mais aussi de sa première dérivée et de la partie droite de l'équation 
différentielle sans l’introduction des points fictifs et aussi d'un resserrement particulier du 
maillage. Les exemples présentés montrent la précision des résultats et la simplicité de 
l'algorithme utilisé. Les réponses dynamiques des plaques orthotropes reposant sur la surface 
des milieux élastiques modélisées numériquement à l'aide de la méthode de quadrature 
différentielle sont étudiées par M.-H. Hsu [8]. Dans cette étude, la technique de quadrature 
différentielle est utilisée pour transformer des équations différentielles partielles en un 
problème de valeurs propres discrètes. Les résultats numériques obtenus démontrent que la 
rigidité de la fondation à un impact significatif sur le comportement dynamique des plaques 
orthotropes. 


Une solution exacte est établie dans le travail de M. À. Foyouzat et al [9] concernant le problème 
de la vibration libre non amortie d'une plaque mince circulaire reposant sur la surface du milieu 
élastique de type de Winkler. La procédure de résolution est démontrée par un exemple 
illustratif, dans lequel l'équation de fréquence générale est dérivée pour deux conditions aux 
limites différentes. La solution est également vérifiée à l'aide des résultats disponibles dans la 
littérature et la méthode de solution proposée est directement applicable aux problèmes des 
plaques circulaires reposant sur les surfaces des milieux élastiques de type Pasternak à module 
variable. L’étude apportée par Z. Lu et al. [10] concerne les caractéristiques vibratoires d'une 
plaque reposant sur la surface d’un milieu élastique à deux paramètres et soumise à des charges 
rectangulaires mobiles avec des vitesses variables. Dans cette étude, la solution générale pour 
les déflexions dynamiques de la plaque est dérivée en utilisant la transformée double de Fourier. 
Les effets de la vitesse de charge initiale, de l'accélération de la charge, de la décélération de la 
charge et de la résistance horizontale au fond de la plaque sur la déflexion dynamique sont 
discutés. La vitesse initiale, l'accélération et la décélération de la charge rectangulaire 
influencent sur la réponse dynamique, et en plus la déflexion dynamique de la plaque à la vitesse 
critique diminue considérablement. Dans l’article de M. Dehghan et al [11] une analyse des 
vibrations et de la stabilité de plaques rectangulaires épaisses reposant sur la surface des milieux 
élastiques est donnée. Dans cette étude, un modèle à deux paramètres (Pasternak) est considéré 
pour décrire la fondation élastique. Le problème des valeurs propres dans le domaine 3-D est 
résolu numériquement par une combinaison de la méthode des éléments finis et de la quadrature 
différentielle (DQM). La précision des résultats obtenus est validée en comparant avec quelques 
solutions analytiques de la littérature. Dans le travail de B. Akgoz et al [12], l'analyse des 
vibrations libres géométriquement non linéaires des plaques minces stratifiées reposant sur la 
surface des milieux élastiques non linéaires est faite. Le modèle de milieu élastique considéré 
est de Winkler-Pasternak. La méthode de convolution singulière discrète est utilisée pour 
obtenir les équations discrétisées du mouvement des plaques. Les effets de la géométrie des 
plaques, des conditions aux limites, des propriétés des matériaux et des paramètres du milieu 
élastique sur le comportement vibratoire non linéaire des plaques sont présentés aussi dans cette 
étude. 
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Une solution sous forme de série pour l’étude de la vibration transversale des plaques 
rectangulaires Mindlin avec des supports ponctuels élastiques autour de son contour est 
apportée dans le travail de F. Pang et al [13]. Dans cette étude, la solution en série du problème 
est obtenue en utilisant une méthode améliorée de série de Fourier, dans laquelle les 
déplacements de vibration et les rotations transversales du plan médian sont représentés par une 
double série de cosinus de Fourier et quatre fonctions supplémentaires. La convergence, la 
précision, la stabilité et l'efficacité de la méthode proposée ont été examinées à travers une série 
d'exemples numériques. Quelques exemples numériques sur la fréquence non dimensionnelle 
et les formes de mode de plaques rectangulaires Mindlin avec différentes conditions aux limites 
sont donnés. Le problème de contact concernant l'action de deux plaques de Kirchhoff semi- 
infinies reposant sur la surface de milieu élastique de type multicouches élastiques est considéré 
dans le travail de V. Babeshko et al [14]. Le problème est d'étudier les concentrations de 
contraintes de contact entre les plaques et la couche sur laquelle se repose la plaque. Cette étude 
permet aussi d’induire le paramètre de la destruction des matériaux au cas du séisme. M. Sobhy 
[15] a traité la vibration et le comportement au flambage d'une plaque sandwich de matériau à 
gradient exponentiel (EGM) reposant sur la surface des milieux élastiques avec diverses 
conditions aux limites. Le milieu élastique est modélisé comme type de Pasternak qui peut être 
isotrope ou orthotrope et si le cisaillement est négligé, la solution converge vers le type de 
Winkler. Les équations régissant la réponse dynamique des plaques composites non homogènes 
sont déduites en utilisant diverses théories des plaques de déformation par cisaillement. Les 
résultats numériques pour les fréquences propres et les charges de flambement critiques des 
plaques sandwich symétriques sont présentés. La validité de la solution est démontrée par 
comparaison avec les solutions disponibles dans la littérature. Dans l’article de E. V. 
Barmenkova et al [16], il est donné le calcul des plaques reposant sur la surface d’un milieu 
élastique de type Winkler d’un coefficient de réaction de fondation variable et constant. Le 
calcul des déflexions des plaques est fait par la méthode des éléments finis. Les résultats du 
calcul obtenus sont comparés avec ceux de la même plaque reposant sur milieux élastiques de 
différents modèles. 


La réponse dynamique d'une plaque élastique mince, isotrope et linéaire reposant sur la surface 
du semi-infini élastique et sollicitée par une charge mouvante sur sa surface est obtenue 
analytiquement dans le travail Y. Chen et al. [17]. Dans cette étude, il a supposé que la charge 
est répartie sur une longueur finie et se déplace à vitesse constante et la charge mobile est 
considérée sous une forme complexe de séries de Fourier impliquant la coordonnée horizontale, 
le temps et la vitesse. La solution trouvée est validée en comparant ses résultats à d'autres 
solutions analytiques se référant à des cas plus simples. Enfin, des études paramétriques sont 
menées pour évaluer l'effet d'anisotropie sur la réponse du système au chargement du 
mouvement pour différentes valeurs de porosité, perméabilité et vitesse de charge. L’instabilité 
dynamique due aux résonances paramétriques de plaques rectangulaires minces reposant sur la 
surface d’un milieu élastique sous l’effet des masses mobiles successives est étudiée dans le 
travail de E. Torkan et al [18] et considéré comme un problème linéaire périodique. Dans cette 
étude, l'équation différentielle partielle est transformée en un ensemble d'équations 
différentielles ordinaires en utilisant la méthode de Galerkin. Il a conclu que les résultats de son 
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étude sont utiles pour la conception de chaussées routières reposant sur des sols, des fondations 
des ponts et des chaussées des aéroports. La réponse dynamique des plaques à gradation 
fonctionnelle (FG) dans un environnement thermique sous une charge mobile est étudiée dans 
le travail de P. Malekzadeh et al [19]. La formulation est dérivée sur la base de la théorie de la 
déformation par cisaillement du premier ordre (FSDT), qui inclut les effets des contraintes 
thermiques initiales induites par l'environnement thermique, ainsi que du milieu élastique. La 
formulation et la méthode de solution sont validées par l’étude de leur comportement de 
convergence et par comparaison avec les résultats existants dans la littérature. Enfin, les 
influences de l'élévation de température, de l'indice gradué du matériau, de la vitesse de la 
charge mobile, des paramètres de fondation élastiques et des conditions aux limites sur le 
comportement dynamique des plaques FG dans un environnement thermique et soumises à une 
charge mobile sont présentées aussi dans cette étude. Le problème des réponses dynamiques 
d'un semi-infini poroélastique insaturé et soumis à une charge rectangulaire mobile est étudié 
de manière analytique et numérique dans le travail de Z. Lu et al [20]. En utilisant la transformée 
de Fourier, la solution générale des équations est dérivée dans le domaine transformé, puis un 
problème de valeurs limites correspondant est formulé. Les résultats numériques révèlent que 
les effets de paramètres comme saturation du sol, vitesse de chargement et fréquence 
d'excitation sur la réponse dynamique du sol non saturé sont significatifs. 


Le semi-infini élastique de modèle de Pasternak est utilisé par R. Buczkowski et al [21] pour 
étudier des plaques épaisses en statique et pour déterminer ses fréquences propres en utilisant 
la méthode des éléments finis. Dans le cas de l'analyse des fréquences propres, l'équation du 
mouvement est dérivée par l’application du principe de Hamilton impliquant la variation de 
l'énergie cinétique et potentielle de la plaque et du semi-infini élastique. La théorie de la 
déformation par cisaillement du premier ordre de la plaque est utilisée. Les résultats numériques 
de cette analyse sont vérifiés par l’utilisation des éléments minces basés sur la théorie classique 
de Kirchhoff et des éléments de plaques épaisses. La solution exacte concernant la vibration 
libre de la plaque circulaire élastique multicouches reposant sur la surface d’un milieu élastique 
de n’importe quel modèle est donnée dans le travail de A. Ghannadiasl et al [22]. L'effet de la 
translation de la plaque sur la surface du milieu élastique, ainsi que le coefficient du milieu 
élastique de type de Winkler et d'autres paramètres sont tenus en compte. Quelques exemples 
numériques sont présentés dans ce travail afin de présenter l'efficacité et la simplicité de la 
fonction de Green utilisée dans la formulation finale. Dans le travail de A. E. Musa et al [23] 
est présentée l'application de la méthode de Ritz pour déterminer les déflexions d’une plaque 
mince reposant sur la surface d’un milieu élastique non homogène. Une approche est élaborée 
aussi pour dériver des fonctions d'essai de Ritz qui satisfont les conditions aux limites de ce 
type de structures. Trois problèmes sont résolus numériquement pour vérifier l'exactitude de la 
méthode proposée. Dans le travail de S. A. Eftekhari et al [24] est présenté une application 
combinée de la méthode de Ritz et de la méthode de quadrature différentielle (DQ) pour étudier 
les vibrations des plaques rectangulaires. Dans cette étude, les dérivées partielles spatiales sont 
discrétisées à l'aide de la méthode de Ritz et ensuite le système résultant des équations aux 
dérivées partielles et des conditions aux limites associées est discrétisé sous forme forte en 
utilisant la méthode DQ. Des comparaisons sont faites avec des solutions analytiques et 
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numériques existantes dans la littérature et les résultats numériques prouvent que cette méthode 
convient très bien au problème considéré en raison de sa simplicité et de sa grande précision. 
L’étude apportée dans le travail de Q. Wang et al [25] présente une procédure de solution unifiée 
basée sur la théorie de la déformation par cisaillement du premier ordre pour l'analyse des 
vibrations libres de plaques rectangulaires orthotropes moyennement épaisses reposant sur la 
surface d’un milieu élastique. Tous les coefficients de dilatation inconnus sont traités comme 
coordonnées généralisées et sont déterminés en utilisant la méthode de Raleigh-Ritz. 
L'excellente précision et la fiabilité de la solution trouvée sont démontrées par des exemples 
numériques et des comparaisons avec les résultats disponibles dans la littérature. L’article de 
M. Dehghany et al [26] présente la solution exacte pour l'analyse des vibrations libres de 
plaques rectangulaires reposant sur la surface d’un milieu élastique. Le milieu élastique est 
décrit par le modèle Pasternak (à deux paramètres). Les résultats numériques obtenus sont 
comparés avec des résultats disponibles dans la littérature et aussi les paramètres importants 
concernant les théories des plaques rectangulaires reposant sur les surfaces des milieux 
élastiques sont discutés. 


Le comportement vibratoire de grandes amplitudes des vibrations libres d'une plaque carrée 
isotrope mince reposant sur la surface d’un milieu élastique de type Winkler est donné dans le 
travail de J. B. Gunda [27]. Ce travail est exprimé sous la forme de plusieurs solutions simples 
en utilisant la méthode de Rayleigh-Ritz (R-R) basée sur les champs des déplacements couplés 
(CDE). La précision et la robustesse des solutions proposées sont comparées avec les résultats 
de base de la méthode des éléments finis disponibles. Une analyse des vibrations libres des 
plaques orthotropes non homogènes reposant sur la surface d’un milieu élastique de type 
Pasternak est investie dans le travail de A. Rahbar-Ranji et al [28]. L'exactitude et l'applicabilité 
de la méthode sont examinées par comparaison des résultats pour différentes conditions aux 
limites et différents types de matériaux avec ceux disponibles dans la littérature. Ils ont conclu 
que la fréquence propre des plaques est plus sensible au coefficient de couche de cisaillement 
plutôt qu'au coefficient de Winkler. Il est investigué dans le travail de V. C. Haciyev et al [29] 
les effets de l'inhomogénéité d’une plaque orthotrope reposant sur la surface d’un milieu 
élastique sur les valeurs de ses fréquences propres. L'équation de base est dérivée en appliquant 
la théorie de Donnell-Mushtari et elle est résolue en utilisant la méthode de Galerkin. Les effets 
de l'inhomogénéité de la plaque et du milieu élastique de type de Winkler sur les fréquences 
propres sont étudiés avec détail. 


L'objectif principal de l’étude apportée dans le travail de A. I. Karakas et al [30] est d'appliquer 
le modèle de fondation Vlasov modifié pour analyser les vibrations libres des tours de 
refroidissement de forme hyperbolique reposant sur la surface de milieux élastiques. Les 
paramètres des fréquences propres circonférentiels, latéraux et extensifs des tours de 
refroidissement sont présentés sous forme de tableaux et graphiques pour évaluer les effets des 
paramètres géométriques et de l'interaction sol-structure. Des exemples numériques sont 
considérés, où il est conclu que le comportement interactif sol-structure conduit à une 
diminution remarquable du paramètre des fréquences du système par rapport à la condition de 
base fixe. Dans l’article de B. F. Apostol [31], il est apporté une solution générale formelle des 
vibrations générées dans un semi-infini élastique homogène et isotrope de modèle de Lamb 
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causées par des forces concentrées. La méthode utilisée est basée sur les transformées de 
Fourier temporelles, les transformées de Fourier spatiales par rapport aux coordonnées 
parallèles à la surface et les équations des ondes généralisées incluant les valeurs de surface des 
fonctions et leurs dérivées. Des résultats explicites sont donnés pour le cas d’une force 
surfacique perpendiculaire à la surface du semi-infini élastique en un point de la surface. Enfin, 
dans le travail de M. Emami et al. [32], un tableau résumant quelques articles relatifs à des 
solutions des différents problèmes de contact liés au problème de Lamb est donné. 


L3. Généralités sur les plaques : 


Les plaques minces ou épaisses de n’importe quelle forme géométrique font partie des éléments 
structurels les plus utilisés dans plusieurs domaines tels que : les structures aéronautiques et 
automobiles, d’ingénierie navale et civile, etc. De ce fait, leurs importances et leurs rôles ne 
cessent de grandir avec le développement technologique. Pour de meilleures performances de 
conception, de résistance, de légèreté etc., il est nécessaire d’investir d’avantage dans le 
domaine des études approfondies des plaques afin d’aboutir à une meilleure structure soit sur 
le plan de sécurité soit sur le plan économique ou autre. 


L3.1. Définition des plaques : 


La plaque est une structure tridimensionnelle solide ou flexible limitée par deux plans parallèles 
(les faces de la plaque) et par un bord perpendiculaire aux faces (l’épaisseur de la plaque). La 
surface moyenne située à m1-distance entre le plan inférieur et le plan supérieur de la plaque, 
est connue sous le nom de la surface moyenne ou "feuillet moyen". L’épaisseur de la plaque est 
définie comme étant la distance entre ses faces. L’épaisseur de la plaque notée généralement 
par le symbole h est supposée petite devant ses autres dimensions qui représentent sa longueur 
et sa largeur. Les plaques peuvent être classées en trois groupes [33] : les plaques minces à 
petites déformations, les plaques minces à grandes déformations et les plaques épaisses à 
déformations négligeables. Le critère souvent utilisé pour définir une plaque mince est le 
rapport entre l'épaisseur et la plus petite dimension transversale. Pour qu'une plaque soit 
considérée comme mince, ce rapport ne doit pas dépasser 7/20 


L3.2. Les déflexions des plaques : 


Les déflexions des plaques se réfèrent aux déplacements de leur plan moyen 
perpendiculairement à leurs plans inférieurs et supérieurs causés par l'action de forces et de 
moments extérieurs [33]. La quantité de ces déflexions peut être déterminée en résolvant les 
équations différentielles appropriées à la plaque étudiée. Les contraintes dans la plaque peuvent 
être calculées à partir de ces déflexions et une fois les contraintes sont connues, les théories de 
rupture peuvent être utilisées pour déterminer si une plaque en mesure de supporter les charges 
agissant sur elle. 
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(a) (b) 


Fig. L 1 : Schéma d'une plaque mince avant et après déformation : 


(a) avant déformation  (b) après déformation 
L.3.3. Comportement général des plaques : 


Les hypothèses fondamentales (ou hypothèses de Kirchhoff) [33] de la théorie de flexion des 
plaques minces avec petites déformations sont : plaques isotropes, plaques homogènes et 
plaques élastiques et qui ce sont toutes basées sur la géométrie des déformations. Ces 
hypothèses peuvent s'énoncer en détail comme suit : 


1- La flèche de la surface médiane w est très faible par rapport à l'épaisseur de la plaque 
h. La pente de la surface fléchie est par conséquent très petite et le carré de la pente 
est une quantité négligeable ; 


2- Le plan moyen reste indéformée latéralement pendant la flexion, c’est-à-dire : 
u(x, y,0)=v(x, y,0) =0 ; 


3- Une section de la plaque normale à son plan moyen reste plane après déformation. 
Cela veut dire que les déformations du cisaillement vertical £. et €. sont 


négligeables ; 


4- La contrainte normale ©, sur tout plan parallèle au plan moyen peut être négligée. 


L.3.4. Champs des déformations (Relations cinématiques) : 


Considérons maintenant l'état de contrainte dans une plaque avec une petite flexion arbitraire 
w(x, y), figure [.2. 
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(a) (b) 


Fig. L 2 : Schéma de déformation d'une plaque mince [33] 


(a) avant déformation  (b) après déformation 


Le plan médian se trouvant équitablement entre le plan inférieur et le plan supérieur de la plaque 
est un plan neutre et donc nous nous concentrerons sur l'état de déformation d’un plan situant à 
une distance z par rapport à ce plan médian. Les pentes dans un plan (x, y) situant à une distance 
z du plan médian ce sont donc : ra) et 2) les déplacements et v de ce plan (x, y) sont 
donnés par [33], [34] : 


= 
ox 
Li 
0 (LD) 
y 
Donc les déformations dans le plan (x, y) sont données par 
LU ON. 
F OL 
Ov O°w 
==: (12) 
©y ©y 
Ou Ov O°w 
” y x Ôxy 
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Un état de contrainte existant dans le plan (x, y) est considéré de sorte que les déformations (€, 


€,» E;y) SONT liées aux contraintes (©, O,,0,,) par les relations suivantes [33], [34] : 


Il 
£. ——(O —vOo),; 
x FA x 


1 
; Pre vo.) (L3) 
1 2(1+v) 
_ - 2. 


Ex 6 XY E xy ? 


Les équations (1.2) et (1.3) peuvent être combinées pour donner (1.4) 


2 2 
6, = . (Tr +) (LA) 
| —v x 
Ez ©w 
T,=-— : 
à 1+v Ox0y 


Ces contraintes varient linéairement par rapport à la distance z et sont équivalentes à des 
moments par unité de longueur agissant sur un élément de la plaque, comme le montre la figure 
IL.1. (b). Donc ; 


h 

2 2 2 
M,= [20,4 =- ci SsÈ ; 

: ox Oy 

2 

: 

2 2 2 
M,= [rod =-D| 4 © | (LS) 

ur ©y ox 

72 

h 

2 2 
M, = [25 dz =-D(1-v) _ ; 


” x Ôy 


DIE 


Ici le paramètre D est défini comme étant la rigidité cylindrique de la plaque et donnée par : 


3 
p=-E 
12(1-v ?) 


Avec E,v: module d’élasticité et coefficient de Poisson de la plaque respectivement. 


= —————— 
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Les composantes de contraintes peuvent s'écrire en fonction des moments de la manière 
suivante : 


I, et 1, : moments d’inertie de la section de la plaque par rapport aux axes x et y 


respectivement ; 


T, et 7, : contraintes de cisaillement. Généralement ces contraintes de cisaillement sont 


négligeables devant les contraintes normales ©, et ©. 


L3.5. Equations d'équilibre d'une plaque : 


Les composantes de différentes contraintes varient généralement d'un point de la section de la 
plaque chargée à un autre point de la même section. Ces variations sont gouvernées par les 
équations d'équilibre statique de la plaque. Considérons un élément de la plaque uniformément 
chargée de dimensions dx et dy, figure L.2 : 

Puisque l’élément est en équilibre statique, alors la somme des forces dans la direction z est 
nulle ce qui permet d’écrire [33], [34] : 


0Q, 
20, dxdy + 9, dxdy + q dxdy = 0; (L.6.a.1) 
Ox ©y 


La formule (1.6.a.1) se simplifie en : 


00. 00, 
. Ye L.6.a.2 
ne += -4; (1.6.a.2) 


De la même manière pour les moments autour de l’axe x on trouve : 


M. OM, 
— dxdy +—— dxdy —Q, dxdy =0; (L.6.a.3) 
x ©y 


La formule (1.6.a.3) se simplifie en : 


om. 2m, _,  . 
ir 5. 6.a. 
Ôx Ôy d 


De la même manière pour les moments autour de l’axe y on trouve : 


CM, ©M 
© + 


= oi (L6.a.5) 
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Finalement, nous introduisons les expressions des efforts tranchants Q, et Q, données par les 


formules (1L.6.a.4) et (1.6.a.5) dans la formule (1.6.a.2), nous obtenons : 


+ 
ox” ©x0y 0y” 


2 © M. OM 
CELP à Le) = _q (L.7) 


Cette expression représente l'équation différentielle de l'équilibre statique des plaques minces. 
Les expressions des forces de cisaillement peuvent alors s'écrire en fonction de la déflexion w, 
en se basant sur les équations (1.6.a.4) et (1.6.a.5) : 


Q.=-D Ô ao) D @) (vw): 


ox\ &x* y ox 
o(®w. 0 Ô ju 
Q,=-D + |=-D(V°w); 
‘ Oy| Oy” x Oy 
Avec 
2 2 
v? — 07 + 07 
x? Oy° 
Rapportons l’équation (1.8) dans l’équation (1.6.a.2) on obtient l’équation (1.9) 
VV2w= À; (L.9) 
D 


Cette expression représente l'équation différentielle des déflexions de la plaque rectangulaire 
trouvée par la savante française Sophie Germain en (1812). 


L4. Modèles des milieux élastiques : 


Quand on parle des problèmes de contact, c’est-à-dire interaction entre structure et un milieu 
élastique, il faut distinguer deux parties différentes, la structure et le milieu élastique. 

Dans les différentes études relatives au domaine des problèmes de contact, la structure est 
toujours modélisée séparément, alors que le milieu élastique est modélisé aussi séparément par 
différents modèles, dont le modèle le plus simple est celui de Winkler. Ensuite, l’analyse finale 
de la structure étudiée se fait par le couplage des deux différentes modélisations faites 
séparément. 


1.4.1. Modèle de Winkler : 


Ce modèle est proposé par le savant allemand Winkler en (1867). 


Le modèle de Winkler repose sur les trois hypothèses suivantes [35] : 


= ————— 
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1. La charge appliquée à la surface du milieu élastique produit des déplacements de la surface 
du milieu uniquement sous la charge appliquée et ne produit pas de déplacements ni de 
contraintes en dehors de la zone chargée. 


2. Le milieu élastique peut résister aux contraintes de la compression, ainsi qu'aux contraintes 
de la traction. 


3. La forme et la taille de la structures en interaction avec le milieu élastique n'affectent pas 
le déplacement de la surface du milieu élastique. 


Ces hypothèses ne sont pas toujours exactes et ne reflètent pas la réalité, car il est bien remarqué 
que dans de nombreux cas une charge appliquée, par exemple au sol, produit des déplacements 
non seulement sous la charge appliquée, mais également aux alentours de la zone chargée. En 
plus, les déplacements de la surface du milieu élastique ne dépendent pas seulement du type du 
milieu élastique, mais aussi des caractéristiques géométriques et mécanique de la structure. 


Le principe de ce modèle de Winkler prévoit que la relation entre la pression exercée sur la 
surface du milieu élastique et son déplacement résultant est exprimée par la formule [35] : 


W, (x, y)=P(x >)/Kk, (L.10) 
Avec : 
W, (x ») : déplacement de la surface du milieu élastique (m) au point d'application de la charge 


extérieure ; 


P(x ») : charge extérieure (N), appliquée sur la surface du milieu élastique ; 


k, : coefficient de proportionnalité qui dépend uniquement des caractéristiques physiques du 


N ; 
milieu élastique | Ve ) ; 


QI? 


Fig. L 3 : Illustration du modèle de Winkler sur une plaque mince reposant sur la surface d’un milieu 
élastique 


q(x ») : réaction du milieu élastique ne +): due à l’action de la charge extérieure PÜE ») ; 


= —— " ————— 
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L'équation différentielle des flèches verticales w(x, ») d'une plaque mince de rigidité 


cylindrique D et soumise à une charge verticale uniformément répartie sur toute la plaque 
P(Xx, y) est donnée par [35] : 


DV'w(x, ») = q(x, »)— PC, y) Œ11) 


En combinant (1.10) et (1.11), on obtient l’équation différentielle du quatrième ordre des 
déflexions des plaques reposant sur la surface du milieu élastique de type Winkler sous la forme 
suivante : 


DV'w(x, y) +kw, (x y) = q 0, >) (L.12) 


Certaines études dans le domaine de l'interaction sol-structure ont été menées sur la base de 
l'hypothèse de Winkler pour sa simplicité. Le problème fondamental de l'utilisation de ce 
modèle est de déterminer la rigidité des ressorts élastiques utilisés pour remplacer le sol sous 
les fondations. Le problème se double puisque la valeur numérique du coefficient de réaction 
de la structure étudiée dépend non seulement de la nature de cette dernière, mais également des 
dimensions de la zone chargée. 


En effet, les chercheurs ont commencé à insérer plus de paramètres et d'éléments dans le modèle 
de Winkler afin de présenter un comportement de la structure plus précis et plus proche de la 
réalité. 


L.4.2. Modèle de milieu élastique à deux paramètres : 


Comme le milieu élastique peut présenter une action d'interaction considérable entre ses 
différentes zones, les ressorts considérés par le modèle de base de Winkler ne tiennent pas 
compte cette interaction entre différents ressorts. Pour surmonter ce problème, des modèles à 
deux paramètres permettant une interaction entre les ressorts ont été proposés. On peut citer 
parmi eux, les modèles suivants : 


1.4.2.1.Modèle de Filonenko Borodich : 


Pour obtenir un certain degré d'interaction entre les éléments des ressorts, les extrémités 
supérieures des ressorts sont reliées à une membrane élastique étirée soumise à un champ de 
tension constante T comme le montre la figure [.4. [36] : 


2 ————— 
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Membrane élastique 


P(x.t) 
TL 7 


k, 


T0 
Fig. L 4 : Modèle de Filonenko Borodich 


La condition d'équilibre dans la direction verticale d'un élément de membrane donne la relation 
de déplacement de la charge [36] : 


P(xt)=kn,(x y)-TV°w, (xt) (1.13) 


Où : 
2 L CCZ Ê 2 ps 
V” : est l'opérateur différentiel de Laplace en coordonnées cartésiennes ; 


T :est la force de traction de la membrane. 
1.4.2.2. Modèle de Pasternak : 


Parmi les modèles du milieu élastique à deux paramètres les plus connus dans la théorie 
d’élasticité on trouve le modèle de Pasternak [37]. Dans ce modèle, Pasternak suppose 
l'existence d'interaction de cisaillement entre les éléments des ressorts. Ceci peut être accompli 
en reliant les extrémités des ressorts à une poutre ou une plaque constituée des éléments 
verticaux incompressibles qui se déforment uniquement par cisaillement transversal. 

La relation entre la force extérieure et le déplacement résultant pour ce modèle est donnée par 


[37] : 
P(xy)=kw (x y)-G,V'w, (x, y) (I. 14) 


Où : 

G,, : est le paramètre d'interaction physique. 

Ce paramètre représente l'interaction due à l'action de cisaillement entre les ressorts. On peut 
facilement remarquer que T de l’équation (1.13) est remplacé par G, dans l’équation (1.14). 
L'équation différentielle d'une plaque mince d'épaisseur constante soumise à une charge 
verticale uniformément répartie sur toute ou sur une partie de la plaque res ») est : 


ee ») (1.15) 


5 ———— 


G, k 
V'w(x, y)— D VW (x, y)+ rs w, (x, y) = 
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L.4.2.3. Modèle de Hetenyi : 


Ce modèle prévoit que la relation entre la charge extérieure P ( ») et la déflexion de la plaque 


W, (x ») est donnée par la formule suivante [38] : 
P(xr)=kwm(x y)-DV'w (x, y) (L16) 


On constate bien la différence entre ce modèle et le modèle de Pasternak réside dans la rigidité 
cylindrique D et le paramètre d'interaction physique G,. 


De plus, d'autres modèles ont été proposés par Reissner (1958), Kerr (1964) et Vlasov & 
Leentiev (1966) dans le cadre des modèles à deux paramètres et d’autres modèles que l’on 
trouve dans le livre de TSUDIK. E (2012) [35]. 


Comme conclusion concernant cette famille de modèles, c’est que l'inclusion du deuxième 
paramètre peut améliorer la précision de la modélisation. En contrepartie, les modèles à deux 
paramètres aboutissent à des formules mathématiques compliquées ce qui rend la détermination 
des paramètres du milieu élastique plus difficile. Enfin, cette famille de modèles basant sur le 
principe des ressorts reste encore loin de la réalité par rapport aux autres modèles. 


1.4.3. Modèle du semi-infini élastique ou modèle de Boussinesq : 


Les charges appliquées sur la surface d’un milieu élastique produisent généralement des 
déplacements non seulement sous la charge, mais également en dehors de la zone chargée, 
figure [L.5. Le semi-infini élastique, comme le cas du sol, se comporte comme un matériau 
élastique suivant les règles de la théorie de l'élasticité. On comprend généralement, qu'un semi- 
infini élastique est un matériau élastique homogène infiniment grand dans les deux directions 
horizontales (x, y) et aussi dans la direction verticale (z) [35] (0 < z<o0,—00( x {00,—00( y (oo). 


Fig. L 5 : Déformation de la surface du semi-infini élastique 


a) selon le modèle de Winkler ; b) selon le modèle de Boussinesq. 


———————— 
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Soit une charge verticale concentrée P est appliquée à la surface du semi-infini élastique, 
comme le montre la figure L.6. 


L 


Fig. L 6 : déformation de la surface du semi-infini élastique : 


1-sous l’effet d’une charge concentrée ; 2- sous l’effet d’une charge uniformément repartie 


Le déplacement vertical de tout point j appartenant à la surface du semi-infini élastique et situé 
à une distance r du point d’application de la charge est donné par la formule de base de 
Boussinesq suivante [39] (proposée en 1885) : 


Le Ft) (L.17) 
rE,r 


Où 
E, : module d'élasticité du semi-infini élastique ; 
v, : coefficient de Poisson du semi-infini élastique ; 


r : distance entre le point où la charge est appliquée et le point où le déplacement se détermine. 


On constate bien que mathématiquement, le déplacement donné par (1.17) au point de 
l’application de la charge P, lorsque r =0, est égal à l’infini. Cependant, en réalité les 
déplacements sous mêmes les charges concentrées sont limités. En plus, généralement les 
charges appliquées sont toujours considérées appliquées d’une manière uniforme sur une 
certaine zone de la surface du semi-infini élastique, comme par exemple pour le cas du sol. 


Dans ce cas, la formule (1.17) peut être réaménagée sous la forme suivante : 


(1-v5) P [I dédn (1.18) 
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Avec : 


b et c dimensions de la zone chargée ; 
x et y coordonnées du point où se détermine le déplacement ; 


& et 7 coordonnées du centre de la zone chargée ; 


Le modèle de Boussinesq, présenté par la formule (1.18), concerne un semi-infini élastique 
linéaire isotrope homogène sollicité par une charge uniformément répartie sur une zone de sa 


surface. Ce modèle représente la base des autres modèles à savoir : semi-infini élastique 
anisotrope, semi-infini élastique non homogène et le semi-infini élastique en couches etc. 


L.4.3.1.Modèle de semi-infini élastique non homogène : 


La non-homogénéité des matériaux est un phénomène courant et très observé dans de nombreux 
matériaux naturels. Comme exemple, le sol est le milieu élastique ayant de fortes non- 
homogénéités en raison du processus naturel de sa formation. Le phénomène de non- 
homogénéité du matériau présent en géomécanique a des applications dans de nombreux 
secteurs technologiques d'importance comme les constructions du génie mécanique ou du génie 
civil. De nombreux efforts ont été fournis pour tenter de comprendre l'influence de la non- 
homogénéité du semi-infini élastique sur son comportement suite à son exposition aux diverses 
charges extérieures. 

Le déplacement vertical de la surface du semi-infini élastique non homogène est donné par 
[40] : 


w(x, ») = ——— (19) 


Avec : 


P(x, y) : charge extérieure appliquée sur la surface du semi-infini élastique ; 


+] 
or | | - 
D È eu inf À) . Eg ; B= 
a 


2 Ja+1 


de) 
24* 2 s) 
He 0 T(3+@) 


F : la fonction gamma ; 


Le module de Young du milieu non homogène peut être donné par [40] : 


E(z) = M2" (O(a{1) (1.20) 
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Où : 
m, : constante égale au module de Young lorsque z =]; 


a : indice de non-homogénéité. Quand & = 1 , le module de Young augmente linéairement avec 
la profondeur et un tel semi-infini élastique est appelé un semi-infini élastique de Gibson. 


Il convient de noter que lorsque & =1, le semi-infini élastique devient homogène et dans ce cas 
2 

on trouve que H=1,F,,=—. Lorsque & =0, on trouve que v, =0,5, B=0 et F,, =— 
T 

L.4.3.2.Modèle du semi-infini élastique transversalement anisotrope : 


Selon Nayak [41], le déplacement de la surface du semi-infini élastique transversalement 
anisotrope w(x, ») dû à une charge ponctuelle P est donné par : 


ONE Jr (121) 


r = 4]x° + y” : distance entre le point où se détermine le déplacement et le point de l’application 
Y P P P PP 
de la charge ; 


(rm, +ms5)(1-as?)-(m, +ms$)(1-as;) 


1, = (1.22) 
2n(1-°)]s.(c-ds?)(1-as2)-s,(c-ds>)(1-as?)] 
m=-cn-(1+b)y,; m,=dn+2au, ; 
11)2 112 
(a+b)+[(a+c) -44 P (a+b)-[(a+c) -44 P (2) 
S, — ; S = AZ —— | 
2d 24 n— 
,-U-28)4,-B(tn)a.+nu, . _- BAT)4. un, + B(+n) - d= 1} 
nu} | ni | n— Li 


L.4.3.3.Modèle du semi-infini élastique en couches : 


Le modèle du semi-infini élastique en couches horizontales ayant une large utilisation dans les 
problèmes de contact. Ce modèle permet d’avoir la réponse du semi-infini élastique à de 
différents niveaux de sa profondeur comme illustré sur la figure L.7. 


= ———— 
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Fig. L 7 : Modèle du semi-infini élastique en couches [42] 


Solon Strukelj et al. [42], l’expression des déplacements verticaux de la surface du semi-infini 
élastique en couches w, (a,0) est donnée par : 


P(œ)1- 
w (a, 0) = 2 = 1 — . (23) 


W, (a,0) : déplacements de la surface de la première couche du semi-infini dus à la force 


ponctuelle verticale P ; 


: . | ; @r L, 
a : coordonnée sans dimension et donnée par : a=— avec c;;, = |— ; 


, 
Cr i 


2, : masse volumique de la couche i du semi-infini élastique ; 


v, : coefficient de Poisson de la première couche (couche superficielle) ; 


“ 


: distance entre le point où se détermine le déplacement et le point où la charge est appliquée ; 


: rapport du module de cisaillement ; 


8 *& 


: fréquence d’excitation ; 


v : rapport du coefficient de Poisson entre la dernière couche et le semi-infini élastique. 
1,= [B, (7) J,(an)dn (L.24) 
0 


B,(r)=n VA (aa aa (00 


JA 
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J: (a n) fonction de Bessel d'ordre 0. 
n : coordonnée radiale sans dimension ; 


|A] déterminant de la matrice [A] définissant les paramètres du modèle utilisé. 


1 Boy kKn u. 


Fig. L 8 : Propriétés mécaniques et géométriques de la première couche du semi infini élastique [42] 


Plus de détails concernant la définition des autres paramètres utilisés par ce modèle du semi- 
infini élastique en couches pour la détermination des déplacements verticaux de la surface 
supérieure dus à l’application d’une force ponctuelle verticale sont donnés dans la référence 


[42]. 


L.4.3.4.Modèle du semi-infini élastique avec une cavité cylindrique : 


Le modèle du semi-infini élastique contenant une cavité doit être présenté d’une manière 
différente. De ce fait, de nombreuses études ont fait l’objet de ce cas de semi-infini élastique 
que l’on trouve principalement dans certains secteurs comme le cas de : pipelines souterrains, 
tunnels et les pipelines sous-marins etc. La réponse dynamique du modèle de semi-infini 
élastique homogène et isotrope avec une cavité cylindrique avec section circulaire à une 
profondeur donnée est analysée dans le travail [43]. La cavité est supposée infiniment longue 
et placée horizontalement de façon parallèle à la surface et soumise à une pression harmonique 


uniformément répartie sur sa surface intérieure en contact avec le semi-infini élastique. 
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Fig. L 9 : semi-infini élastique avec une cavité cylindrique [43] 


La charge harmonique appliquée sur la paroi intérieure de la cavité est donnée par [43] : 
P(?) De : 
La réponse dynamique de la surface libre du semi-infini élastique est donnée par : 


u, =u,CoSÜ—u, Sin 0 (1.26) 

Avec 

u, = ÿ|4 É ACOELAUR un) nr: É Y (&r)-KkY.. un) LC TI NID Ga eos n0 
hn=Ù r r r r 


u,= ÿ|-4. LL &n-B,? Y(&r-C, É ACIDE AU &n) D, É Y. (or) RY.. an) sinnO 
r r r r 


J,et Ÿ, : fonctions de Bessel du premier et du second type et d'ordre n ; 


À,, B,, C, et D, : constantes arbitraires à déterminer à partir des conditions aux limites ; 


k =æo/p/(2+ 2u) et k.=œ\plu : paramètres définissant les ondes longitudinales et 
transversales respectivement ; 


© : fréquence circulaire ; 
À et 1 : constantes de Lame ; 


p : masse volumique du semi-infini élastique. 


L.4.3.5.Modèle du semi-infini élastique homogène à deux tunnels : 


Ce modèle de semi-infini élastique est largement utilisé pour l’analyse des vibrations du sol et 
les bruits induits par le trafic des véhicules et de chemins de fer souterrains. Dans le travail [44], 
il est présenté le modèle du semi-infini élastique homogène à deux tunnels circulaires parallèles 
se trouvant à une profondeur donnée de la surface extérieure. 
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Dans cette étude, le semi-infini élastique homogène avec deux cavités cylindriques est simulé 
comme un solide élastique, isotrope, homogène et tridimensionnel (3D), figure L.10. 


Fig. L. 10 : Modèle du semi-infini élastique à deux tunnels parallèles de sections circulaires [44] 


Les déplacements du sol dans le repère cylindrique local (r,,6,,z) du tunnel 1 sont donnés 
par [44] : 


ET #1 + ns Zaire 


m=-N j=1 
D frire SES T° BE £T2B= |4E 127 
L = [Ts ar > (5 7 i à) ; 7) 
—0 JL CR OMS 
s D SAÛ(n. m,k,,@) 
m=-N 


Les déplacements du sol dans le repère cylindrique local (7, ,6, ,z ) du tunnel 2 peuvent être 
exprimés par [44] : 


= D >} # B° . + > HGuSRei, + 


m=-N j=l n=- 
+ frir es Si | te B° + Ty B; )dk (1.28) 
= = k,, er OR : 
N A 
= SU, m,k, ,@) 
m=-N 


m zm ? m Om ?°°m rm  ? 


> —————— 


ÜrL = a) , Ür£ URL VrL URL R.L RL : 
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ñ,ÿ et w composantes de déplacement de la coque cylindrique dans les directions : z, 0 et 
r respectivement, figure .[.10. 


L.4.3.6 Modèle du quart-infini élastique : 


Le modèle du quart-infini élastique trouve certaines applications dans différents domaines de 
l’ingénierie et de la physique. Certains cas des problèmes de contact font recours à ce type de 
modèle pour étudier certains cas spécifiques où le milieu élastique est effectivement décrit par 
le quart-infini élastique plutôt que le semi-infini élastique. Les études apportées dans les 
références [guenfoud et al [45], Hanson et Keer [46], Aleksandrov et Pozharskogo [47]] font la 
lumière sur ce modèle et ses applications que l’on peut trouver dans le domaine des problèmes 
de contact. Le principe du quart-infini est également utilisé dans la science du cosmos. Par 
exemple, Schultz et Finn-Foley [48] ont décrit les utilisations, les avantages et les inconvénients 
des expériences du quart-infini par rapport aux expériences de semi-infin1 lors de l'examen de 


la dynamique d'impact. 


Fig. L 11 : Modèle du quart-infini élastique [45] 


Dans le domaine des problèmes de contact, l'expression des déplacements de la surface du 
quart-infini élastique due à l'action d'une force concentrée P, est donnée par [45] : 
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P(1-v) Î 
V{x, +B 1—-2v) B 
(x y) 2E a 7 ( ) B, +R 
(1.29) 
2 2 2 
ux|u—Tux+x" +7y 
B,—> B;7 5 ) 
R R 
4 
R={u-x) +y° ; R=J{u+x) +y° LT +4 , y _150091 : B,=1,14597 : 
B, =0,00812 
Où : 


E : module d'élasticité du quart-infini élastique ; 
v : coefficient de Poisson du quart-infini élastique ; 
P : force extérieure appliquée sur la surface du quart-infini élastique ; 


x et y : coordonnées du point M où se détermine le déplacement de la surface du quart-infini 
élastique ; 


u : distance entre l'origine et le point d'application de la force, figure I.11. 
L4.4. Modèle de Lamb : 


Le modèle de Lamb concerne la réponse du milieu élastique due à une excitation extérieure 
verticale ou horizontale appliquée sur la surface ou à l’intérieur d'un milieu élastique remonte 
au début du XXe siècle. L’article relatif à ce modèle est apparu pour la première fois en 1904 
dans la revue scientifique « Philosophical Transactions of the Royal Society of London ». Dans 
ce travail, il a tenu en compte du phénomène ondulatoire des déplacements générés sur la 
surface du semi-infini élastique, isotrope et homogène dus à une charge appliquée sur la surface 
ou à l'intérieur du semi-infini élastique. Le principe de ce travail s’est rapidement propagé et 
appliqué dans plusieurs domaines à savoir les sciences de : tremblement de terre, géo- 
mécanique, fondations, anisotropie générale, problèmes de contact, géodynamique, mécanique 
de contact, élasto-dynamique et autres sciences des problèmes complexes de la physique des 
solides [32]. Le problème de Lamb ensuite est considéré comme référence essentielle pour de 
nombreuses recherches scientifiques vu sa large application dans différente domaine 
d’ingénierie puisque il reflète presque la réalité des problèmes étudiés. 


La fonction de Green définissant les déplacements verticaux de la surface du semi-infini 
élastique de propriétés inertielles, modèle de Lamb, est donnée par [49] : 


v= __ +iruk J(4r) (L.30) 
0 


2x G 


Où : 
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v : déplacements verticaux de la surface du semi-infini élastique de propriétés inertielles ; 
P : amplitude de la charge harmonique ; 
© : fréquence d'excitation ; 


r : distance entre le point où la force dynamique est appliquée et le point où le déplacement de 
la surface du semi-infini élastique est déterminé ; 


x : paramètre ayant la valeur : 1,07236 ; 


2 


26407 -6) 
_ -F()f() 


J,(ér) : fonction de Bessel du premier type ; 


a=\Jé -hR . 


Selon Guenfoud et al. [50], la formule (1.30) devient : 


_kPe° | TK J,(xr 
CC ere | (1.31) 


Où : 
k : paramètre définissant le caractère ondulatoire des déplacements et donné par : 
2 
© X 
L- [exe 
G, 
p : masse volumique du semi- infini élastique ; 


G, 2) avec E, et v, module d'élasticité coefficient de Poisson du semi-infini 
+" 


élastique, respectivement. 


L,;, L, 1, L, : intégrales définies dans la référence [501]. 


Considérons que la charge extérieure Pe“” = Pcos (œt) et en tenant compte que de la partie 


réelle de l'équation (1.31), cette dernière prend la forme suivante [50] : 
(2(r)+2(r)) (1.32) 


Les expressions des intégrales Z,, (r) et Z, (r) sont données avec plus de détail dans le chapitre 
IT. 
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LS. Différents types de structures reposant sur milieux élastique : 


Les différentes structures (1D, 2D ou 3D) utilisées dans ce vaste domaine des problèmes de 
contact trouvent leurs applications dans plusieurs domaines d’ingénierie, de physique, secteur 
de l’aéronautique etc. Le choix du modèle du milieu élastique avec lequel la structure est en 
interaction dépond de plusieurs paramètres à savoir : la complexité de la solution, la couverture 
au maximum de la nature du problème étudié, la précision des résultats recherchés et leur impact 
lors de la mise en service de la structure étudiée etc. 


L6. Méthodes de solution : 


L'analyse statique ou dynamique des poutres, plaques et coques en interaction avec des milieux 
élastiques est basée sur différentes méthodes de résolution. Ces dernières sont issues du 
développement d’anciennes méthodes déjà existantes depuis longtemps et certaines sont 
nouvellement proposées pour répondre aux différentes carences remarquées des méthodes déjà 
existantes, ainsi que pour surmonter certains défis rencontrés lors de leurs applications. En vue 
d’une tendance vers une analyse plus parfaite, plus rationnelle et plus proche de la réalité, le 
besoin du développement de ces méthodes existantes ou de la proposition de nouvelles 
méthodes s’impose. 

En effet, malgré les hypothèses adoptées et les exigences imposées, certaines simplifications 
sont inévitables pour aboutir à la formulation des différentes équations mathématiques pour que 
leurs solutions soient possibles et réalisables. De ce fait, les deux types de résolution connus 
sont bien évidemment la résolution analytique (généralement limitée) et la résolution 
approximative. 


L.6.1. Méthodes analytiques : 


Les méthodes analytiques sont préférables, car elles fournissent des solutions exactes et plus 
proche à la réalité exprimée par le modèle approprié et sélectionné. Néanmoins, les méthodes 
analytiques sont généralement très limitées et difficiles à appliquer pour des structures de 
formes compliquées et constituées de matériaux non-homogène. De ce fait, plusieurs 
hypothèses simplificatrices sont nécessaires. C’est pourquoi, un autre type de méthodes est 
apparu appelées analytiques approximatives ou méthodes semi-analytiques. Ce nouveau type 
de méthodes a réussi à offrir certains avantages par rapport aux méthodes purement analytiques. 
Ces méthodes semi-analytiques généralement présentent la réponse du problème étudié sous 
forme d'une série qui converge vers la solution exacte. Dans cette catégorie de méthodes semi- 
analytiques, on peut citer : méthodes itératives variationnelles (VIM) décrite dans la référence 
[51] ; méthode d'analyse par Homotopie (HAM) utilisée dans l’analyse des structures 
viscoélastiques décrite dans la référence [48] ; méthode de décomposition Adomienne (ADM) 
utilisée pour l’étude des structures immergées dans les fluides décrite dans la référence [49] ; 
techniques de perturbation et la méthode de transformation différentielle sont introduites 
décrites dans les référence [50, 51] respectivement. 
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L.6.2. Méthodes numériques : 


Les différentes méthodes numériques connues sont utilisées pour résoudre les problèmes 
complexes lorsque la solution analytique ou semi-analytique est difficile à exploiter ou 
lorsqu’elle n’existe pas. À cet égard, les différentes méthodes numériques sont employées et 
leurs applications sont inévitables. Les méthodes numériques les plus connues et les plus 
utilisées dans le calcul des problèmes de l’ingénierie et de la physique sont : la méthode des 
différences finies ; la méthode des éléments finis ; la méthode des éléments aux frontières ; la 
méthode des volumes finis. Le détail de ces différentes méthodes numériques utilisées dans le 


domaine des problèmes de contact est donné dans la référence [52]. 
L7. Applications pratiques : 


Dans la section précédente, de nombreuses applications pratiques des problèmes de contact 
comme structures en interaction avec des milieux élastiques ont été présentées. 


L7.1. Applications pour les structures utilisées dans le domaine du génie mécanique ou 
du génie civil : 


Les applications les plus connues des problèmes de contact concernent le domaine du génie 
mécanique et génie civil. Les voies ferrées, les ponts, les chaussées, les réservoirs d’eau et les 
interactions de quelques pièces de machines entre-elles ce sont quelques applications parmi des 
dizaines d’autres que l’on trouve dans ce domaine où la structure est en interaction avec le sol 
ou en interaction avec une autre structure de matériau différent. Plus de détail concernant les 
structures de ce domaine se trouve dans la référence [35]. 


L6.2. Applications aérospatiales : 


Les autres applications des problèmes de contact se trouvent dans le domaine aérospatial où des 
éléments de structures secondaires s’ajoutent pour la suppression des grandes vibrations et leurs 
effets pour prévenir les phénomènes indésirables sur les structures aérospatiales comme le 
mouvement de flottement et le flambage thermique. Ces structures subissent une combinaison 
de forces élastiques, de forces thermiques et des forces aérodynamiques où la nécessité 
d’introduire des éléments de structures légers et rigides en flexion comme le cas des plaques et 
des poutres sandwich ou en matériaux composites. Plus de détail concernant les structures de 
ce domaine se trouve dans la référence [35]. 


L6.3. Applications biomécaniques : 


L'analyse des vibrations du corps humain est importante du point de vue Hygiène, Sécurité et 
Environnement (HSE). II est également très nécessaire pour la conception d'organes artificiels 
pour remplacer certaines parties du corps humain malades d’être représentées par les modèles 
mentionnés précédemment. Par exemple, le crâne humain peut être modélisé par une 
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membrane, les tissus et les vaisseaux sont modélisés par une ficelle et la peau est modélisée par 
une couche mince sur un substrat mou. Plus de détail concernant les structures de ce domaine 
se trouve dans la référence [53]. 


L7. Conclusion : 


Ce chapitre est consacré aux trois points essentiels sur lesquels se base ce travail entrant dans 
le cadre de la préparation de cette thèse de doctorat. Ces points essentiels représentent : (1) l’état 
de l’art de cette thèse, y compris une synthèse bibliographique ; (ii) des généralités et des 
notions de base sur les plaques ; (iii) des généralités et des définitions de base des différents 
modèles des milieux élastiques. Il est indispensable d’effectuer la recherche bibliographique 
pour mettre la lumière sur les différents travaux de recherche similaires et surtout réalisés 
récemment liés à ce domaine pour savoir les différentes réalisations et innovations apportées à 
ce vaste domaine des problèmes de contact. Les définitions de base des différentes structures 
qui peuvent être, lors de leur mise en service, en interaction avec des milieux élastiques de 
différents types sont aussi si importantes, non seulement, pour faciliter leurs études statiques 
ou dynamiques, mais aussi de les concevoir de manière plus rationnelle pour être plus 
sécurisées, plus économiques, plus esthétiques etc. Les définitions de base des différents 
modèles définissant les milieux élastiques sont aussi si importantes pour faciliter 
l’accomplissement des études statiques ou dynamiques des problèmes de contact. Le savoir de 
de ces trois points essentiels apportés avec plus de détail dans ce chapitre est la base de ce travail 
de thèse de doctorat pour faire face aux différents défis de modélisation et de calcul que nous 
définissons dans les chapitres qui suivent ce chapitre. 
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CHAPITRE II : Détermination des déflexions de la plaque 
IL1 Introduction : 


Les plaques de différentes formes, de différentes dimensions et de différents matériaux font 
partie des éléments structurels les plus utilisés dans différents domaines de : construction, 
industrie, aéronautique etc. et leur utilisation ne cesse de grandir avec le développement 
technologique. Cependant, leur étude n’est pas faite complètement à nos jours ce qui nécessite 
l’élaboration de nouvelles méthodes de calcul plus précises ou d’améliorer certaines méthodes 
approchées déjà existantes. Dans ce contexte, on propose une méthode semi analytique basant 
sur la solution générale de Clebsch des équations différentielles des déformées des plaques et 
complétée par l’application du principe de la méthode énergétique de Ritz pour déterminer les 
déflexions d’une plaque mince isotrope de forme rectangulaire et de conditions aux limites : 


libre sur tout son contour et reposant sur la surface d’un milieu élastique. 


Pour déterminer les déflexions de la plaque étudiée, l’équation différentielle des déformées de 
cette dernière est résolue par l’adoption de la solution générale de Clebsch. Cette solution est 
accompagnée de constantes qui se déterminent selon la spécificité de la plaque, à savoir sa 
forme géométrique, ses conditions aux limites, type de son matériau etc. Pour notre cas, on a 
choisi l’application du principe de la méthode énergétique de Ritz pour déterminer ces 
constantes. Ce principe se base sur la combinaison de la variation de l’énergie de déformation 
de la plaque avec le travail de déformation de la plaque lors du passage de cette dernière de 
l’état déformé à l’état initial. Ce principe nous a permis de déterminer les constantes associées 
à la solution de Clebsch et par conséquent avoir la formule mathématique finale représentant la 
solution de l’équation différentielle de la plaque étudiée. La précision des résultats obtenus est 
très satisfaisante et la formule finale est donnée sous une forme très cohérente afin d’être 
compatible lors de son application dans le système général étudié qui est la plaque reposant sur 


la surface d’un milieu élastique (problèmes de contact). 


IL2. Problématique du chapitre : 


Ce chapitre est focalisé sur la détermination des déflexions d'une plaque mince de forme carrée, 
plaque de Kirchhoff, faisant partie du système général étudié, figure. IL.1, qui est plaque 
reposant sur la surface d'un semi-infini élastique (problèmes de contact). 
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Fig. IL.1 : Plaque carrée reposant sur la surface d'un semi-infini élastique 


L'équation différentielle des déflexions d’une plaque en coordonnées cartésiennes est donnée 
par [54], [33], [34] : 


AAW (x, y)= a») (IL.1) 


Avec, 


o © 
A= e + æ : Opérateur de Laplace ; 
x 


En tenant compte de l’expression de l’opérateur de Laplace, l’équation (11.1) prend la forme 

suivante : 

oW (x, oW(xy) 2W(x 

En, MC PE) e05) 2) 
x ©x y Ôy D 


W(x. y) : Fonction définissant les déplacements verticaux de la plaque ; 
q(x y) : Fonction définissant la charge extérieure appliquée sur la plaque ; 


Exh° 


= ————; Rigidité cylindrique de la plaque ; 
12(1-v°) en D 


(x ») : Les coordonnées d’un point de la plaque où se détermine le déplacement ; 
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Pour résoudre l'équation différentielle d'ordre 4, donnée par la formule (IL2), nous adoptons la 
solution générale de Clebsch donnée par la formule suivante [55], [45], [56] : 


co 


W(xy)=W(xy)+X AW, (x) (IL3) 
n=l 

Pour pouvoir étudier ce problème de contact, une plaque rectangulaire reposant sur la surface 
du semi-infini élastique, on propose une approche faisant combinaison de trois différentes 
études. La première étude qui fait l’objet de ce chapitre II est la détermination des déflexions 
de la plaque carrée ; la deuxième étude qui fait l’objet du chapitre III est l’étude de la fonction 
de Green définissant les déplacements verticaux de la surface du semi-infini élastique et la 
troisième étude qui fait l’objet du chapitre IV est la combinaison des deux premières études. La 
plaque étudiée sera remplacée par une plaque de contour libre avec un encastrement fictif en un 
point quelconque de la plaque, fig. IL2. Pour faciliter la présente tâche, on a choisi 
l’emplacement de l’encastrement fictif en un point coïncidant avec son centre. 


Fig. IL. 2 : Plaque carrée de contour libre et encastrée au niveau de son centre 


Par conséquent, les termes de l'expression (1.3) satisfaisant les nouvelles conditions aux limites 
de la plaque sont [57], [56], [58] : 


W(x, ») =W, (x y)+AW, Ce »)+B,W, (Gé y)+A,W, (x, y)+B,W, (x, y)+-. (L.4) 
Avec, 


W, (2, ») : Solution particulière appropriée à la plaque considérée. Elle est donnée par [59], 
[56], [S8] : 
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AT Pab É .) (2 :) ” É .) (2 z) L 
167 D a a b b a a b b 
ee 2 + 2 
(#42) RAT PE MU A Se (IL.5) 
a b b 
1? z° t z? x? ÿ” x y 
EAU ET In rt . 


P : Charge extérieure ponctuelle appliquée en un point de la plaque : 


a et b : demi-longueur et demi-largeur de la plaque respectivement ; 


t et z : coordonnées du point ( j ) où la charge extérieure est appliquée ; 


x et y : coordonnées du point ( É ) où se détermine le déplacement vertical de la plaque causé 


par la charge appliquée au point ( j ) | 


Selon [60], [45], [61] : 


x ÿ 
2xy 
W. , DT 
20 v) ab 
x (IL.6) 
CON ESES 
2 2 
45e) 
Eh A : 
D = ————— : rigidité cylindrique de la plaque : 
12(1-v°) 


E,v : module d'élasticité et coefficient de Poisson du matériau de la plaque ; 


A, B;;, A, et B;; : coefficients à déterminer par l'application du principe de la méthode 
énergétique de Ritz. 


IL3. Application du principe de la méthode énergétique de Ritz : 


Les coefficients À,,, B,,, A3, et B:, associés à la solution de Clebsch sont déterminés par 


l'application du principe de la méthode énergétique de Ritz. Autrement-dit, à partir de la 
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fonctionnelle totale de l’énergie composant de la variation de l'énergie de déformation de la 
plaque et du travail de déformation de la plaque lors de son passage de l'état déformé à l'état 


initial [62], [45], [56]. 


L'énergie de déformation de la plaque U est donnée par [63], [45], [61]: 


D y=b  x=a 
1e, : | w(W(x»))dxdy (1.7) 
Avec, 
 #YŸ Je & (& | 
an 2 de) AE (2) | 
B=2(1-v) ; 


v : coefficient de Poisson du matériau de la plaque. 
Le travail de déformation de la plaque lors de son passage de l'état déformé à l'état initial est 


donné par [45], [64], [65] : 


II=-A,W (#, z) —B,,W, (r, z) — AW, (4, z) —B,W, (r, z)— …. (IL.8) 


Avec, 


W(t,z), W,(,z), W,(,z) et W,(r,z) fonctions définies par les mêmes formules des 
fonctions : W.(x, y), W,(x, y), W, (x, y) et W, (x, y) (expressions IT.6) où les variables (x, y) 


sont remplacés par les variables (£, z) ; 


La fonctionnelle totale de l'énergie de déformation de la plaque est la combinaison des formules 


(IL.7) et (IL8) : 


>=U(x, y)+11(4,z) (IL.9) 


Enfin, les coefficients À,,, B,,, À, et B;, sont déterminés par le système suivant [56], 


[66] : 
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© 3 0. 
PA; | 
© > —_ D: 
©B;; ? 
ss _. (IL.10) 
OA;, ? 
© > —_ D: 
©B;, ? 


La résolution du système matriciel (IL. 10) permet de déterminer les expressions des coefficients 
À. 1. À: B;;, Par conséquent, la formule finale donnant les déflexions de la plaque 
étudiée est déterminée. Mais néanmoins, on veut apporter ici quelques détails concernant les 
étapes de cette résolution pour montrer sa complexité mathématique et les efforts importants 


fournis et à l’aide de Mathematica nous avons surmonté cette difficulté mathématique pour 
arriver enfin à la formule des déflexions de la plaque étudiée sous sa forme la plus simplifiée. 


Pour but de simplification des formules mathématiques finales, on considère que la plaque est 
carrée (a=b) et après une série de différents calculs, transformations et simplifications 
mathématiques nous nous sommes arrivés au système sous la forme matricielle suivante : 


b a P 7 
”_ 0 0 0 -ffotsreaée( Ê | 
a a 
16DB À, cf a : >: (2e) 
oo e ° . B,, | [fers] Loin 
-32D(B-8) À, E b a P+z 
0 0 3 0 B, -ffotsrraae (2) 
-32D(B-8) . 
ù à 0 7 22 t+z 
3a° IDC EI a | 
Où 
: > y f—Xx 
Q = B| 2x2+2x (21 +x)arctan) = |-2ÿ(y-22)arctan| = |+2x(—21+ x)arctan : 
8a 7 y = — 
(84522 }arean( LE] (84 (5-9 Jarn[ 2€ |) 
IX t—x 
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OA -0-2) 2x (ire jaran] à feastirs jaran] 2 frésearian] 2] 
y y 


8a 7 * 
#7 log| x? + y ]-y(y-22)log] x? + y? ]+x(r—iz)log| x + y? |+x(r+iz)log( x + y |+ 


+(r-2) togf 8 2042 + (5-2) ]+ (9 loefr 2m + + (5-2) 


9, 


1 | 6m -3202+16222+ 4078 +828 -4x22p 2( 46 ( 4+ B)+3(x + y°)(-8+38))arctan LA PE 
24a° x y 


12x89) (46 art 16 art | 48ty° arcan| 2 £ £ h96nartan| 27 £)- 
Y—Z Ÿ t—x + 


48 areun| = £ hear pare] = Efur patan] 27 2 |-s6ncpareun| 2 Es arctan] 2 | 
t—x t T—x 


1=X 


2 
—16y° af fra prof 2 [1265 log | x° + y° ]-24x log] x + y |-24y*zlog] x° + » ]- 
>? y 


2, 2. 2 
> RCE 13 |: 
a 


+248 ylog[ 2 2x + + (5-2) ]-85 lo8[ 2042 + (5-2) |- 241 *loe[ —2x+x2+(y-2) J* 


-6y" Blog| x” + »° ]+9x°2Blog] x + y |+95*z8log] x° + DETTE ss 


+4 zlog[ P 2 + x +(5-2) |+245/2108| ? —2n+x + (5-2) |-24ye log] ? — 2x +x° +(»-2) |+ 
+82 log 2x4 +(y-2) |-6Py8lo8[ ° 21 + x +(y-2) [+29 Blog] 2042 +(y-2) |+ 
+37 2Blog[ ° -2x+x +(5-2) | +62 Blog? 2x + x +(5-2) |-928108[ —2x+x2+(y-2) |- 


93 zBlog[ # 242 + (5-2) |+12y2/Flog( 204 x + ( +(y- a) |- 52 Blog[  -2m+ x +(y-2) |- 


2 Dx+x +(y-2) P—Dx+x +(y-2) 
“asste) x _ (> z) frérspts x 7 (y z) } 


a 


1 
Dar 


Q, = (-167 x -8n° +16xyz +16x7° +10P/x8 +27 B—AxyzB +2x7° B+2(3x°7(8-38)+33°2(8-38)+43° (-4+8))x 


ea arn[ 2 |-agreartan] 2£ | 48 zarctan( 2 | 


ÊX t—x Lx 


ST 
+ 
= 


xarctan | 
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IL4. Vérification du nombre de termes pris en compte : 


Pour justifier le nombre des termes de la solution de Clebsch retenus pour notre cas, nous 
procédons à la vérification suivante : 


Premier cas, tenons compte de la solution particulière appropriée à la plaque étudiée et 
uniquement les deux premiers termes de la solution de Clebsch. Deuxième cas, tenons compte 
de la solution particulière et uniquement les quatre premiers termes de la solution de Clebsch. 
Les graphes des déplacements des points situant sur la diagonale de la plaque de caractéristiques 


géométriques et mécaniques suivantes (a=b=1m;v=1/3;E=2x10"N/m ;h=0.01m) 


pour les deux cas sont présentées sur la figure IL.3. 


(i) Lorsque la plaque est sollicitée par une charge, d’intensité g=1KXN /m°, uniformément 
répartie sur toute sa surface ; 


(ii) Lorsque la plaque est sollicitée par une charge concentrée, d’intensité P =14N , appliquée 


au niveau de son centre. 


Il est bien clair et d’après les graphes, les déplacements des points situant sur la diagonale de la 
plaque des deux cas les valeurs des deux cas se convergent bien entre eux. Autrement dit, la 
différence entre les valeurs des deux cas est pratiquement négligeable. Ceci justifie bien le 
nombre de termes de solution pris pour déterminer les déflexions de la plaque étudiée. 


Deux termes de Quatre termes de solution Deux termes de solution 


Quatre termes de solution solution \ 
| ARÈN 


F4 


{ 
j # 
# 
| l ru 


2x 1076 


\ 1.5x 1076 / 


, Lx 106 f 
\ 
À F à 
À 5.x 1077 / 
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Fig. IL 3 : Comparaison des résultats obtenus en tenant compte de différents nombres de termes de la 
solution de Clebsch 
a. Plaque sollicitée par une charge uniformément répartie sur toute sa surface 
b. Plaque sollicitée par une charge concentrée appliquée au niveau de son centre 
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ILS. Exemple d’application : 


Comme exemple d’application, prenons une plaque de caractéristiques géométriques et 
mécaniques suivantes : 


(a=b=1m;v=1/3;E=2x1I0" N/m ;h=0.01m ;). 

Les formes générales des déflexions de la plaque en 3D dues à deux types de charges extérieures 
sont illustrées sur la figure IL.4. Il s’agit des types de chargements suivants : casl, plaque 
sollicitée par une charge, d’intensité q =1N /m°, uniformément répartie sur toute sa surface ; 


cas2, plaque sollicitée par une charge concentrée, d’intensité P =1KN , appliquée au niveau de 
son centre. 


(b) 


Fig. IL 4 : Formes des déflexions de la plaque dues à deux types de charges extérieures : 
a. Plaque sollicitée par une charge uniformément répartie sur toute sa surface ; 
b. Plaque sollicitée par une charge concentrée appliquée au niveau de son centre. 
Les projections des graphes de la figure IL.4 sur le plan Oxy, ce qu’on appelle aussi les schémas 


isovaleurs, sont illustrées sur la figure ILS. 


Fig. IL. 5 : Présentation de la projection des déflexions de la plaque (isovaleurs) sur le plan Oxy dues 


à deux types de charges extérieures : 
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a. Plaque sollicitée par une charge uniformément répartie sur toute sa surface 
b. Plaque sollicitée par une charge concentrée appliquée au niveau de son centre 


IL.6. Résultats et discussion : 


L’illustration des résultats de calcul, concernant les déflexions de la plaque étudiée donnés par 
la formule finale (IL.4) représentant la résolution de son équation différentielle, est faite sur les 
figures IL.4 et IL.S. Il est montré clairement, que les déplacements de la plaque sont en parfaite 
concordance avec les déflexions naturelles de la plaque étudiée sous l’effet de ces types de 
chargement. Pour le cas où la plaque est sollicitée par une charge concentrée appliquée en son 
centre, il est bien remarqué que le déplacement maximal se trouve au centre de la plaque et plus 
on s’éloigne du centre de la plaque, les déplacements se minimisent. Il est aussi bien clair que 
les déplacements deviennent un peu plus grands aux niveaux des coins de la plaque, puisque la 
plaque est libre sur tout son contour. Pour le cas où la plaque est sollicitée par une charge 
uniformément répartie, on remarque bien que les coins de la plaque possèdent les déplacements 
les plus faibles. Ceci est en parfaite concordance avec l’état naturel de ce type de chargement 
pour la plaque étudiée. 


IL7. Conclusion : 


L’étude apportée dans ce chapitre est consacrée à la détermination des déflexions d’une plaque 
mince de forme carrée de n’importe quelles conditions aux limites en partant de l’équation 
différentielle générale des plaques en coordonnées cartésiennes. Parmi les solutions connues 
dans la littérature de cette équation différentielle d’une manière générale, on note celle de 
Clebsch qu’on a adopté. L’introduction des conditions aux limites appropriées à notre cas nous 
a permis de spécifier et simplifier la solution recherchée sous forme de combinaison de 
plusieurs termes respectant les conditions aux limites de la plaque étudiée et qui sont aussi 
accompagnés de coefficients à déterminer (11.4). La détermination de ces coefficients nécessite 
encore l’application du principe de la méthode énergétique de Ritz qui a abouti à un processus 
de calcul mathématique très compliqué. Cette complexité est surmontée grâce aux efforts de 
calcul importants fournis, ainsi qu’à l’aide du logiciel de calcul formel Mathematica et 
finalement la solution est donnée sous une forme simplifiée permettant son application aisément 
dans les calculs ultérieurs. Il est démontré aussi que le nombre de termes de la solution de 
Clebsch pris pour notre cas est largement suffisant, car la différence entre les résultats trouvés 
avec trois termes et ceux trouvés avec cinq termes est au voisinage de zéro. 
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CHAPITRE III : Etude de la fonction de Green définissant les déplacements verticaux de 
la surface du semi-infini élastique de propriétés inertielles 


IIL.1. Introduction : 


Dans ce chapitre on présente une approche semi-analytique pour l’évaluation de certaines 
intégrales de la fonction de Green définissant les déplacements verticaux de la surface du semi- 
infini élastique de propriétés inertielles (modèle de Lamb) sur lequel repose une plaque 
rectangulaire qui fait l’objet d’une étude dynamique. Le grand challenge de cette étude 
dynamique de la plaque est de surmonter le problème de singularité associée à ce problème de 
contact très compliqué. 

L’approche semi analytique proposée pour résoudre ce problème de contact compliqué se base 
sur la discrétisation du système étudié (plaque rectangulaire reposant sur la surface du milieu 
élastique de propriétés inertielles) ce qui conduit à une solution numérique finale sous forme 
matricielle. Bien évidemment, tous les termes de la matrice sont doublement indicés et la 
singularité se présente seulement au niveau des termes ayant le même indice. C’est pourquoi, 
une attention particulière est tournée vers les termes de la matrice de mêmes indices, c’est-à- 
dire les termes diagonaux où la singularité apparaisse. Cette situation nous a imposé à surmonter 
ce problème de singularité par n’importe quel moyen et il s’est avéré que le seul moyen possible 
est l’évaluation analytique ou semi-analytique des intégrales de la fonction de Green associées 
à ce problème. Par contre, les intégrales de la fonction de Green des autres termes de la matrice 
d’indices différents sont évaluées par de méthodes numériques aboutissant à de solutions 


approchées de très bonne précision. 


IIL2. Problématique du chapitre : 


Ce chapitre est focalisé sur l’étude de la fonction de Green. Autrement-dit, la détermination des 
déplacements verticaux de la surface du semi-infinie élastique de propriétés inertielles que l’on 
trouve aussi dans la littérature sous le nome modèle de Lamb. 

L’approche proposée consiste à discrétiser le système étudié, plaque rectangulaire reposant sur 
la surface du semi-infini élastique de propriétés inertielles, en un nombre d’éléments de 
dimensions identiques. Donc l’étude dynamique sera accomplie sur un système discret et par 


conséquent, le point d’application de la force causant les déplacements verticaux de la surface 
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du semi-infini élastique et le point de détermination de l’un de ces déplacements sont considérés 
au niveau des centres des éléments. De ce fait, on se trouve devant deux cas de situations : cas 
(1) lorsque le point d’application de la force coïncide avec le point de détermination du 
déplacement, fig. IIL1 ; cas (2) lorsque le point d’application de la force et le point de 
détermination du déplacement sont différents, fig. IIL.2. La singularité concerne uniquement le 
cas (1) et donc tous les efforts de résolution s’accentuent sur ce cas (1) afin de surmonter ce 
challenge de singularité qui s’est avéré ni possible que par la résolution analytique ou semi- 
analytique de toutes les intégrales de la fonction de Green de ce problème de contact compliqué. 
En revanche, les intégrales de la fonction de Green concernant le cas (2) ne nécessitent pas une 


résolution analytique, juste une évaluation numérique avec une précision très suffisante. 


Fig. IL 1 : Géométrie discrétisée de la zone de contact de la surface du semi-infini élastique 


IL3. Cas (1): 


La fonction de Green définissant les déplacements verticaux de la surface du semi-infini 
élastique de propriétés inertielles (modèle de Lamb) due à l'action de la charge externe est 
donnée par [49], [67] : 


. —kcos(ot) 
| 27rG, 


2 ———— 


(Z(r)+B(r)) (III.1) 
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k : paramètre définissant le caractère ondulatoire des déplacements et donné par : 
2 
© X 
FE 
G, 
& : fréquence d'excitation ; 
t : temps d'excitation ; 
p : densité du semi-infini élastique ; 


E, 


G, : module d’élasticité transversale du semi-infini élastique : G, = ——— 
2(1+v,) 


E, et v, : module d'élasticité et coefficient de Poisson du semi-infini élastique respectivement ; 
r : distance entre le point d’application de la charge et le point où se détermine le déplacement ; 


Les expressions de Z,,(r)et Z,(r) sont données par [50] : 


I — F,(m1m+1;-Kkr/4)-4" E. A (IIL2) 
21 = Om 152 us) ° 1 2 1 "16 : 
I. =d, {(x 22) ,F,(1/2:13/2;-k °° 14)}+ 

+d, {Gr /2){7,(xkr)+H,(xkr)]- 


_ D (tm) T 1F,(A+m)/2:1(8+m)/2;-k Fe. 14) (UIL.3) 


m=0 
+d, 28" | -K ‘r°,F(L12,3,3-4°r7/4)+ 
+16(4-k%r2(2-27+m4)+k rt mfe2r2)) |]. 
Où : 
di, d,, d,, b,, b,,..…., b,, : coefficients déterminés par la méthode des moindres carrés, 


Demidovich et al [68] ; 


H,(@) et Y,(a) : fonction Struve et fonction de Bessel du second type, Gradshteyn et Ryzhik 
[69], 


en (a;b: d) : fonction hypergéométrique généralisée, Gradshteyn et Ryzhik [69] ; 
7 50.577216, constante d’Euler ; 


selon [50], 7 =1.07236, pour v=1/3. 
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Il est important de noter que l'application de l'approche proposée nécessite l’intégration de la 


formule (IL.1) sur toute la surface de l'élément chargé de dimensions : b=2xb2 et c=2xc2 


(figure IIL.1). Donc, la variable r devient r = /x? + y° : 


—kcos(ot) "FF — — 
ee], Le (dre)+ (Ve) Jen (IIL.4) 
0 


y=—b2 x=-c2 


L'expression (IIL.4) est divisée sur la surface de l’élément chargé bc, car le chargement est 


considéré comme étant uniformément réparti sur toute la surface de l’élément chargé. 


Fig. LIL 2 : Géométrie de l'élément chargé montrant la zone sur laquelle l'intégration doit être effectuée 


Dans le cas (1), le point à coïncide avec le point j (Fig. IIL2), alors l'équation (IIL.4) doit être 


intégrée sur la toute la surface de l'élément chargé par une charge considérée uniformément 


répartie. C’est-à-dire : 


et | in -8ÿ +65) }+2 (2) + (5m) Jésan (LS) 


—b2 -c2 


Où : 


r=(x-€) +(5-7) : 


x et y : coordonnées du point i où se détermine le déplacement ; 


ë et 7 : coordonnées du point j où la charge est appliquée ; 
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Certaines intégrales de J,, et Z, de la formule (II.S) sont évaluées sans difficulté, maïs d'autres 


intégrales où apparaissent les fonctions spéciales ne sont pas évaluées en raison de leur 
complexité. Pour faire face à ce challenge, l'approche de Johnson [70] est appliquée, dans 
laquelle un changement de variables du système de coordonnées cartésiennes en un système de 
coordonnées polaires est exécuté. Cette technique permet l’évaluation de toutes les intégrales 
de la formule (IIL.S), et par conséquent la détermination de l’expression finale du déplacement 


y. 


Après avoir introduit l’approche de Johnson, l’expression (IL.S) prend la forme suivante : 


b c 
—kcos(wt) NC Vas D. 


ee ] ] LL, (rp)+1(r.p)] trdrdp (IIL.6) 


IIL3.1 Exemple explicatif : 


Pour plus de détail concernant cette technique de changement de variable utilisée, principe de 


Johnson, pour évaluer les intégrales de la formule (IIL.6) du cas (1), une intégrale de Z,, est 


choisie pour expliquer la manière utilisée pour évaluer ces intégrales compliquées : 


r= Ge) + (5-7) : 


n=b2 E=c2 


[ [ (a)déan- > Be TG F,(m1,m+1-r /4))dédm- a TC (mme 


m=-b2 £=-c2 m=-b2 E=-c2 m=-b2 E=-c2 


2, 
0 


Choisissons cette intégrale, 


(2) c2 
g=arctan| — | r= 
€ 


n=b2 E=c2 R2 >] cos(o) -j2,2 
[ Î (mine) dédn=8x | | (mime) rdrdg 


n=-b2 £=-c2 p=0 r=0 
A l’aide du logiciel de calcul formel Mathematica, on trouve : 


On calcule l'intégrale par rapport à r, on trouve : 


c2 


TT cos(p) 2 2 2 $ 
Î LE fmmer je ef moment | I 
16 p) ( 


16(cos Ep) cos Ep) 


On calcule l'intégrale par rapport à o, on trouve : 


Faisons encore le changement de variable suivant : f = 


, On trouve : 
cos @ 
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k(c2) , 


Faisons approcher la fonction ,F, Ê ARR |: sous forme d’une série à cinq 


termes, on trouve : 


k° ae ni Dem +0[r 
m 


Ce qui permet d’écrire : 


4 
=arctan| LE =arctan b2 272 1 
Ci) | “ Ê 2,m+l; k (2) | l le S (2) | ES l ® . 1 


16(cosp) J(cosp) 2 (cosp) 32(1+m) 


__b2(96—(-96+((b2)° +3(c2) )k° )m)c2 
192(m +1) 


Le résultat final de cette intégrale est donc : 


n=b2 ë=c2 


&((x-8) +(5-7)) d£an 8x 22067 -96+ (2) +3(c2) me? 


El ml m+l; 
16 192(m +1) 


n=-b2 £=-c2 
__b2(96—(-96+((b2)° +3(c2)°)k° )m)c2 
E 24(m +1) 


Enfin, cette procédure d’évaluation de l’intégration est suivie pour le reste des intégrales de la 
formule (IIL.6), ce qui a permis de définir la formule finale du déplacement d’un point de la 


surface du semi-infini élastique lorsque la charge coïncide avec ce même point, le cas (1). 
IHL3.2 Vérification : 


Pour montrer la fiabilité et l’exactitude de la méthode d’évaluation utilisée de ces intégrales, on 


procède à la vérification suivante : 


Prenons une intégrale de la formule (IIL.6) dont son évaluation ordinaire est possible, c’est-à- 


dire on peut l’évaluer classiquement, sans recours à l’approche de Johnson : 
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E (227 +62 -b2* log[ b2° ]-c2*log[ c2° ]+2c2/ log [2 +42 +02 | +2b2° log] c2+ Vb22+02° l 


Si on utilise l’approche, décrite ci-dessus, on trouve : 


n=b2 &=c2 n - p=arcan( 2) re 
(x-6) +(»-7) Jd£an=8x Î Î (r)rdrdg 
m=-b2 £=-c2 L he 


On calcule l'intégrale par rapport à r , et à l’aide du logiciel Mathematica, on trouve : 


c2 


E _. c2 
2 3(cos p) 


On calcule l'intégrale par rapport à o, et à l’aide du logiciel Mathematica, on trouve : 


b2 
p=artan( 2) : 3 œ 
Î _ de D Des +2c2arctanh| tan garant À 
go | 3(coso) 6 c2 2 c2 | 
Le résultat final de cette intégrale est donc : 
n=b2 £=c2 2 
2 2 1 b2 1 b2 
= +(y— dédn =8x| —c2"| b2,[1+—— +2c2arctanh| tan| —arctanh| — 
nel (x é) (» 1) san ie | | > c2arcta un] aa 21) 
2 
— 4» Pa [ + + 22h tan Ê arctanh 21) 
3 c2 2 C2 


En tenant compte des valeurs suivantes : b2=0.1; c2=0.1, les deux formules obtenues par 
l’application de deux différentes méthodes d’évaluation des intégrales donnent le même 
résultat. C’est-à-dire : 


y=0.1 x=0.1 


[ [ ()d= = (4202 b27 +027 -b2" log[ b2° ]-c2" log c2° |+ 


y=—0.1 x=—0.1 
+202 log |b2 +022 +02? | + 2b2° log |e2 + 22 +02? ) = 0.00306078 


y=0.1 x=0.1 2 
( Î (r)dxdy = 4» b2, fi + +2c2arctanh| tan Ê arctanh ee) = 0.00306078 
y=—0.1 x=—0.1 3 c2 2 C2 
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IIL3.3 Précision de la solution obtenue en considérant seulement cinq termes des séries 
d’approximation des fonctions : 


Pour s’assurer de la précision des résultats obtenus en appliquant l’approche proposée pour 

l’évaluation des intégrales, apportons comme exemple la comparaison graphique de la fonction 

hypergéométrique ,F, avec son approximation sous forme d’une série de seulement cinq 

termes. Considérons les caractéristiques mécaniques du semi-infini élastique de propriétés 

inertielles suivantes : 

E, =250x10° N/m’; 

V, = 0.3; 

p=1500 Kg /m’; 

Le Il 

cos(p 

T . 
4’ 


| 


Pix = — VME [151 1]. 


1F2//Series — La fonction approchée 1F2//Series — La fonction approchée 
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—— La fonction exacte PS 


500 1000 1500 2000 n) (Hz) 
(a) (b) 


Fig. IL. 3 : Comparaison de la fonction hypergéométrique ,F, avec sa série d’approximation de 


seulement cinq termes en fonction de : 


(a) fréquences d'excitation © lorsque c2=0.1m 
(b) distance c2 lorsque © = 250 Hz 


La distance c2, définit la moitié de la longueur de l'élément chargé, varie en fonction des 
dimensions de la zone de contact étudiée et aussi en fonction de la taille du maillage de sa 
discrétisation. Pour notre cas, nous avons choisi une zone de contact carrée de côtés de 2 met 
de maillage moyen. Donc la valeur de c2 ne dépassera pas 0.1 m pour un gros maillage. 


D’après les graphes (a) et (b) de la fig. IIL.3, on remarque bien que l’erreur entre la fonction 


PF; et la série de cinq termes de son approximation est pratiquement zéro. Cela signifie que 
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l’approche de résolution proposée basant sur le principe de Johnson aboutit, sans aucun doute, 
à une solution exacte. 


Enfin, et par la même manière on trouve l’évaluation des restes des intégrales de la formule 
(IL.6) ce qui conduit à la formule finale de la fonction de Green définissant les déplacements 
verticaux de la surface du semi-infini élastique de propriétés inertielles (modèle de Lamb), 
lorsque le point de l’application de la force j coïncide avec le point i où se détermine le 
déplacement, le cas(1). 


(este a 


2bcrG, 


(©) c 
arctan| — 
c2) 2cosp 


A=8 Î Î LL, (r,0)+L (r,@) |rdrdg = 


0 0 


ï 1 [(b2c2(24-(-244(b22 +3027)£°}m) b2c2(96-(-96+(52° +3c2°)k°}m) 


= b X 4" : 
DE 6(1+m) 24(1+m) 
16b 2° Pan] an] arctan( 2 |] 
+dÿx| x (FE ma (px 432 )et)) 4 
k b2 18 
1+—— 
c2 


| Da (2+c2k xx) 23-50 1 124. b2 2.2 
“apr (| TE Pe 24(24-12c27k ? 37 +c2*k* 7" )arctan HULL ER 


x|127(-24+ (62° +302 )k 22° )+144(3+l08 [4])-# 222 x(b27 (19 +10g[4096])+ 


+362 (23+l0g[4096])) +6 (-24 +(b2° 2) (PE |a(-as (2 +362 )k °° Jlog[c24 a) 
CG 


F3 2 
ARR pa f+t2 (60 (202 +502 )k°? 4 )+302(-20 +024? 3° )x log| cos parctn( 2 sin rarcan] 2 | 
457 c2 2 c2 2 c2 


| 


6(3+n) 


paf eo[ Lurun(22]] ei) ; _—. 


4 
+(d,)x128 x 2 [-241920(352" +10b2*c2? +15c2*)+1080(5b2° +21b2*c2? +35b2?c2* +350 2° )k°? - 
2 2721600 


—(35b2* +180b2%c2? +378b2*c2* +420b2°c2° +315c2*)k ‘) 16/4 É (b2* +3b2c2°)k° 2-27 + Log[4]) + 


A3 {rscrtaun[22] 12202 7b2? -33c2° -6(b2° ar )uaf(tr ea ))] 
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Avec, selon [50] : 


4 =1.0723562676808107 ; 


d, =—0.638753325488385 ; d, =0.14255586708614507 ; d, =—0.45897252064027794. 


IL.4. Cas (2): 


Dans ce cas le point d’application de la charge j et le point où se détermine le déplacement à 

sont différents, fig. IHIL.1. Ce cas (2) est caractérisé par l'absence totale de la singularité puisque 
la distance r entre les points i et j est différente du zéro, c’est-à-dire r Æ0. Donc la formule 
finale des déplacements de la surface du semi-infini élastique de propriétés inertielles pour ce 
cas (2) est évaluée numériquement et elle est donnée tout simplement par l’expression suivante : 


een DA CET Ga +6-n))] (LT) 


S, n=]l 


Le paramètre s, paru dans l’équation (IIL.7) exprime le nombre de sous-éléments auxquels 


l’élément chargé est divisé. Pour une précision testée très suffisante, on a opté pour un nombre 
de sous-éléments égale à 16 , c’est-à-dire s, =16, (fig. IIL.4). 


Et 
(Il 
AS 


Fig. IL. 4 : Division de l'élément chargé en 16 sous-éléments 


IL.4.1 Vérification de la précision de l’évaluation numérique des intégrales : 


L'évaluation numérique des intégrales de la fonction de Green (eq. IIL.7) concernant le cas (2) 


nécessite impérativement de tester sa précision en comparant les résultats numériques trouvés 


a —— 
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avec les résultats analytiques disponibles. Pour ce faire, on procède au processus de 


comparaison de l’intégrale suivante : 


Es ï ' oaus-[ DT T (x-5) #(5-n) Jasan = 


y=—b2 x=-c2 


(te Meg ren ho em 
x(b2-y)-2(c2-x)\/(c2-x) +(b4- y) (b4-y)+2(c4-x)J(c4-x) +(b4- 5) (b4-»)-(b2- 5) x 
«iog| c2+x+4(c2-x) +(b2- y) |+(b2 »)'te| c4+x+4/(c4-x) +(b2 >) |+(b4 DE (LS) 
«iog| c2+x+4(c2-x) +(b4 ” CA | c4+x+4/(c4-x) +(b4 ” | 
(24) el b2+4/(c2- x) +(b2- y) -» |+(c4 s)'te b2+/(c4-x) +(b2- 5) » + 


(c2 s)'t| ba+f(c2-x) +(b4- y) + » cs) og] ba+(c4-x) +(b4-5) »])} 


| 


(b4) y=7 


Fig. IL 5 : Numérotation des sous-éléments de l'élément chargé 


L’application de l’équation (IIL.7), en choisissant 16 sous-éléments numérotés comme la montre 


la figure IILS, aboutit à l’expression suivante : 
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NT 2 
GB) T | 


(IIL.9) 


Le tableau IL. 1 compare les résultats donnés par les deux formules (IIL.8) et (IIL.9), ainsi que 


l’erreur du calcul entre elles en fonction des valeurs de b et c : 


Valeurs de cet b 


Résultats donnés par 


Résultats donnés par 


Erreur du calcul 


(IL.8) (IL.9) 
c=b=02m 3.3946035676 3.3945728732 3x10% 
c=b=0.1im 1.6973017838 1.6972864366 1.5x10* 
c=b=0.0666667 m 1.1315345225 1.1315242910 1.02x10° 


Tab. IL. 1 : Comparaison des résultats donnés par (LS) et (HIL9) pour différentes valeurs des 


dimensions 


de l’élément chargé b et c 


IIS. Résultats et discussion : 


Pour plus de détails concernant le domaine d’application de la solution trouvée dans le domaine 


des problèmes de contact en général, deux exemples sont présentés. 
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IIL5.1 Variation des déplacements de la surface du semi-infini élastique de propriétés 
inertielles : 


Les formules finales trouvées par cette étude montrent que les déplacements verticaux de la 
surface du milieu élastique de propriétés inertielles causés par une charge ponctuelle sont 
obtenues en fonction de deux paramètres : les fréquences d'excitation « et les dimensions de 
l'élément chargé supposé carré de côté b respectivement. 


v (10m) 


100 200 300 40 5 (Hz) 


(a) (b) 


Fig. IL. 6 : Variation des déplacements verticaux de la surface du semi-infini élastique en fonction de: 
(a) fréquence d'excitation © et dimensions de l'élément carré chargé b (3D) ; 
(b) fréquence d'excitation © ; 
(c) dimensions de l'élément carré chargé b. 


On constate bien que les graphes (a), (b) et (c) de la figure IIL.6 sont en proportionnalité inverse. 
C’est à dire les déplacements de la surface de la zone contact du semi-infini élastique de 
propriétés inertielles tendent vers zéro lorsque les valeurs des fréquences d'excitation et les 
dimensions de l'élément chargé augmentent. 
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IIL5.2 Présentation en 3D des déplacements verticaux de la surface du semi-infini 
élastique dus à deux types de charges extérieures : 


La présentation graphique des solutions trouvées des cas (1) et cas (2) apportée sur la figure 
(IL.7) en 3D illustre la distribution des déplacements verticaux de la surface du semi-infini 
élastique de propriétés inertielles dus à deux types de charges extérieures. 


(a) (b) 


Fig. IL 7 : Distribution des déplacements verticaux de la surface du semi-infini élastique causés par : 
a) une force ponctuelle appliquée sur la surface du semi-infini élastique ; 
b) deux forces ponctuelles appliquées sur la surface du semi-infini élastique. 


On constate bien de la figure (IIL.7) et pour les deux cas de chargement, que la valeur maximale 
du déplacement est obtenue au point où la force ponctuelle est appliquée. En outre, plus on 
s'éloigne du point d'application de la force ponctuelle, les valeurs des déplacements diminuent 


et convergent vers ZÉrO, y — 0]. . Ceci est en parfaite concordance avec la logique et avec la 


réalité des choses dans ce domaine des problèmes de contact. 
IIL6é Comparaison des résultats obtenus : 


Pour valider les résultats obtenus via cette étude de la fonction de Green définissant les 
déplacements verticaux de la surface du milieu élastique de propriétés inertielles (modèle de 
Lamb), alors une étude comparative s’impose. Pour ce faire, on considère un milieu élastique 
de propriétés distributives (modèle de Boussinesq) défini par la formule suivante [71] : 


2 
y, = re F. (IIL.10) 
bcrE, 
n=b2 E=c2 
Le dédr (I.11) 


m=-b£=-02 (x 6 ÿ + (> — n) 


5 —————— 


CHAPITRE IN : Etude de la fonction de Green définissant les 2021 
déplacements verticaux de la surface du semi-infini élastique de propriétés inertielles 


Avec: 
E, : module d'élasticité du semi-infini élastique de propriétés distributives ; 
v,, : coefficient de Poisson du semi-infini élastique de propriétés distributives ; 


b et c longueur et largeur de l'élément chargé respectivement. 
IIL6.1. La surface du semi-infini élastique est sollicitée par une charge statique : 


Considérons un élément carré de dimensions, b2= c2 (Fig. HL.1), chargé entièrement par une 
charge uniformément répartie. L’évaluation de l’intégrale de la formule (HIL.11) définissant 
l’expression du déplacement vertical en un point i de la surface du semi-infini élastique de 


propriétés distributives (modèle de Boussinesq) F;; est: 


F,=(2-phoefheee frame GE | 

een re |. 
-(b2+ log) -b24x+ 2 -2b2x 4 (624) |+ 
+(82-x)log| 02 + + /b2 202048 (6245) |+ 
+(b24»)lo8 5244 02 42244284 (0243) }+ 
+(b2+)log) 52494 b2+2b2x 4 +(b24 3) |+ 
+(-b2+ y)log]-62+2+ 42027 +47 +37 -202(x4 5) |+ 
+(-b2+x)log]-b2+ 3 +202 +2 + 372b2(4+ 5) 


Pour pouvoir comparer les valeurs des déplacements verticaux de la surface du semi-infini 
élastique donnés par les deux modèles (Lamb et Boussinesq), on fixe la valeur de la fréquence 
d'excitation pour le modèle de Lamb à une valeur moyenne © = 250 Hz et on utilise les mêmes 


(III. 12) 


caractéristiques mécaniques du semi-infini élastique pour les deux modèles : 
E,=250x10$ N/m°; v, =0.3; p =1500 Kg /m'; 


Après le calcul des déplacements verticaux de la surface du semi-infini élastique par les 
formules des deux modèles comme il est illustré graphiquement sur la figure IIL&, on constate 
bien que l’allure des déplacements des deux modèles est la même, mais celle de Lamb donne 
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une valeur de déplacement assez grande au point d'application de la force par rapport à celle de 
Boussinesq. 


-- 0.00001 


Fig. IL. 8 : Les épures des déplacements verticaux de la surface du semi-infini élastique données par : 
(D modèle de Lamb 
(ID) modèle de Boussinesq 


IIL6.2. La surface du semi-infini élastique est sollicitée par une charge dynamique : 


Comme indiqué précédemment que les expressions finales des déplacements de la surface du 
semi-infini élastique de modèle de Lamb pour le cas (1) ou pour le cas (2) sont trouvées en 
fonctions de deux variables qui sont la fréquence d’excitation et les dimensions de l’élément 
chargé. Cela indique bien que ce modèle tient compte du caractère ondulatoire des 
déplacements de la surface du milieu élastique sous l’effet d’une charge dynamique ce qui le 
rend le modèle le plus proche de la réalité. Les schémas (a) et (b) de la fig. IIL.9 montrent la 
propagation des déplacements sur la surface du semi-infini élastique de propriétés inertielles 
dans toutes les directions, causés par l'effet de : (a) une charge dynamique ; (b) deux charges 
dynamiques de même intensité. 

Les valeurs des caractéristiques mécaniques du semi-infini élastique et celles des charges 
dynamiques prises sont : 


E,=250x10° N /m°; v, =0.3; p=1500Kg /m°; P = Fcos(œt); F =1000 N; © =250 Hz; 
t=0.195s; 


P : charge dynamique extérieure appliquée sur la surface du semi-infini élastique ; 
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F : amplitude de la charge dynamique appliquée ; © : fréquence d'excitation ; f : temps 
d'excitation. 


Fig. IIL. 9 : Caractère ondulatoire des déplacements de la surface du semi-infini élastique à { =0.19s 
a) causés par une charge dynamique ; 
b) causés par deux charges dynamiques. 


Grace au logiciel du calcul formel Mathematica, on peut même présenter la projection 
instantanée du caractère ondulatoire des déplacements de la surface du semi-infini élastique de 
modèle de Lamb causés par une charge dynamique et par deux charges dynamiques. Cette 
présentation est 1llustrée sur la figure IL. 10. 


(a) 6) 


Fig. IIL 10 : Présentation de la projection des déplacements verticaux de la surface du semi-infini 


élastique (isovaleurs) sur le plan Oxy causés par : 


a) une charge dynamique 
b) deux charges dynamiques 
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Enfin, la figure IIL.11 montre une superposition des formes et des amplitudes des déplacements 
de la surface du semi-infini élastique données par le modèle de Lamb et par le modèle de 


Boussinesq causées par une charge dynamique extérieure P=F cos(œt), Où : 
E, =250x10° N/m°; v, =0.3 p =1500Kg/m°; F =1000 N; © =250 Hz. 


En comparant l’allure des déplacements donnés par les deux modèles, on constate clairement 
qu’uniquement le modèle de Lamb tient compte du caractère ondulatoire des déplacements 
causés par la charge dynamique. Cette qualité possédée uniquement par le modèle de Lamb le 


rend le modèle le plus proche de la réalité. 


Fig. IIL 11 : Allure des déplacements verticaux de la surface du semi-infini élastique due à une charge 
dynamique donnée par : 


(D modèle de Lamb 


(ID) modèle de Boussinesq 
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IIL7. Conclusion : 


Ce chapitre est consacré à l’étude de la fonction de Green définissant les déplacements verticaux 
de la surface du semi-infini élastique de propriétés inertielles (modèle de Lamb) causés par une 
charge statique ou dynamique. L’étude de la fonction de Green est essentielle pour l’étude 
dynamique des différents problèmes de contact. Cette étude de la fonction de Green vise en 
particulier l’étude dynamique d’une plaque rectangulaire reposant sur la surface du semi-infini 
élastique de propriétés inertielles qui sera apportée dans les chapitres qui suivent ce chapitre. A 
noter que le problème traité dans ce chapitre représente le plus grand challenge de cette thèse 
de doctorat puisqu'il est lié à des difficultés mathématiques majeures. La difficulté de l’étude 
de la fonction de Green s’est compliquée encore lors de l’apparition du problème de singularité 
lors de l’évaluation de certaines intégrales associées à cette étude. Cette situation nous a obligé 
en premier lieu à surmonter ce problème de singularité qui s’est avéré n’est possible que par la 
solution analytique ou semi-analytique de certaines intégrales contenant des fonctions 
spéciales. Ensuite, exécuter le reste des calculs mathématiques et surmonter aussi leurs 
difficultés et enfin de compte construire des codes de calcul appropriés pour arriver aux résultats 
recherchés. Tous les efforts fournis ont abouti à la résolution des différents problèmes 
rencontrés et finalement les résultats trouvés sont très satisfaisants et de très bonne précision et 
qui sont tous en totale concordance avec la réalité des choses dans ce domaine. En outre, les 
formules finales obtenues ont été arrangées de façon être plus compatibles pour différentes 
applications d’ingénierie, non seulement dans le domaine des problèmes de contact, mais aussi 
dans d’autres domaines comme le cas de la propagation des ondes dans différents milieux 


continus due à plusieurs facteurs extérieurs. 


Enfin, la comparaison des résultats obtenus effectuée avec un autre modèle a montré que le 
modèle de Lamb étudié est le modèle le plus proche de la réalité en raison de certaines qualités 
qu’il possède. Par exemple, il est le seul modèle qui tient en compte du caractère ondulatoire 
des déplacements ce qui lui classe comme étant le meilleur parmi les autres modèles existants 


dans ce domaine. 
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CHAPITRE IV : Présentation de l’approche de calcul utilisée dans l’étude dynamique des 
plaques reposant la surface du semi-infini élastique de propriétés 
inertielles (modèle de Lamb) 


IV.1. Introduction : 


L’étude des problèmes de contact comme le calcul des plaques en interaction avec des milieux 
élastique en statique ou en dynamique est liée à de difficultés mathématiques majeures. 


Ce chapitre est donc consacré à la présentation d’une approche de calcul proposée pour l’étude 
dynamique d’une plaque rectangulaire reposant sur la surface du semi-infini élastique de 
propriétés inertielles (modèle de Lamb). L'approche proposée se base sur la méthode des forces 
connue dans le domaine des sciences liées à la résistance des matériauxe RDM qui est connue 
aussi dans certaines littératures sous le nom de la méthode de Zhemochkin[72]. L'application 
de cette approche nécessite en premier lieu la discrétisation du système étudié plaque reposant 
sur la surface du semi-infini élastique de propriétés inertielles en un nombre fini d’éléments 
identiques. Donc le contact continu entre la plaque et la surface du semi-infini élastique est 
remplacé par un contact partiel assuré par des liaisons se trouvant au niveau des centres des 
éléments de discrétisation. Cette opération aboutit à un système d'équations canoniques avec 
plusieurs paramètres. Certains de ces paramètres sont connus apriori (généralement les 
caractéristiques mécaniques et géométriques du système étudié) et d’autres sont à déterminer 
par différents calculs. Les deux paramètres du système d’équations canoniques à définir et qui 
représentent le défi de cette étude sont : l’étude de la fonction de Green, représentant les 
déplacements verticaux de la surface du semi-infini élastique, et la détermination des déflexions 
de la plaque. 


Enfin, la résolution du système d’équations canoniques résultant permet la détermination des 
forces réactives dans la zone de contact située entre la plaque flexible et la surface du semi- 
infini élastique de propriétés inertielles. Ceci permettra ensuite de déterminer le reste des autres 
différentes entités physique, à savoir, les déplacements verticaux de la zone de contact, les 
efforts internes dans les sections de la plaque, les valeurs des fréquences propres et les modes 
propres de la plaque étudiée, la réponse dynamique de la plaque etc. 


IV.2. Problématique du chapitre : 


Le problème posé de cette étude est la formulation de l’approche proposée qui sera utilisée dans 
l’étude dynamique d’une plaque rectangulaire reposant sur la surface du semi-infini élastique 
de propriétés inertielles (modèle de Lamb). Cette étude dynamique comporte la détermination 
des valeurs des fréquences propres et les modes propres de la plaque, ainsi que sa réponse 
dynamique due à de différents types de charges extérieures harmoniques appliquées sur la 
plaque. La méthode de résolution choisie, nécessite la discrétisation du système étudié, c’est-à- 
dire le passage du continu au discret est inévitable puisque la résolution continue est 
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pratiquement impossible pour ce genre de problèmes de contact. La structure étudiée est donc 
discrétisée en un nombre d’éléments rectangulaires identiques dont le nombre est proportionnel 
à la précision des résultats souhaitée. Cette opération de résolution, basée sur la discrétisation 
du système étudié, aboutit enfin de compte à un système d'équations canoniques de plusieurs 
paramètres à déterminer apriori avant sa résolution finale. La détermination des différents 
paramètres du système d’équations canoniques nécessite la combinaison de plusieurs études 
séparées liées à chaque paramètre et qui sont enfin de compte couplées pour formuler le système 
final. On a choisi que le système final soit formulé de telle sorte que la première solution donne 
les valeurs des forces réactives dans la zone de contact entre la plaque et la surface du semi- 
infini élastique de propriétés inertielles et ensuite la détermination des fréquences propres et les 
modes propres de la plaque, ainsi que sa réponse dynamique. 


IV.3. Discrétisation du système étudié (plaque reposant sur la surface du semi-infini 
élastique de propriétés inertielles) : 


Soit une plaque mince de masse m et de rigidité cylindrique D repose sur la surface du semi- 
infini élastique de propriétés inertielles comme la montre la figure IV.1 suivante : 


Fig. IV. 1 : Plaque reposant sur la surface du semi-infini élastique de propriétés inertielles 


Divisons la plaque en un nombre n d’éléments rectangulaires de dimensions identiques, et au 
centre de chaque élément on met une liaison rigide à travers laquelle se réalise le contact de la 
plaque avec la surface du semi-infini élastique. C’est-à-dire, le contact continu entre la plaque 
et la surface du semi-infini élastique est remplacé par un contact partiel accompli au niveau des 
haisons se trouvant aux centres des éléments. Admettons que la masse de chaque élément est 
concentrée en son centre (fig. IV.2). 
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Fig. IV. 2 : Plaque discrétisée en un nombre d’éléments rectangulaires identiques reposant sur la surface 
du semi-infini élastique 


Pour des raisons de simplification du calcul, la numérotation des éléments de la plaque 
discrétisée est choisie de la manière suivante : 


Fig. IV. 3: Choix de la numérotation des éléments de la plaque discrétisée 
Où : 

r : nombre d’éléments selon l’axe Ox ; 

s : nombre d’éléments selon l’axe Oy. 
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IV.4. Formulation de l’approche : 


L’approche proposée pour l’étude dynamique de la plaque de forme rectangulaire reposant sur 
la surface du semi-infini élastique de propriétés inertielles se base sur le principe de la méthode 
des forces connue dans la science de la RDM. Cette méthode connue aussi sous le nom de la 


méthode de Zhemochkin, et son principe est illustré sur la figure IV.4. 


Xa()  Xa(t) Xuu(t  X,(t) 


Fig. IV. 4 : Discrétisation du système étudié (plaque en interaction avec la surface du semi-infini 
élastique) 


A noter que notre approche suppose que les forces d'inertie J; (1) vibrant les masses sont 
appliquées uniquement sur la plaque, alors que les efforts de liaisons X, (1) sont appliqués sur 


la plaque et sur la surface du semi-infini élastique [73], [72], [67], [74]. 


Le système d’équations canoniques permettant l’étude dynamique de la plaque de forme 
rectangulaire reposant sur la surface du semi-infini élastique de propriétés inertielles est le 


suivant [75], [76], [72], [2] : 
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,+W,)X,()-D 7,048 d,()+6, (1) +u(r)+A% =0 i=L1...,n 


x,(9-7,(11]6,=L4,(0) (IV.1) 


X,(6)-J,(0)]=mu(r). 


Les définitions détaillées de différents paramètres du système d’équations canonique (IV.1) 

sont les suivantes : 

v. : fonction définissant le déplacement vertical de la surface du semi-infini élastique au point 
i dû à la force X, =1 appliquée au point j de la même surface. Pour sa détermination, on 


doit étudier la fonction de Green ayant comme modèle le cas du semi-infini élastique de 
propriétés inertielles (modèle de Lamb). Cette fonction représente le premier grand défi en 
matière de complexité de cette étude dynamique. Le chapitre II est complètement consacré 
à la détermination de ce paramètre ; 

W: fonction définissant la déflexion de la plaque au point à due à la force X; =1 appliquée 


au point ; de la plaque. Pour sa détermination, on doit appliquer la méthode énergétique 

de Ritz adoptant la solution de Clebsch de l’équation différentielle des déformées de la 
plaque. Cette fonction considérée comme étant le second grand défi de cette étude 
dynamique. Le chapitre IT est complètement consacré à la détermination de ce paramètre ; 

X, (r): effort de liaison appliqué sur la plaque et sur la surface du semi-infini élastique 
(inconnue). Ces efforts de liaison représentent les forces éactives dans la zone de 

contact entre la plaque et la surface du semi-infini élastique de propriétés inertielles ; 

À» : fonction représentant la déflexion de la plaque au point : due à l’action de la charge 
extérieure appliquée au point p coïncidant avec l’un des centres de gravité des éléments 


de discrétisation de la plaque. Donc, ce paramètre n’est tenu en compte que pour le cas 
des vibrations forcées de la plaque; 


1e (1): force d’inertie vibrant la structure étudiée, elle est supposée appliquée uniquement sur 


la plaque ; 


L.. : bras de levier du centre de gravité de l’élément à par rapport à l’axe Ox ; 
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l,, : bras de levier du centre de gravité de l’élément à par rapport à l’axe Oy ; 
Pr (r), D, (1) : angles de rotation de la plaque par rapport aux axes Ox et Oy (inconnues) ; 


U (1) : déplacement vertical initial de la plaque toute entière, considérée comme étant un seul 
bloc rigide (inconnue) ; 


Lors Lo, : Moments d'inertie de la plaque par rapport aux axes Ox et Oy; 


Ox”? 


m : masse totale de la plaque : 


Por (1) : accélération angulaire de la plaque par rapport à l’axe Ox ; 
D, (1) : accélération angulaire de la plaque par rapport à l’axe Oy ; 


U (1) : accélération verticale de la plaque par rapport à l’axe Oy ; 


Les accélérations angulaires et verticales de la plaque peuvent être exprimées par [74] : 


.e 2 LI] d’ .e 2 
L d'hx. > Pos u (1) = d " | 
dt 


IV.S5. Cas des vibrations libres de la plaque : 


La vibration libre de la plaque est supposée sous forme harmonique, alors on pourra écrire [77], 
[78] : 


v,()= De 

W,()=We; 

X ,()=X;e""; 

BD = he". (V2) 
hO = he": 

u,(t) =ue"”; 

JO =Je". 


Du fait que la plaque est réellement en parfait contact avec la surface du semi-infini élastique 
de propriétés inertielles, alors le déplacement de ce dernier est considéré égal à la déflexion de 


la plaque (W, = v,). Donc, la force d'inertie (J ;) peut être donnée par l'expression suivante 


[79] : 
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d’W, (t d°v,(t 
#, — M, mn ) _ | 1 ) _ M ;œ° ve" (IV.3) 


Reportons (1V.2) et (IV.3) dans le système (IV.1) et considérons A, =0, ce dernier devient : 


>= 


nm 1! IX 


(v, +W, X, D JL: +10, +04 +u, =0; 


Xe 10; AV4) 


> X;—J;] Ly =, Lys 


n .e 


X,-J,|=mu,, 


En tenant compte de toutes les définitions précédentes de tous les paramètres du système 
d’équations canonique (IV.3) apportées dans ce chapitre ou dans les chapitres II et IIL et après 
plusieurs transformations et simplifications mathématiques, le système peut être représenté sous 


la forme matricielle suivante : 


Ki; ea 7 bi 1 X, [O0 
Ki Km Cm fm 1 NX, TUTO vs 
Kia Kiun A 0 0 |#, | |0 
K,42,1 F0 4 0 p, 0 
Le Ka, 0 0 mo JL% | [0] 


Les différentes transformations et simplifications mathématiques effectuées et grâce à 
l’utilisation du logiciel Mathematica, les termes de la matrice du système (1V.5) peuvent être 


définis par les expressions suivantes : 
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K,, =W,; var, -MaxY(F,) 1) ie[Ln|]; jef[in|:; 
p= 


cal Man x DL, XE,.) je[La]: 


r 
p=l 


Kyo, 


=(,-Mao x Ÿ (1, XF) jelLnl: 


p=1l 


fe _1-Maw xŸ(F,.,) jelbaf 


p=l 


À, a 
3 
Il 2 2 
À,=-a mo"; 
3 


m : la masse total de la plaque : 
M: la masse élémentaire correspondant à l’élément ; ; 


—k il 


KES 
2rG. 
Q, : la surface de l’élément à ; 


k : paramètre définissant le caractère ondulatoire des déplacements est donné dans le chapitre 


IT ; 


E, 


G, : module d’élasticité transversale du semi-infini élastique et défini par : G, = Zn) 
+" 


L.) 


E, et v, : module d'élasticité et coefficient de Poisson du semi-infini élastique respectivement ; 


IV.6. Fréquences propres de la plaque étudiée : 


Le système d’équation sous la forme matricielle (1V.5) issu du système d’équations canonique 
(IV.1) après avoir subi plusieurs transformations mathématiques permet de déterminer les 
fréquences propres de la plaque étudiée. Le processus de détermination des fréquences propres 
de la plaque nécessite la détermination des racines de l’équation du déterminant de la matrice 


[A] en égalant son expression à zéro, c’est-à-dire : 


det[ A]=0 (IV.6) 
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La détermination des racines de l’équation du déterminant est l’opération de calcul de cette 
étude dynamique de la plaque la plus compliquée et la plus lente vu la complexité de son 
expression ce qui nécessite beaucoup de temps de calcul et aussi une machine de calcul très 
puissante surtout lorsque la taille de discrétisation est grande. Plus de détail sur cette opération 
de calcul sera apporté dans le chapitre V. 


IV.7. Modes propres de la plaque étudiée : 


Après avoir déterminé les valeurs des fréquences propres de la plaque étudiée, on passe 
automatiquement à la détermination de ses modes propres. Chaque mode propre correspond à 
une valeur de fréquence propre. Plus de détail sur cette opération de détermination des modes 
propres de la plaque sera apporté dans le chapitre V. Enfin, le but de détermination des 
fréquences propres et des modes propres de la plaque étudiée est de savoir ses propriétés 
approfondies et par conséquent éviter l’apparition des phénomènes de résonnance lors de sa 
mise en application. 


IV.8. Réponse dynamique de la plaque étudiée : 


L’étape finale de cette étude dynamique de la plaque rectangulaire reposant sur la surface du 
semi-infini de propriétés inertielles est la détermination de sa réponse dynamique due à l’action 
de plusieurs types de charges extérieures harmoniques. Autrement-dit, savoir son 
comportement et les déplacements verticaux de la plaque étudiée au moment de l’action de la 
charge extérieure. Plus de détail sur la détermination de la réponse dynamique de la plaque due 
aux différents types de charges harmoniques extérieures sera apporté dans le chapitre V. 


IV.9. Conclusion : 


L'approche proposée, dans le cadre de la préparation de ce travail de thèse de doctorat, pour 
l’étude dynamique d’une plaque de forme rectangulaire reposant sur la surface du semi-infini 
élastique de propriétés inertielles (modèle de Lamb) a finalement permis la détermination de la 
solution recherchée. En se basant sur la discrétisation du système étudié et en application du 
principe de la méthode des forces connue dans la littérature et après introduction de plusieurs 
études intermédiaires définissant certains paramètres, l’étude finale s’est formulée sous forme 
d’un système d’équations sous forme matricielle. Ce système matriciel final a permis l’étude 
dynamique comportant le calcul des forces réactives dans la zone de contact et puis la 
détermination des valeurs des fréquences propres de la plaque, la détermination de ses modes 
propres et de sa réponse dynamique due aux différents types de charges harmoniques 
extérieures. 
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CHAPITRE V : Calcul dynamique d’une plaque rectangulaire reposant sur la surface du 
semi-infini élastique de propriétés inertielles (modèle de Lamb). 


V.1 Introduction : 


Les plaques de n’importe quelle forme et constituées de n’importe quel matériau ce sont des 
éléments structurels très utilisés dans de nombreux domaines du génie mécanique, du génie 
civil, de l’aéronautique, de l’industrie automobile etc. Donc, leur utilisation ne cesse de grandir 
avec le développement technologique, chose qui impose la bonne maitrise de leur 
comportement en statique ou en dynamique. Les plaques de formes rectangulaire en interaction 
avec des milieux élastiques fait partie de cette problématique et donc l’approfondissement de 
leurs études dynamiques est nécessaire pour, non seulement, maitriser parfaitement leur 
comportement dynamique, mais aussi pour aboutir à une conception pour une application plus 
adéquate lors de leur mise en service. Ce type de plaques intervient dans le domaine des 
problèmes de contact qui n’est pas étudié complètement à nos jours vu sa complexité et son 
accompagnement de difficultés mathématiques majeures. 


Dans ce contexte, ce chapitre sera focalisé sur le calcul dynamique d’une plaque de forme 
rectangulaire reposant sur la surface du semi-infini élastique de propriétés inertielles (modèle 
de Lamb). En vue de contribuer à résoudre ce problème de contact non traité suffisamment, une 
approche semi-analytique, bien décrite dans les chapitres précédents, sera utilisée pour son 
analyse dynamique. Cette étude dynamique se base sur la discrétisation de la structure étudiée, 
plaque rectangulaire reposant sur la surface du semi-infini élastique de propriétés inertielles, en 
un nombre d’éléments identiques fini. Les deux études auxiliaires essentielles faisant partie de 
cette étude dynamique sont : la méthode énergétique de Ritz et l’étude de la fonction de Green. 
Le couplage de ces deux études auxiliaires essentielles avec d’autres études secondaires est 
accompli par la méthode des forces connue dans la littérature et connue aussi sous le nom la 
méthode de Zhemochkin. 


L’application de l’approche sous sa forme finale permet la détermination des forces réactives 
dans la zone de contact entre la plaque et la surface du semi-infini élastique ce qui permettra 
ensuite de déterminer les valeurs des fréquences propres de la plaque, ses modes propres 
appropriés à chaque fréquence propre, ainsi que sa réponse dynamique due à l’action de 
différents types de charges extérieures harmoniques. 


V.2. Problématique du chapitre : 


Le problème posé de cette étude dynamique est de discrétiser la structure étudiée, une plaque 
mince isotrope de forme rectangulaire reposant sur la surface du semi-infini élastique de 
propriétés inertielles (modèle de Lamb), pour déterminer les valeurs de ses fréquences propres, 
ses modes propres associés à chaque fréquence propre, ainsi que sa réponse dynamique due à 
l’application de différents types d’excitations extérieures harmoniques. L’amortissement et les 
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forces de frottement dans la zone de contact entre la plaque et la surface du semi-infini élastique 
sont négligés. 


V.3. Discrétisation du système étudié : 


L’application de l’approche utilisée nécessite le passage du système continu au système discret, 
c’est-à-dire le système étudié, plaque mince de forme rectangulaire reposant sur la surface du 
semi-infini élastique de propriétés inertielles (modèle de Lamb) doit être discrétisé en un 
nombre n d’éléments rectangulaires identiques. 


Donc le remplacement du calcul continu par un calcul fini implique la détermination des forces 
de réaction dans la zone de contact uniquement aux points situant au niveau des centres des 
éléments auxquels la plaque est divisée. Le calcul des autres entités physique, à savoir les 
valeurs des fréquences propres et les modes propres de la plaque, ainsi que sa réponse 
dynamique, se fait par l’application des différentes formules de la théorie d’élasticité. 


Pour simplifier un peu la tâche du calcul très compliqué de ce problème de contact, on considère 
une plaque mince homogène de forme carrée. La numérotation choisie des éléments de la 
plaque carrée discrétisée est illustrée sur la figure (V.1). On constate bien que tous les éléments 
de discrétisation de la plaque sont identiques et leur nombre total considéré est égal à n. Ceci 
implique automatiquement que le nombre d’éléments dans les deux directions est le même, 


c’est-à-dire r=s=4V/n. 


Fig. V. 1 : Numérotation des éléments de la plaque discrétisée 
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V4. Formulation du système d’équations canoniques : 


Le système d’équations canoniques permettant l’étude des vibrations libres de la plaque mince 
carrée reposant sur la surface du semi-infini élastique de propriétés inertielles (modèle de 
Lamb) [72] est bien décrit dans le chapitre (IV). Sa forme matricielle finale vaut : 


Ki Kn “ A 1 X, [0] 
K,1 nn É. ke Il X, : (V 1) 
n+1,1 Ki n À, 0 0 ps E ; | 
Ku, K,42 n 0 À, 0 p, 
lé Ky43 n 0 0 mo JL Un | 0] 


Tous les paramètres du système d’équations sous forme matricielle (V.1) sont définis dans le 
chapitre (IV). 


V.5. Détermination des valeurs des fréquences propres de la plaque reposant sur la 
surface du semi-infini élastique de propriétés inertielles : 


La détermination des valeurs des fréquences propres de la plaque mince de forme carrée 
reposant sur la surface du milieu élastiques de propriétés inertielles est le point le plus 
compliqué à réaliser de cette étude toute entière. Vu le long temps de calcul exigé, la nécessité 
de l’utilisation d’une machine de calcul très puissante est donc inévitable pour accomplir cette 
tâche de calcul. A noter aussi que plus le maillage de la plaque est fin, plus de le temps de calcul 
devient important ce qui en parfaite concordance avec le principe général des méthodes 
numériques. Autrement dit, la précision des résultats de calcul est proportionnelle avec la taille 
de discrétisation de la structure étudiée. C’est-à-dire, plus la taille de discrétisation de la plaque 
n est grande, plus la précision elle est bonne, mais au prix d’un effort et de temps de calcul 
plus importants. 


Un autre challenge s’est fait face lors de la détermination des valeurs des fréquences propres de 
la plaque, c’est la convergence de la valeur de chaque fréquence propre vers une valeur stable. 
Pour surmonter ce chalenge, on est contraint de finir le maillage de la plaque jusqu’à ce que la 
stabilité de toutes les valeurs des premières fréquences propres soit observée. Pour atteindre la 
stabilité recherchée, le calcul donc est accompli pour différentes tailles de maillage allant de 
100 éléments jusqu’à 900 éléments. Cette opération est exécutée grâce au développement des 
codes de calcul spécifiques à l’aide du logiciel de calcul numérique et formel « Mathematica » 
et grâce aussi au Plateau Technique des Calculs Intensifs de l’université 8 Mai 1945 Guelma. 
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La détermination des valeurs des fréquences propres de la plaque étudiée se fait par la manière 


suivante : 


Le système étudié, plaque mince carrée reposant sur la surface du semi-infini élastique de 
propriétés inertielles, est discrétisé en un nombre n d’éléments carrés identiques, figure (V.1). 
Tenons compte des valeurs suivantes des propriétés géométriques et mécaniques de la plaque 


et du semi-infini élastique 


a=b=1m : longueur et larguer de la plaque ; 

v =1/3 : coefficient de Poisson de la plaque 

E=2x10"N/m° : module d'élasticité de la plaque ; 

h=0.01m : épaisseur de la plaque ; 

E,=25x10" N/m° : module d'élasticité du semi-infini élastique ; 
v, =1/3 : coefficient de Poisson du semi-infini élastique ; 

p =1500 Kg / m° : densité du semi-infini élastique. 


Les valeurs des fréquences propres « de la plaque carrée reposant sur la surface du semi-infini 


élastique de propriétés inertielles représentent les racines de l’équation du déterminant de la 


matrice [A] du système d’équations (V.1) : 
det[ A]=0 (V2). 


Vu la complexité et la longueur excessive de l’expression du déterminant de la matrice [A] 


dépassant plusieurs dizaines de pages, on a fait recours à l’application de la méthode dichotomie 
(bissection) [66] et à l’aide de l’utilisation du logiciel du calcul numérique et formel 


« Mathematica » les valeurs recherchées sont finalement trouvées. 


Les convergences des valeurs concernant les six premières fréquences propres vers des valeurs 


stables en Hertz (Hz) sont illustrées dans le tableau (V.1) suivant : 
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2021 


ie ch 0 © 2 © 3 © 4 © 5 © 6 
100 52,30528 | 64,44392 | 70,24914 | 74,16060 | 126,44535 | 250,12467 
121 | 52,23230 | 64,38295 | 69,99672 | 73,979044 | 126,18498 | 249,61034 
144 | 5217675 | 64,34161 | 69,80591 | 73,80554 | 125,71228 | 249,32589 
169 | 52,13349 | 64,30583 | 69,65812 | 73,69786 | 125,54611 | 249,02114 
196 | 52,09916 | 64,27999 | 69,54128 | 73,59239 | 125,26489 | 248,82994 
225 | 5207145 | 64,25722 | 69,44730 | 73,52313 | 125,15276 | 248,63502 
256 | 52,04877 | 64,24002 | 69,37056 | 73,45441 | 124,97232 | 248,50234 
289 | 52,02996 | 64,22465 | 69,30709 | 73,40721 | 124,89319 | 248,37029 
324 | 5201420 | 64,21262 | 69,25398 | 73,35996 | 124,77071 | 248,27514 
361 | 52,00086 | 64,20176 | 69,20910 | 73,32633 | 124,71281 | 248,18161 
400 | 5198946 | 64,19303 | 69,17082 | 73,29247 | 124,62597 | 248,11131 
441 | 5197966 | 64,18507 | 69,13792 | 73,26766 | 124,58234 | 248,04268 
484 | 5197116 | 64,17854 | 69,10942 | 73,24257 | 124,51859 | 247,98939 
529 | 5196374 | 64,17253 | 69,08457 | 73,22374 | 124,48489 | 247,93755 
576 | 5195723 | 64,16752 | 69,06278 | 73,20463 | 124,43675 | 247,89626 
625 | 5195068 | 64,16287 | 69,04356 | 73,19001 | 124,41018 | 247,85615 
676 | 5194638 | 64,15894 | 69,02653 | 73,17512 | 124,37294 | 247,82353 
729 | 51,93584 | 64,15427 | 69,01136 | 73,16353 | 124,35163 | 247,79186 
784 | 51,93088 | 64,14673 | 68.974790 | 73,13294 | 124,28095 | 247,71898 
841 | 51,93086 | 64,14668 | 68,97465 | 73,13287 | 124,28135 | 247,71889 
900 | 51,93084 | 64,14663 | 68,97462 | 73,13283 | 124,28131 | 247,71888 


Tab. V. 1 : Convergence des valeurs des six premières fréquences propres de la plaque étudiée vers des 
valeurs stables 


Il est important de mentionner que pour arriver à la convergence des valeurs des fréquences 
propres de la plaque étudiée, il nous a fallu environ sept mois de calcul continu sur le Plateau 
Technique des Calculs Intensifs de l’Université. 


V.6. Détermination des modes propres de la plaque étudiée : 


La détermination des modes propres de la plaque étudiée se fait de la manière suivante : 
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Remplaçons dans le système d’équations de forme matricielle (V.1) la fréquence « par l’une 


de ses valeurs propres. Ensuite, le système matriciel obtenu sera réduit à la taille 
(n+2)x(n+2) au lieu de sa taille initiale (r+3)x(n+3) par la suppression de la première 


ligne du système matriciel et le transfert de la première colonne de la matrice au vecteur libre 


du même système matriciel en normalisant l'inconnue X, = 1. Le système résultant vaut : 
1 


K2 KE, 2y 1 X, Ki 
K,,2 one K, n ee es 1 X, Kii 
1 0 = (V.3) 
K,412 Koun x 0 bn Kyu 
K,4202 ie Kyo n 0 À, 0 b,, K,,21 
2 
Êe Kii3n 0 0 mao | L Un | Le | 


La résolution du système d’équations sous forme matricielle résultant (V.3) permet de 
déterminer les valeurs des autres inconnues X;. Cependant, la détermination des déplacements 
verticaux de la plaque étudiée définissant la forme propre de la plaque étudiée appropriée à 


chaque valeur de sa fréquence propre se fait par l’application de la formule (V.4) suivante [67], 
FA 


= — x—DX, F;; (V4) 


v, : déflexion de la plaque étudiée au point i formant avec ses autres déflexions sa forme 
propre; 


k : paramètre définissant le caractère ondulatoire des déplacements verticaux de la surface du 


2 
D 1 - OX 
semi-infini élastique, donné par : k = FL : 
0 


& : fréquence d'excitation ; 
p : densité du semi-infini élastique ; 


E, 
2(1+v,) 


? 


G, : module d’élasticité transversale, défini par : G, = 


E, et v, : module d'élasticité et coefficient de Poisson du semi-infini élastique respectivement ; 


Q, : la surface de l’élément 1 ; 
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F, : fonction de Green définissant les déplacements verticaux de la surface du semi-infini 
élastique de propriétés inertielles (modèle de Lamb) au point i dû à la force X ; appliquée 


au point j de la même surface. Le détail de sa détermination est donné dans le chapitre 
IL. 


Cette opération de détermination des modes propres de la plaque étudiée est répétée pour 
chaque valeur de sa fréquence propre, ce qui a permis la détermination du reste de ses autres 
modes propres. 


La figure V.2 tabule les formes des modes propres correspondant aux seize premières valeurs 
des fréquences propres de la plaque étudiée. À noter que le maillage de la plaque étudiée pris 
en compte dans ce calcul est très fin où n =900 éléments. 


œ 1=51,93083502352 o 2=64,1466333568 


7 “& 


o 3=68,97461734712 o 4=73,13282933831 
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o 5=124,28131334483 o 6=—247,71888233721 


@ 7 = 344,12066601216 œ 8=360,851501002907 


© 9 = 445,963996723294 © 10=451,56233813241 
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o 11 =468,3279841952 @ 12 = 485,244140684604 


© 13 = 647,58467203006 © 14—773,7379338033 


© 15 =829,0718657337 © 16 =965,3014776669 


Fig. V.2 : Les formes des modes propres de la plaque étudiée correspondant à ses seize premières 
fréquences propres 
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On constate bien, d’après la figure V.2, que la forme des trois premiers modes propres de la 
plaque étudiée représentent les phénomènes du pompage (translation de la plaque 
verticalement) et le tangage (rotation de la plaque par rapport aux axes des coordonnées Ox et 


Oy). Il est aussi bien clair que les autres formes des modes propres de la plaque changent de 
signe par rapport au plan de coordonnées Oxy et que ce changement de signe, dans son 


ensemble, est proportionnel avec l’augmentation de la valeur de la fréquence propre 


correspondante. Ceci est en parfaite concordance avec la théorie de la dynamique des structures. 


V.7. Vérification des valeurs des fréquences propres par le principe de la méthode de 


superposition modale : 


Afin de s’assurer de la certitude des résultats obtenus concernant les valeurs des fréquences 
propres de la plaque étudiée, on procède à leur vérification en appliquant le principe de la 


méthode de superposition modale. Pour ce faire, on doit déterminer le signal T ;(r) représentant 


la somme des amplitudes des seize premiers modes propres en un point donné j de la plaque. 


La formule du signal pour notre cas est donnée par [80]: 


n=16 
TO = Ÿ v,,xcos(27wf) (V5) 


i=1 
v;,: amplitude du mode propre i de la plaque au point j ; 
©, : fréquence propre de la plaque ; 

t : temps de réponse modale ; 


n=16 : nombre de modes propres pris ; 


L'application de la formule (V.5) sur un intervalle de temps donné permet de tracer le graphe 


du signal recherché au niveau du centre de l’élément de la plaque portant le numéro j = 807 


choisi arbitrairement. 
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v (mm) 


0.002 


0.001 


LL. 


-—0.001 | 


-—0.002 


Fig. V.3 : Signal des amplitudes modales en un point de la plaque donné par la formule (V.5) 


L’étude du problème inverse de ce signal, par l’application de la transformée de Fourrier via un 
code de calcul approprié élaboré sur Mathematica, nous a permis d’aboutir au spectre localisant 


les valeurs des fréquences propres de la plaque en fonction de la magnitude de la FFT. 


La figure (V.4) représente la transformée de Fourrier du signal des amplitudes modales. La 
lecture du spectre obtenu permet la localisation des seize premières fréquences propres de la 


plaque étudiée. 


Magnitudes de la FFT 


À 
0.3 
02F 


Les cinq premières fréquences propres 


01+ 


L L 1 L 1 L L JL 
200 400 600 800 1000 
Fréquences (Hz) 


Fig. V. 4 : Spectre localisant les seize premières fréquences propres de la plaque étudiée 
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Les magnitudes de la FFT des premières fréquences propres de la plaque ayant des amplitudes 
très grandes par rapport aux autres. De ce fait, les piques des premières fréquences propres de 
la plaque sont bien clairs, comme le montre la figure (V.4), et donc leur localisation sera facile 
à faire. Le reste des piques des autres fréquences propres ne figurent pas sur le spectre de la 
figure (V.4), c’est pourquoi on procède au zoom de tout le spectre sur différentes bandes de 
fréquences afin de les montrer. 


Magnitudes de la FFT Magnitudes de la FFT 


250 300 350 400 450 
Fréquences (Hz) 


40 60 80 100 120 140 
Fréquences (Hz) 


A B 


Magnitudes de la FFT Magnitudes de la FFT 
o, =485 Hz ©; =647H LUE =829 


0.00020 Ÿ | 0.00020 


0.00015 0.00015 


0.00010 0.00010 


0.00005 0.00005 


450 500 550 600 650 700 700 750 800 850 900 950 1000 


Fréquences (Hz) Fréquences (Hz) 


Fig. V. 5 : Zoom spectral pour la localisation des fréquences propres de la plaque sur différentes bandes 
de fréquences 
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Le zoom du spectre sur différentes bandes de fréquences nous a permis l’illustration de tous les 
piques des fréquences propres de la plaque et par conséquent la localisation de leurs valeurs. 
Le pique correspondant à la fréquence propre @,, = 451 Hz est mal apparu sur la figure (V.5.C), 


du fait que la différence des amplitudes des magnitudes de la FFT des fréquences @,et @,,est 


très grande ce qui a rendu le pique de la fréquence @,, noyé dans le pique de la fréquence @,. 


Magnitudes de la FFT 


440 445 450 


455 Fréquences (Hz) 


Fig. V. 6: Zoom spectral pour la localisation de la fréquence propre @,, = 451 Hz 


Enfin, cette vérification des valeurs des fréquences propres de la plaque étudiée par le principe 
de la méthode de superposition modale nous a confirmé la fiabilité de l’approche utilisée et par 
conséquent la certitude des résultats obtenus. 


V.8. Vibrations forcées de la plaque rectangulaire reposant sur la surface du semi-infini 
élastique de propriétés inertielles dues aux différents types de charges extérieures 
harmoniques : 


L’approche utilisée permet, non seulement, l’étude des vibrations libres de la plaque qui est la 
détermination de ses fréquences propres et ses modes propres, mais aussi l’étude de ses 
vibrations forcées. Donc, la détermination de la réponse dynamique de la plaque rectangulaire 
reposant sur la surface du semi-infin1 élastique de propriétés inertielles due aux différents types 
des excitations extérieures harmoniques est accomplie à l’aide de cette approche. La réponse 
dynamique de la plaque consiste à déterminer ses déplacements verticaux causés par l’action 
d’une ou de plusieurs excitations extérieures harmoniques durant le temps d’excitation. 
Autrement dit, les déplacements de la plaque varient avec le temps d’excitation et plusieurs 
facteurs comme : le nombre d’excitation agissant sur la plaque, l’amplitude et la fréquence 
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d’excitation etc., influent considérablement sur leur variabilité. Plusieurs cas de la réponse de 
la plaque ont été étudiés et les détails de quelques cas sont présentés dans les sections suivantes. 


V.8.1. Réponse de la plaque étudiée due à une excitation extérieure harmonique : 


Considérons les mêmes valeurs des caractéristiques géométriques et mécaniques de la plaque 
et du semi-infini élastique données dans ce chapitre. La plaque étudiée considérée carrée et elle 
est discrétisée en n =900 éléments carrés identiques. 


Fig. V.7 : Numérotation des éléments de discrétisation d’une plaque carrée reposant sur la surface du 
semi-infini élastique de propriétés inertielles sollicitée par une seule excitation extérieure 


harmonique verticale 


Cas 8.1.1 : Plaque sollicitée par une seule excitation extérieure harmonique 


La plaque est sollicitée par une seule excitation extérieure harmonique verticale appliquée au 
centre de l’élément portant le numéro 205, figure (V.7). Supposons que l’excitation extérieure 
harmonique se varie selon la loi suivante : 
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P, = Pos = Fcos(2xxx xt) (V.6) 
Où : 

P, : amplitude de l’excitation ; 

© : fréquence de l’excitation ; 

27@ : pulsation de l’excitation ; 

p : numéro de l’élément où l’excitation extérieure est appliquée, dans ce cas p = 205 ; 


Î : temps de l’excitation. 


L’étude des vibrations forcées de la plaque nécessite la tenue en compte d’un autre paramètre 
À,, dans le système d’équations sous forme matricielle. Ce dernier prend la forme suivante : 


Ki nr de & À | #Æ [A,, | 
Ke “a D da 4 1 X, L A,, . 
Kunst Kuun À 0 0 |, 0 | 
Kiyor Kyo 0 À, 0 p, 0 
Le EU 0 0 mo || ü | L 0 | 


Dans ce cas, le paramètre A,, est déterminé par la formule suivante (V.8) [67], [72] : 


À, => WP (V8) 


p=l 


Dans le cas d’une seule excitation extérieure harmonique appliquée au centre de l’élément 
n=900 


portant le numéro 205, (V.8) devient : À; = > WP = W 05 Pos 
p=l 


» : déflexion de plaque au point { causée par l’excitation extérieure P, appliquée au point 
p de la plaque (dans ce cas, p = 205, donc l’excitation extérieure P, est appliquée au 
centre de l’élément portant le numéro 205). Le détail de détermination des déflexions de 


la plaque est donné dans le chapitre I. 


La résolution du système d’équations sous forme matricielle (V.7) permet la détermination des 
valeurs des inconnues X,, i=1,:::,900. Ces inconnues X; représentent les forces de réaction 
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dans la zone de contact qui de leur tour varient avec le temps d’excitation. Finalement, la 
réponse de la plaque à chaque instant d'excitation est déterminée par (V.9) : 


_k 1 n=900 
y, x— SX (:)F,;i=1..,n (V.9) 
27G, Q, = | 


Considérons les valeurs suivantes de l’amplitude d’excitation : Æ, =10000 N et de la fréquence 
d’excitation : © = 400 Hz. 


La figure (V.8) représente la réponse de la plaque en (2D) au niveau du point d’application de 
l’excitation extérieure en fonction du temps d’excitation sur un intervalle de temps donné. Le 
temps d’excitation est exprimé par : 


t,=t,+ jxAt 
t, =0S : temps initial ; 


At =0,000555S$ : pas du temps considéré. 


v (mm) 


t(s) 


0.0015 


Fig. V. 8 : Réponse de la plaque au niveau du point d’application de l’excitation extérieure durant un 
intervalle de temps de son excitation 


Cas 8.1.2 : Comparaison de la réponse de deux différents points de la plaque : 


La figure (V.9) représente la comparaison de la réponse de la plaque au niveau du point 
d’application de l’excitation extérieure avec un autre point de la plaque portant le numéro 697 
durant un intervalle de temps de son excitation. 
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v (mm) 


t(s) 


——— la réponse du point d’application de l’excitation extérieure 


______ la réponse de l’élément portant le numéro 697 


Fig. V. 9 : Comparaison de la réponse du point où l’excitation extérieure est appliquée avec un autre 
point de la plaque durant un intervalle de temps de son excitation 


Cas 8.1.2 : Variation de la réponse de la plaque en fonction de l’amplitude de l’excitation 
extérieure : 


La figure (V.10) représente la variation de la réponse de la plaque au niveau du point 
d’application de l’excitation extérieure en fonction de la variation de son amplitude durant un 
intervalle de temps de sa sollicitation. 


v (mm) 


t(s) 


Fig V. 10 : Variation de la réponse de la plaque au niveau du point d’application de l’excitation 
extérieure en fonction du changement de l’amplitude de l’excitation 


89 


CHAPITRE V : Calcul dynamique d’une plaque rectangulaire reposant sur la 2021 
surface du semi-infini élastique de propriétés inertielles (modèle de Lamb). 


Cas 8.1.4 : Variation de la réponse de la plaque en fonction de l’amplitude et de la 


fréquence de l’excitation extérieure : 


La figure (V.11) représente la variation de la réponse de la plaque au niveau du point 
d’application de l’excitation extérieure en fonction du changement de deux paramètres de 
l’excitation extérieure qui sont l’amplitude de l’excitation et la fréquence de l’excitation durant 


un intervalle de temps de son excitation. 


v (mm) 


t(s) 


—— P=1000N @ =150 
—— P=15000 N; &, =250 H 
— B=20000 N; « =400 Hz 


Fig. V. 11 : Variation de la réponse de la plaque au niveau du point d’application de l’excitation 
extérieure en fonction de deux paramètres : l’amplitude de l’excitation et la fréquence de 
l’excitation 


Cas 8.1.5 : Visualisation (3D) de la réponse de la plaque toute entière : 


La figure (V.12) représente la réponse de la plaque toute entière (3D) sur un intervalle de temps 


de son excitation en considérant le même pas du temps. 
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CHAPITRE V : Calcul dynamique d’une plaque rectangulaire reposant sur la 2021 
surface du semi-infini élastique de propriétés inertielles (modèle de Lamb). 


t=t,+6At t=t,+7A t =t, +8At 


Fig. V. 12 : Visualisation (3D) de la réponse de la plaque toute entière due à une seule excitation 
harmonique appliquée au centre de l’élément portant numéro 205 


V.8.2. Réponse de la plaque étudiée due à deux excitations extérieures harmoniques de 


mêmes et/ou de différentes amplitudes et ayant les mêmes fréquences d'excitation : 


Dans ce cas la plaque est sollicitée par deux excitations extérieures harmoniques de mêmes 
et/ou de différentes amplitudes et de mêmes fréquences d'excitation et qui sont appliquées en 


deux points différents de la plaque, figure (V.13). 
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CHAPITRE V : Calcul dynamique d’une plaque rectangulaire reposant sur la 2021 
surface du semi-infini élastique de propriétés inertielles (modèle de Lamb). 


pe... 


‘ 


‘ F,cos(27 


Fig. V. 13 : Plaque carrée discrétisée en 900 éléments identiques reposant sur la surface du semi-infini 
élastique de propriétés inertielles sollicitée par deux excitations extérieures harmoniques 
verticales 


Cas 8.2.1 : Plaque sollicitée par deux excitations extérieures harmoniques : 


Considérons la même structure étudiée dans la section précédente, sauf que la plaque est 
sollicitée par deux excitations extérieures harmoniques de mêmes amplitudes et de mêmes 
fréquences d'excitation et qui sont appliquées en deux différents points de la plaque portant les 
numéros 218 et 683. Les deux excitations agissant sur la plaque sont exprimées par : 


Pis = Fcos(2x7xwx1) (V.10) 
Ps = P,cos(2xrx xt) | 


Ici : 


F, : amplitude de l’excitation ; 
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CHAPITRE V : Calcul dynamique d’une plaque rectangulaire reposant sur la 2021 
surface du semi-infini élastique de propriétés inertielles (modèle de Lamb). 


© : fréquence de l’excitation ; 


t : temps de l’excitation. 


On considère les mêmes valeurs de l’amplitude Æ, 10000 N et de la fréquence d’excitation 


© = 400 HZ pour les deux excitations. 

Le paramètre A,, dans ce cas de chargement est déterminé par la formule suivante [67], [721 : 
n=900 

À, — D WP, D Wos Pis +W 683 Fes (VID) 
p=l 

La réponse de la plaque à chaque instant de son excitation est déterminée par la même formule 

donnée par l’équation (V.9). 


La figure (V.14) représente la réponse de la plaque au niveau de tous les points situant tout au 
long de sa diagonale où se trouvent les deux points d’application des excitations extérieures 
ayant les mêmes amplitudes et les mêmes fréquences en fonction du temps d’excitation en (2D) 
pour différents instants de temps exprimés par : f, =1, + jxAf. Considérons le même pas du 


temps Af =0,000555$ et le même temps initial #, =0. 


Row(2ra) er) 


v (mm) 


0.6 
i 04 


2 


-04 te ! 


Fig. V. 14 : Réponse de la plaque tout au long de sa diagonale où se trouvent les deux points 
d’application des excitations extérieures en fonction du temps de l’excitation 
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CHAPITRE V : Calcul dynamique d’une plaque rectangulaire reposant sur la 2021 
surface du semi-infini élastique de propriétés inertielles (modèle de Lamb). 


Cas 8.2.2 : Visualisation de la réponse de la plaque des points situant sur sa diagonale : 


La figure (V.15) représente la réponse de la plaque au niveau de tous les points situant tout au 
long de sa diagonale où se trouvent les deux points d’application des excitations extérieures de 
différentes amplitudes et ayant les mêmes fréquences en fonction du temps d’excitation en (2D) 
pour différents instants de temps. F =10000 N ; P. =15000 NW. 


By = PR cos(2xX 7x œxf) (V.12) 
Fix = PL, COS(2X 7 x xt) | 
P, cos(27 wt) 
F cos(27 wt) — 
0.6 
nS 
t is 
oh 04 13 ÈS 
© de : 
t LE 
de 0. 
t 4 
=10 = 0.5 0.5 1.0 


= 06 


Fig. V. 15 : Réponse de la plaque tout au long de sa diagonale où se trouvent les deux points 
d’application des excitations extérieures de mêmes fréquences et de différentes 
amplitudes en fonction du temps de l’excitation 


Cas 8.2.3 : Visualisation (3D) de la réponse de la plaque toute entière : 


La figure (V.16) représente la réponse de la plaque toute entière en (3D) due à ce type de 
chargement extérieur composé de deux excitations harmoniques de mêmes amplitudes et de 
mêmes fréquences pour différents instants de temps de son excitation. 
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CHAPITRE V : Calcul dynamique d’une plaque rectangulaire reposant sur la 2021 
surface du semi-infini élastique de propriétés inertielles (modèle de Lamb). 


t =t, +6At t=t,+7At t =t,+8At 


Fig. V. 16 : Visualisation (3D) de la réponse de la plaque toute entière due à deux charges extérieures 
harmoniques appliquées aux centres de deux différents éléments 


V.8.3. Réponse de la plaque étudiée due à deux excitations extérieures harmoniques de 
mêmes amplitudes et de différentes fréquences d'excitation : 


Prenons maintenant le cas de la plaque sollicitée par deux excitations extérieures harmoniques 
de mêmes amplitudes et de différentes fréquences d'excitation et qui sont appliquées en deux 
différents points de la plaque portant les numéros 218 et 683. Les deux excitations sont 
exprimées par : 


Pie = F,coS(2X 7 x ©, Xf) (V3) 
Ps = PF, COS(2X 7 x ©, Xf) | 
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CHAPITRE V : Calcul dynamique d’une plaque rectangulaire reposant sur la 2021 
surface du semi-infini élastique de propriétés inertielles (modèle de Lamb). 


F, : amplitude de l’excitation ; 


@, : fréquence d’excitation de la charge extérieure appliquée au centre de l’élément de la plaque 


portant le numéro 218 : 


©, : fréquence d’excitation de la charge extérieure appliquée au centre de l’élément de la plaque 


portant le numéro 683 : 


œ : fréquence équivalente des deux fréquences de l’excitation des deux charges extérieures. 
 _ D 2 2, 
Elle est exprimée par : © = /@° +0 ; 


t : temps de l’excitation. 


Cas 8.3.1 : Plaque sollicitée par deux excitations de mêmes amplitudes et de différentes 


fréquences d'excitation : 


Considérons les valeurs suivantes de l’amplitude PF, =10000 N et des fréquences de 


l'excitation &, =150 Hz ; ©, = 400 Hz : d’où ©= 150? +400? = 427,2 Hz ; 


Le paramètre A,, est déterminé par la même formule donnée par l’équation (V.11). Tenons 


compte que les excitations extérieures harmoniques sont appliquées aux centres des éléments 


de la plaque portant les numéros 218 et 683. 


De la même manière on procède à la résolution du système d’équation de forme matricielle 
(V.7) donnant les valeurs des efforts de liaison X, (1) dans la zone de contact en fonction du 
temps de l’excitation et ensuite la réponse de la plaque à chaque instant de l'excitation sera 


déterminée par la même formule (V.9). 


La figure (V.17) représente la réponse de la plaque en (2D) tout au long de sa diagonale où se 
trouvent les deux points d’application des excitations extérieures de même amplitudes et de 


différentes fréquences durant un intervalle de temps de son excitation. 
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CHAPITRE V : Calcul dynamique d’une plaque rectangulaire reposant sur la 2021 
surface du semi-infini élastique de propriétés inertielles (modèle de Lamb). 


P cos(27 ©) f) 


v (mm) 


Fig. V. 17 : Réponse de la plaque tout au long de sa diagonale où se trouvent les deux points 
d’application des excitations extérieures de mêmes amplitudes et de différentes 
fréquences en fonction du temps de l’excitation 


Cas 8.3.2 : Visualisation de la réponse de la plaque aux points d’application des excitations 


extérieures : 


La figure (V.18) représente la comparaison des réponses de la plaque en (2D) au niveau des 
deux points d’application des excitations extérieures de mêmes amplitudes et de différentes 
fréquences en fonction du temps d’excitation. On constate bien l’influence de différentes 
valeurs des fréquences d’excitation sur l’amplitude, le signe et la période de la réponse. 
Autrement dit, une proportionnalité inverse entre la valeur de la fréquence d’excitation et le 
temps de la période de réponse. C’est-à-dire, plus que la valeur de la fréquence d’excitation est 


grande, le temps de la période de réponse est court. 
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CHAPITRE V : Calcul dynamique d’une plaque rectangulaire reposant sur la 2021 
surface du semi-infini élastique de propriétés inertielles (modèle de Lamb). 
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La réponse au niveau de l’élément portant le numéro 218 


La réponse au niveau de l’élément portant le numéro 683 


Fig. V. 18 : Comparaison de la réponse de la plaque au niveau des points d’application des deux 
excitations extérieures de mêmes amplitudes et de différentes fréquences en fonction du 
temps de l’excitation 


Cas 8.3.3 : Etude comparative de la réponse de la plaque aux points d’application des 
excitations extérieures 


La figure (V.19) représente une étude comparative de la réponse de la plaque en (2D) au niveau 
des deux points d’application des excitations extérieures en fonction du temps d’excitation et 
pour différents types d’excitations appliquées sur la plaque. 


CHAPITRE V : Calcul dynamique d’une plaque rectangulaire reposant sur la 2021 
surface du semi-infini élastique de propriétés inertielles (modèle de Lamb). 


A.1 : Réponse de la plaque au niveau du point portant le numéro 218 lorsque la plaque sollicitée par : 


(1) 
(D) 


une seule excitation extérieure appliquée au niveau du point portant le numéro 218 
deux excitations extérieures de mêmes amplitudes et de différentes fréquences 
appliquées au niveau des points portant les numéros 218 et 683 


A.2 : Réponse de la plaque au niveau du point portant le numéro 218 lorsque la plaque sollicitée par : 


() 
(D) 


une seule excitation extérieure appliquée au niveau du point portant le numéro 218 
deux excitations extérieures de mêmes amplitudes et de mêmes fréquences appliquées 
au niveau des points portant les numéros 218 et 683 


A.3 : Réponse de la plaque au niveau du point portant le numéro 218 lorsque la plaque sollicitée par : 


(1) 


(D) 


deux excitations extérieures de mêmes amplitudes et de mêmes fréquences appliquées 
au niveau des points portant les numéros 218 et 683 

deux excitations extérieures de mêmes amplitudes et de différentes fréquences 
appliquées au niveau des points portant les numéros 218 et 683 


A.4 : Réponse de la plaque au niveau du point portant le numéro 218 lorsque la plaque sollicitée par : 


(1) 
(D) 


(I) 


une seule excitation extérieure appliquée au niveau du point portant le numéro 218 
deux excitations extérieures de mêmes amplitudes et de mêmes fréquences appliquées 
au niveau des points portant les numéros 218 et 683 

deux excitations extérieures de mêmes amplitudes et de différentes fréquences 
appliquées au niveau des points portent les numéros218 et 683 


B.1 : Réponse de la plaque au niveau du point portant le numéro 683 lorsque la plaque sollicitée par : 


(1) 
(D) 


une seule excitation extérieure appliquée au niveau du point portant le numéro 683 
deux excitations extérieures de mêmes amplitudes et de différentes fréquences 
appliquées au niveau des points portant les numéros 218 et 683 


B.2 : Réponse de la plaque au niveau du point portant le numéro 683 lorsque la plaque sollicitée par : 


(1) 
(D) 


une seule excitation extérieure appliquée au niveau du point portant le numéro 683 
deux excitations extérieures de mêmes amplitudes et de mêmes fréquences appliquées 
au niveau des points portant les numéros218 et 683 


B.3 : Réponse de la plaque au niveau du point portant le numéro 683 lorsque la plaque sollicitée par : 


(1) 


(D) 


deux excitations extérieures de mêmes amplitudes et de mêmes fréquences appliquées 
au niveau des points portant les numéros 218 et 683 

deux excitations extérieures de mêmes amplitudes et de différentes fréquences 
appliquées au niveau des points portant les numéros 218 et 683 


B.4 : Réponse de la plaque au niveau du point portant le numéro 683 lorsque la plaque sollicitée par : 


(1) 
(D) 


ID 


100 


une seule excitation extérieure appliquée au niveau du point portant le numéro 683 
deux excitations extérieures de mêmes amplitudes et de mêmes fréquences appliquées 
au niveau des points portant les numéros 218 et 683 

deux excitations extérieures de mêmes amplitudes et de différentes fréquences 
appliquées au niveau des points portant les numéros 218 et 683 


CHAPITRE V : Calcul dynamique d’une plaque rectangulaire reposant sur la 2021 
surface du semi-infini élastique de propriétés inertielles (modèle de Lamb). 


Cas 8.3.4 : Visualisation (3D) de la réponse de la plaque toute entière : 


ES 


La réponse de la plaque toute entière en (3D) à ce même type de chargement extérieur 


ES 


harmonique est illustrée dans la figure (V.20) à chaque instant de temps de son excitation 
exprimé par : f, =f, +ixAf. Ici le temps initial #, =0 s et le pas du temps Af =0.0002s. 


t =t, +6At t=t,+7At t=t,+8At 


Fig. V. 20 : Visualisation (3D) de la réponse de la plaque toute entière due à deux excitations extérieures 
harmoniques de mêmes amplitudes et de différentes fréquences d'excitation 


Enfin, la projection instantanée de la réponse de la plaque sollicitée par ce type de chargement 
extérieur harmonique sur le plan Oxy (isovaleurs) est illustrée sur la figure (V.21). 
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CHAPITRE V : Calcul dynamique d’une plaque rectangulaire reposant sur la 2021 
surface du semi-infini élastique de propriétés inertielles (modèle de Lamb). 


0.5 


-0.5 0.0 0.5 -0.5 


£ 0.0 
t=0s t=0,000975343s 


-0.5 00 0.5 


t =0,00273096 s t =0,00380384 s 


Fig. V. 21 : Projection instantanée de la réponse de la plaque sur le plan Oxy (isovaleurs) due à deux 


excitations extérieures harmoniques de mêmes amplitudes et de différentes fréquences 
d'excitation 


D'après la figure (V.20) on constate bien que les déplacements maximaux de la plaque sont 
obtenus au niveau des points d’application des excitations extérieures à l’instant initial. Ensuite, 
à chaque instant de l’excitation de la plaque, ces deux points de plaque manifestent de 
déplacements de différentes valeurs de l’amplitude et de différents signes. Ceci explique 
l’importance que peut jouer la valeur de la fréquence d’excitation sur le signe, l’amplitude et la 
période de la réponse de la plaque à chaque instant de son excitation. 
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CHAPITRE V : Calcul dynamique d’une plaque rectangulaire reposant sur la 2021 
surface du semi-infini élastique de propriétés inertielles (modèle de Lamb). 


Cas 8.3.5 : Etude du problème inverse de la réponse dynamique de la plaque : 


L’étude du problème inverse nous permet de déterminer les valeurs des fréquences des 
excitations extérieures agissant sur la plaque étudiée à partir de sa réponse. En effet, lorsque on 
ignore les valeurs des fréquences des excitations harmoniques extérieures et en connaissant la 
réponse de la plaque due à ce type d’excitations extérieures auxquelles la plaque est sollicitée, 
on fait recours à l’étude du problème inverse de ce problème. Cette étude du problème inverse 
nous permet de savoir les valeurs de ces fréquences d’excitations extérieures. Il faut juste faire 
la transformée de Fourrier du signal des déplacements de la plaque, à l’aide des logiciels de 
calcul Mathematica ou Matlab, ce qui aboutit à un spectre fréquentiel. L’analyse du spectre 
fréquentiel obtenu permet la localisation des piques correspondant aux valeurs des fréquences 
des excitations extérieures. 


Les deux spectres présentés concernent les cas suivants : 


Cas 1 : Signal pris au point d’application de l'excitation de fréquence © =150 Hz : 


L’application de la transformée de Fourrier (FFT) pour ce cas 1 est effectuée pour un signal 
représentant la réponse de la plaque au centre de l’élément portant le numéro 218 lorsque la 
plaque est sollicitée par deux excitations extérieures de différentes fréquences. Le spectre 
fréquentiel obtenu de cette FFT, à l’aide du logiciel Matlab, est présenté sur la figure (V.22) où 
l’on voit clairement deux piques de différentes amplitudes. La lecture des fréquences où les 
piques apparaissent indique clairement les valeurs suivantes : ©=150 Hz et © = 400 Hz. Ces 
valeurs des fréquences représentent évidemment les valeurs des fréquences des excitations 
extérieures agissant sur la plaque étudiée. 


1200 


1000 |- 150 Hz : 


400 Hz 


Magnitude de la FFT 
O 
© 
[æ] 
T 
L 


L 
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Fig. V. 22 : Spectre fréquentiel obtenu par la transformée de Fourier (FFT) du signal de réponse de la 
plaque au point portant le numéro 218 
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CHAPITRE V : Calcul dynamique d’une plaque rectangulaire reposant sur la 2021 
surface du semi-infini élastique de propriétés inertielles (modèle de Lamb). 


Cas 2 : Signal pris au point d’application de l'excitation de fréquence © = 400 Hz : 


De la même manière décrite au cas 1, on effectue la transformée de Fourier du signal 
représentant la réponse de la même plaque au centre de l’élément portant le numéro 683. Le 
spectre fréquentiel obtenu de cette FFT, à l’aide du logiciel Matlab, est présenté sur la figure 
(V.23) où l’on voit clairement deux piques de différentes amplitudes. La lecture des fréquences 
où les piques apparaissent indique aussi clairement les valeurs suivantes : ©=150 Hz et 
© = 400 Hz. Ces valeurs des fréquences représentent évidemment les valeurs des fréquences 
des excitations extérieures agissant sur la plaque étudiée. 


1200 


1000 | 400 Hz : 
800 - J 


600 - . 
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Fig V. 23 : Spectre fréquentiel obtenu par la transformée de Fourier (FFT) du signal de réponse de la 
plaque au point portant le numéro 683 


L’analyse spectrale effectuée sur les deux signaux de réponses choisis au niveau des deux points 
d’application de la charge nous a permis la location des piques des deux fréquences des 
excitations extérieures qui sont : 150 Hz et 400 Hz. La différence observée pour chaque cas 
entre les amplitudes de la magnitude de la FFT des deux piques est expliquée par le suivant : le 
point où le signal de réponse est pris indique automatiquement l’amplitude de pique la plus 
grande de l’excitation la plus proche à ce point. 


V.9. Comparaison des résultats : 


La comparaison des résultats est très importante pour s’assurer de la certitude des résultats 
obtenus et par conséquent la vérification de la fiabilité et l’authenticité de l’approche utilisée. 
Le domaine des problèmes de contact ayant une grande importance que ce soit dans les sciences 
de l’ingénierie ou dans les sciences de la physique ou autres, mais ce domaine vaste est 
relativement peu investigué vu la complexité et la difficulté mathématique que présente, comme 
le montre la recherche bibliographique. De ce fait, pratiquement aucune autre étude similaire à 
la nôtre n’est trouvée dans la littérature nous permettant de comparer nos résultats relatifs à 
cette étude dynamique de la plaque carrée reposant la surface du semi-infini élastique de 
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CHAPITRE V : Calcul dynamique d’une plaque rectangulaire reposant sur la 2021 
surface du semi-infini élastique de propriétés inertielles (modèle de Lamb). 


propriétés inertielles (modèle de Lamb). C’est pourquoi, la seule possibilité qu’on dispose pour 


le moment est de faire la comparaison de cette étude avec le cas de la même plaque reposant 
sur la surface du semi-infini élastique de propriétés distributives (modèle de Boussinesq). 


V.9.1. Comparaison des valeurs des fréquences propres : 


Le grand avantage de l’approche utilisée est de choisir le modèle du milieu élastique sur lequel 
la plaque se repose. Pour ce faire, on utilise la même approche, basant sur le principe de la 
méthode de Zhemochkin décrit avec détail dans le chapitre IV, avec la considération du modèle 
de Boussinesq (semi-infini élastique de propriétés distributives). 


Le paramètre v, ; figurant dans le système d’équations canoniques de l’approche utilisée pour 
ce modèle est exprimé par l’expression suivante [71] : 


ue n=b2 ë=c2 
ct déd7 (V.14) 


ae) TES 26 02 (8) +57) 


Tous les paramètres de la formule (V.14) sont définis dans le chapitre III. 


L’étude de la fonction de Green pour le modèle de Boussinesq nécessite l’évaluation de 
l’intégrale de la formule (V.14) ce qui aboutit à l’expression donnant les déplacements de la 
surface du semi-infini élastique de propriétés distributives. Le détail de cette intégration est 
donné dans le chapitre (IT) où le résultat de l’intégration est donné par la formule (III. 12). 


Les changements finaux dans le système d’équations linéaires sous forme matricielle finale 
issus du choix de ce modèle du milieu élastique sont : 


1-v, 1 à ; : : ; 
= x——: ce paramètre æ& apparu au niveau des termes de la matrice [A] et il est issu 
tE, 
de plusieurs transformations et simplifications mathématiques ; 
_ n=b2 £=c2 dédy 


ÿ = : l’expression finale de ce paramètre F; est donnée dans le 


n=-b26=-02 A](x — €) +(y- n) 
chapitre (I). 
Enfin, les valeurs des fréquences propres de la plaque rectangulaire reposant sur la surface du 


semi-infini élastique de propriétés distributives (modèle de Boussinesq) sont déterminées par 
la même manière que pour le cas du modèle de Lamb. 


Considérons les mêmes caractéristiques géométriques et mécaniques pour la plaque et pour le 
semi-infini élastique. L’utilisation du logiciel Mathematica nous a permis de construire un code 
de calcul approprié et son exécution a permis de trouver les résultats concernant les valeurs des 
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CHAPITRE V : 


Calcul dynamique d’une plaque rectangulaire reposant sur la 
surface du semi-infini élastique de propriétés inertielles (modèle de Lamb). 


2021 


six premières valeurs des fréquences propres de la plaque pour différents maillages. Ces 


résultats sont dans le tableau suivant : 


Sea ei © 2 © 3 © 4 © 5 © 6 
100 52,37492 | 64,46812 | 70,30563 | 74,19763 | 126,47649 | 250,13740 
121 52,30285 | 64,40745 | 70,05412 | 7401713 | 126,21661 | 249,62352 
144 52,24808 | 64,36634 | 69,86407 | 73,84385 | 125,74445 | 249,33940 
169 52,20549 | 64,33076 | 69,71694 | 73,73665 | 125,57864 | 249,03499 
196 52,17173 | 64,30510 | 69,60065 | 73,63165 | 125,29782 | 248,84404 
225 52,14453 | 64,28248 | 69,50715 | 73,56276 | 125,18598 | 248,64936 
256 52,12229 | 64,26542 | 69,43085 | 73,49440 | 125,00585 | 248,51687 
289 52,10389 | 64,25016 | 69,36775 | 73,44748 | 124,92692 | 248,38502 
324 52,08849 | 64,23825 | 69,31499 | 73,40052 | 124,80470 | 248,29001 
361 52,07548 | 64,22747 | 69,27041 | 73,36714 | 124,74696 | 248,19665 
400 52,06437 | 64,21883 | 69,23242 | 73,33352 | 124,66033 | 248,12648 
441 52,05483 | 64,21096 | 69,19975 | 73,30890 | 124,61684 | 248,05797 
484 52,04658 | 64,20450 | 69,17149 | 73,28400 | 124,55325 | 248,00479 
529 52,03938 | 64,19855 | 69,14685 | 73,26534 | 124,51968 | 247,95305 
576 52,03308 | 64,19360 | 69,12525 | 73,24641 | 124,47167 | 247,91184 
625 52,02752 | 64,18901 | 69,10621 | 73,23191 | 124,44520 | 247,87183 
676 52,01260 | 64,18512 | 69,08934 | 73,21717 | 124,40808 | 247,83929 
729 52,00823 | 64,18050 | 69,04432 | 73,20571 | 124,35685 | 247,78768 
784 52,00772 | 64,17341 | 69,03808 | 73,17601 | 124,31760 | 247,73581 
841 52,00768 | 64,17300 | 69,03799 | 73,17544 | 124,31689 | 247,73497 
900 52,00765 | 64,17299 | 69,03797 | 73,17534 | 124,31682 | 247,73490 


Tab. V. 2 : Les valeurs des six premières fréquences propres de la plaque étudiée obtenues pour 
différents maillages 


D’après les résultats apportés au tableau (V.2), on constate bien des convergences des valeurs 
des différentes fréquences propres vers des valeurs stables pour chaque fréquence propre. Cette 
convergence est obtenue grâce à l’augmentation de la taille du maillage de la plaque (à n = 900 
éléments) ce qui est en totale concordance avec le principe général des méthodes numériques. 
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La comparaison des valeurs des six premières fréquences propres de la plaque carrée reposant 
sur la surface du semi-infini élastique des deux modèles considérés est donnée dans le tableau 
(V3). 


Fréquence @1 © 2 © 3 © 4 © 5 © 6 
propre 

à de! 5193084 | 6414663 | 68.097462 | 73.13283 | 124.28131 | 247.71888 

Modèle de | 53 00765 | 6417299 | 69.037097 | 73.17534 | 124.,31682 | 247,73490 

Boussinesq 


Tab. V.3 : Comparaison des valeurs des six premières fréquences propres de la plaque étudiée trouvées 
par le modèle de Lamb et le modèle de Boussinesq 


L’analyse des données du tableau (V.3) relatives aux valeurs des fréquences propres de la 
plaque carrée reposant sur la surface du semi-infin1 élastique des deux modèles montre bien une 
petite différence des valeurs données par les deux modèles. Cette différence doit être existée et 
elle est expliquée par la nature de chaque modèle. Le modèle de Boussinesq néglige certains 
facteurs naturels comme le caractère ondulatoire des déplacements que le modèle de Lamb tient 
en compte. Donc le modèle de Lamb est plus proche de la réalité puisque il tient compte de 
plusieurs facteurs naturels influençant sur le comportement dynamique des structures en 
général. De ce fait, les valeurs des fréquences propres de la plaque étudiée considérées plus 
réelles et plus naturelles sont celles données par le modèle de Lamb. 


V.9.2. Comparaison de la réponse dynamique de la plaque : 


Pour s’assurer de la certitude des résultats trouvés concernant la réponse dynamique de la 
plaque rectangulaire reposant sur la surface du semi-infin1 élastique de propriétés inertielles par 
l’approche proposée, on doit les comparer avec ceux de la même plaque reposant sur la surface 
du semi-infini élastique de propriétés distributives. Pour ce faire, la détermination de la réponse 
dynamique de la plaque carrée reposant sur la surface du semi-infini élastique de propriétés 
distributives (modèle de Boussinesq) se fait par la même approche proposée et de la même 
manière décrite pour le modèle de Lamb. Considérons la structure étudiée ayant les mêmes 
caractéristiques géométriques et mécaniques prises pour le cas du modèle de Lamb et que la 
plaque est sollicitée par une seule excitation extérieure harmonique : 


n=900 
À; = >, W,P> =W 305 Pos 


p=l 
F, = 10000 N : l’amplitude d’excitation ; 


© = 400 HZ : la fréquence d’excitation ; 
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205 : numéro de l’élément où l’excitation extérieure harmonique est appliquée. 


Pour ce cas de modèle, la réponse de la plaque carrée se détermine par la formule suivante [71] : 


1-v°, | 
Pal ete 1. V.15 


L’exécution du code de calcul élaboré permet la détermination de la réponse de la plaque 
représentant les valeurs de ses déplacements verticaux au niveau des centres de tous les 
éléments de sa discrétisation. Les déplacements verticaux donnés par les deux modèles au point 
d'application de l’excitation extérieure en fonction du temps de l’excitation sont illustrés sur la 
figure (V.24). 


v (mm) 


0.0015 0.0025 


Fig. V. 24 : Comparaison de la variation du déplacement vertical de la plaque au point d’application de 
l’excitation extérieure sur un intervalle de temps de son excitation donnée par : 


(go) semi-infini élastique de propriétés inertielles (modèle de Lamb) ; 
(D) semi-infini élastique de propriétés distributives (modèle de Boussinesq). 


La remarque que l’on peut soulever via cette comparaison de réponse de la plaque pour les deux 
modèles est que le modèle de Boussinesq donne des valeurs des déplacements de la plaque un 
peu élevées par rapport à celles données par le modèle de Lamb. Cette petite différence des 
valeurs donnée par les deux modèles doit être existée et elle est expliquée par la différence de 
nature des deux modèles. En outre, cette réduction de valeur du déplacement donnée par le 
modèle de Lamb est aussi expliquée par la tenue en compte d’un facteur essentiel qui est le 
caractère ondulatoire des déplacements, ce qui justifie encore la certitude des résultats de ce 
modèle. 
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V.10. Cas particulier : 


Pour montrer l’importance et l’effet du caractère ondulatoire des déplacements sur l’étude 
dynamique des structures en général et d’une manière explicite, considérons la structure étudiée 
précédemment ayant les mêmes caractéristiques mécaniques, mais de différentes 
caractéristiques géométriques. La zone de contact entre la plaque carrée et la surface du semi- 
infini élastique de propriétés inertielles sera élargie considérablement pour montrer 
concrètement ce phénomène. Tenons le même maillage de la plaque et considérons qu’elle est 
sollicitée par deux excitations extérieures harmoniques exprimées par : 


(V.16) 


Pie = F,cos(2X 7 x © Xf) 
Pie = PE, coS(2X 7 x ©, Xf) 


Où : 


L’excitation extérieure est appliquée au centre des éléments portant les numéros p=218 et 
p = 683 ; 


L’amplitude des deux excitations est : P, =10000 N ; 
Les fréquences d’excitation sont : @, =150 Hz ; ©, = 400 HZ ; d’où w= 4150 +400 =427,2Hz : 


Les caractéristiques géométriques et mécaniques de la plaque et du semi-infini élastique sont : 


a=b=10 m : longueur et larguer de la plaque ; 

h=0.2m : épaisseur de la plaque ; 

v =1/3 : coefficient de Poisson de la plaque ; 

E=2x10" N/m° : module d'élasticité de la plaque ; 

E,=25x10 N/m° : module d'élasticité du semi-infini élastique ; 
v, =1/3 : coefficient de Poisson du semi-infini élastique ; 

p =1500 Kg / m° : densité du semi-infini élastique. 


L'application de l’approche proposée pour cette structure pour le modèle de Lamb a abouti aux 
résultats présentés graphiquement en (3D) sur la figure (V.25). 
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t =t, +6At t=t,+7A t=t,+8At 


Fig. V. 25 : Visualisation (3D) de la réponse de la plaque étudiée due à deux excitations harmoniques 
de mêmes amplitudes et de différentes fréquences d'excitation 


D’après la figure V.25, nous pouvons voir clairement le phénomène du caractère ondulatoire 
des déplacements lorsque la structure ayant des dimensions assez grandes. Ceci montre que ce 
phénomène existe toujours et se présente même pour les structures ayant de petites dimensions. 
Autrement dit, plus les dimensions de la structure deviennent grandes, plus le phénomène du 
caractère ondulatoire de ses déplacements devient apparent. 


V.11. Conclusion : 


Ce chapitre est consacré à l’application de l’approche semi-analytique décrite dans le chapitre 
(IV) pour l’étude dynamique d’une plaque rectangulaire reposant sur la surface du semi-infini 
élastique de propriétés inertielles (modèle de Lamb). L’approche utilisée est formulée de telle 
sorte qu’elle permette en premier lieu la détermination des forces de réaction dans la zone de 
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contact qui représentent le phénomène d'interaction entre la plaque et la surface du semi-infini 
élastique. Ensuite, et de façon numérique elle permet la détermination des autres entités 
physiques recherchées en appliquant les différentes lois de la RDM connues dans ce domaine. 
Les entités physiques recherchées dans cette étude dynamique de la plaque en interaction avec 
le semi-infini élastique de propriétés inertielles sont les valeurs des fréquences propres et les 
modes propres de la plaque, ainsi que sa réponse dynamique due à l’action de plusieurs types 
d’excitations harmoniques extérieures. 


La formulation finale de l’approche est faite sous forme matricielle et a nécessité la construction 
des codes de calcul appropriés pour chaque cas de calcul en utilisant le logiciel de calcul formel 
et numérique Mathematica. La précision des résultats nécessite un maillage de la structure très 
fin ce qui rend l’exécution du calcul très difficile et très long. Donc, une autre nécessité 
s’impose, celle d'utilisation des stations de calcul puissantes. 


Tous les défis rencontrés dans cette étude dynamiques ont été surmontés grâce au travail 
laborieux et continu sur l’élaboration des codes de calcul et enfin de compte, cette persévérance 
a abouti à la détermination des valeurs des fréquences propres de la plaque, ses modes propres 
et ses réponses dynamiques dues aux différents types des excitations extérieures harmoniques. 
La comparaison des résultats de calcul avec un autre modèle a montré, non seulement, leurs 
certitudes, mais aussi la fiabilité de l’approche utilisée. La petite différence trouvée entre les 
deux modèles doit être existée vu leur différence de nature et de principe. Le jugement que le 
modèle de Lamb est plus proche de la réalité est fondé sur le fait qu’il tient en compte de 
plusieurs facteurs comme l’inertie du milieu élastique et le caractère ondulatoire des 
déplacements. Tous les résultats de calcul obtenus par l’application de cette approche sont en 
totale cohérence avec la réalité des choses et que leurs précisions sont justifiées et argumentées 
par le maillage fin utilisé. 
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Conclusion générale : 


L'étude apportée dans ce travail de thèse de doctorat concernant le calcul dynamique de plaques 
reposant sur la surface d’un milieu élastique de n'importe quel type est effectuée par une 
approche semi-analytique basant sur la méthode des forces pour le calcul des structures 
hyperstatiques connue dans les sciences de la résistance des matériaux. Cette étude dynamique 
comporte la détermination des fréquences propres et des modes propres de la plaque étudiée, 
ainsi que sa réponse dynamique due à l’action de différents types de charges harmoniques 
extérieures. 


L’application de l’approche proposée pour le calcul dynamique de la structure étudiée (plaque 
rectangulaire reposant sur la surface du milieu élastique de propriétés inertielles) nécessite en 
premier lieu sa discrétisation en un nombre d'éléments rectangulaires identiques. Le nombre 
des éléments de discrétisation de la structure étudiée est inversement proportionnel à l’erreur 
de calcul et donc proportionnel à la précision des résultats recherchés. Le contact continu entre 
la plaque et la surface du semi-infini élastique de modèle de Lamb est remplacé par un contact 
partiel au niveau des points situant aux centres de gravité des éléments de discrétisation de la 
structure étudiée. Ce contact partiel est assuré par des efforts de liaison qui représentent le 
phénomène de l’interaction entre la plaque et la surface du semi-infini élastique de propriétés 
inertielles. Donc, l’approche de calcul proposée est aménagée sous forme matricielle finale dont 
le vecteur des inconnues représente les valeurs des efforts de liaison dans la zone de contact. 
Une fois les forces réactives sont connues, le reste des autres entités physiques, à savoir les 
valeurs des fréquences propres et les modes propres de la plaque et sa réponse dynamique etc., 
se déterminent par l’application des différentes formules de la théorie d’élasticité. Parmi les 
avantages de l’approche proposée pour l’étude statique ou dynamique des différentes structures 
en interaction avec des milieux élastiques, on peut citer l’analyse statique ou dynamique de 
poutres et de plaques de différentes formes géométrique et de différentes caractéristiques 
mécaniques reposant sur la surface d’un milieu élastique de n’importe quel modèle. 


Les deux grands challenges rencontrés lors de l’élaboration de cette approche sont l’étude de la 
fonction de Green définissant les déplacements verticaux de la surface du semi-infini élastique 
de propriétés inertielles et la détermination des déflexions de la plaque. Le premier grand 
challenge est lié au type du modèle du semi-infini élastique choisi qui tient compte de son inertie 
et du caractère ondulatoire des déplacements ce qui a nécessité l’évaluation de certaines 
intégrales de complexité mathématique majeure dont certaines intégrales sont associées à de 
problèmes de singularité. Le second grand challenge est relatif à la résolution de l’équation 
différentielle des déflexions de la plaque étudiée en satisfaisant ses conditions aux limites. Le 
premier challenge est surmonté grâce au changement de variables passant du système cartésien 
au système polaire et aussi par l’approximation de certaines fonctions spéciales par de séries, à 
l’aide du logiciel de calcul formel Mathematica, ce qui a abouti à une solution semi-analytique. 
Le deuxième challenge est surmonté grâce à l’adoption de la solution générale de Clebsch et 
l’application de la méthode énergétique de Ritz pour satisfaire les conditions aux limites de la 
plaque étudiée. Après avoir accompli l’étude de ces deux grandes tâches principales, ainsi que 
d’autres tâches auxiliaires qui sont toutes données sous formes d'algorithmes simplifiés et 
adaptés pour être éventuellement utilisés par les centres de calcul d'ingénieur et de recherche, 
l’approche finale est ainsi formulée sous forme matricielle bien adaptée pour différentes 
applications. 


Les résultats du calcul dynamique obtenus par l’application de l’approche proposée, après avoir 
introduit les caractéristiques géométriques et mécaniques de la plaque et du semi-infini 
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élastique de propriétés inertielles dans le système matriciel final, sont très satisfaisants. A noter, 
qu’un autre défi s’est apparu lors du lancement du calcul des valeurs des fréquences propres de 
la plaque, représentant les racines de l’équation du déterminant de la matrice du système final, 
lié au temps de calcul. Autrement-dit, plus le maillage de la plaque est fin plus le temps 
d’exécution est important et plus la précision de la valeur de la fréquence propre est bonne. Ce 
défi est surmonté grâce à l’utilisation du Plateau Technique des Calculs Intensifs de l’université 
8 Mai 1945 Guelma ce qui a nécessité plus de sept mois de calcul continu. 


Pour s’assurer de l’authenticité et de la certitude des résultats obtenus par l’application de cette 
approche, on a procédé aux vérifications suivantes : 


L comparaison des valeurs des fréquences propres de la structure étudiée (plaque reposant 
sur la surface du semi-infini élastique de propriétés inertielle) avec la même plaque 
reposant sur la surface du semi-infini élastique de propriétés distributives (modèle de 
Boussinesdq) ; 

IL vérification des valeurs des fréquences propres de la plaque étudiée en utilisant le principe 
de superposition de la méthode modale. 

IT. Comparaison de la réponse dynamique de la plaque en considérant la même plaque 
reposant sur la surface du semi-infini élastique de propriétés distributives (modèle de 
Boussinesq) 

IV. Application du principe des problèmes inverses en vue de localiser les fréquences 
d’excitation extérieures de la plaque à partir des signaux de sa réponse dynamique. 

Toutes les comparaisons et les vérifications faites sont satisfaisantes ce qui témoigne de la 

certitude des résultats obtenus et par conséquent la fiabilité de l'approche utilisée. 


Les résultats de recherche de cette thèse de doctorat ont abouti à une publication internationale 
de classe À, à trois communications internationales et une autre publication internationale en 
cours de rédaction. 


Enfin, l'étude apportée dans ce travail de thèse de doctorat peut être utilisée comme base pour 
traiter certains cas plus compliqués qui ne sont pas encore abordés. Il s’agit principalement de 
la prise en compte de l'amortissement dans la zone de contact ; la tenue en compte des forces 
de frottement dans la zone de contact ; l’analyse des structures ayant des formes géométriques 
plus compliquées ; l’analyse des structures constituées de matériaux orthotropes ; l’analyse des 
structures se trouvant à l’intérieur des milieux élastique (en interaction totale) etc. 
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Annexe : O1 : 


Le calcul de l’un des intégrales de la fonction de Green lorsque r =0 (apportée au chapitre ID) 
est détaillée de la manière suivante : 


n=b2 &£=c2 


Î (Lrnnerseomenicrro L{m) ,{L1+m},- : kr ) dédn * 


n=-b2 £=-c2 


J{rsperseomenicPrO| (n) : {1,1 + m} , = x 2 rdr = 


= : r Gamma [1 + m| HypergeometricPFORegularized L{m) ; {1 +m, 2} ,— : kr? | 


2 


O1 [r | = : rGamma [1 + m| HypergeometricPFORegularized L{m) - {1 + M, 2} >— 5 ée : 


Faim O1 ee feuel|- 


cos [e] 


= = HypergeometricPFQ L{m) ,{2,1+m},- : c2?k° sec [el | sec[o] 


Series |HpergeomericPFO L{m) , (2, 1+ m} ,— s ET | f {£, 0, )] = 


(c2 km)" . 
sn À 
__. 
| ] (22) te) d 
cos[p| 8(1+m) 
Lt. c2°k*mtan[p] c2*k*msec| J'tan[ ] 
a un(e ee. . a) 
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ee-1e pe ft re el 


FullSimplify [e arr 2] - 20] x ] _ b2c2[24 - Es ; . +32 )K°) m) | 


n=b2 E=c2 


[ LLrtneneomenricPro| {ri : {1,1 + m} ; = KT” ]) dédn = 


n=-b2 £=-c2 


b2c2(24-(-24+(b2° +3c2)k° }m) 
6(1+m) | 


Le détail du calcul de l’un des intégrales de la fonction de Green lorsque r Æ 0 (apportée au 
chapitre II) s’est fait de la manière suivante : 


Fig.A.1 : Division de l'élément chargé en 16 sous-éléments 
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Où : 
(x; y; ;c2,;c4 ;b2;b4):Les coordonnées de l’élément chargé; 
Avec : 


b b : : 
be PT Fe ee) C2=X, 55 


b 3b 3b b 
ere Dr Ho ha PAN Se 
M2 —1];; 3 476 —Îs; > 4/0 — To; 3 Aa = Miss 
73 =; M = 15; DCE Ms = 13; 

Ta =; Ms 1]; To = Mo; Me — 3° 
C 3c 3c C 
=Cc2+—; =C2+—; = c4——; =c4——; 
a=a+S [rot £-a-À lee 


Es =8 dé =; : = RE 
Lo = 4; 0 = &; di = 35 É = 64 
3 = Gi; ia = Ga; Lis = &35 6 = Sa 


y=b2 x=c2 


Î [rperseomeniePro] {fu Le (6) +0) ) | - 


y=-b2 x=-c2 
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2 2\] 2 2 
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Annexe : 02 : 


Le calcul des coordonnées de chaque élément est très important pour l’application de l’approche 
utilisée. Le calcul des coordonnées est plus exposé aux erreurs surtout lorsque le maillage utilisé 
est plus fin. Si on prend un exemple lorsque la discrétisation de la structure est n = 900 
éléments, alors le nombre de coordonnées nécessaires pour les différentes combinaisons entre 


éléments est : 6x (900) = 4860000 combinaisons. 


C’est pourquoi, il est indispensable de construire un code de calcul pour la détermination des 
coordonnées des éléments pour éviter les erreurs des calculs antérieurs. 


Pour une discrétisation de n éléments, les coordonnées (x; »:) sont calculées de la manière 


suivante : 


ve 


2a 
a =— ; 
Avec : A É 
b=— ; 
s 
Où : 


a, b : la demi-longueur et la demi-largeur de la plaque respectivement ; 
a;, b, : longueur et largeur de l'élément i respectivement ; 


L'indice ( j) est déterminé par la manière suivante : 
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L 
=], ie[t Vh | Vn 
5 ie] Vn+1:24h |; = 2, 14 T<2 
23 ie[2Vn+134h |; L 
n 


Î= 


= Vn, ie[(Nn-1)Vn + n 


L'indice (p) est déterminé par la manière suivante : 


“1 Left n+k2n+tesn-(Vn 1): 
=$, ie[2 Vn+2:2n+2.;n-{Vn-2)| 

p=l=3, Len +3200+%;m- (m3) | 0. 
Fe ie[ Vn;2Vn;3Vn;.;n |; 


Donc, les coordonnées (x;: y,) deviennent : 


Het X =, x F0 
pa [Xp 4 Xr © o 
vi} =; aVr}r2 = 6, Xp = kr 
Era = é 3 he = 3 fesse Xe = Êr 
nr) ET n{r-2) G Xe 
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Vu) 1 fr 
Y = ru Te 1j 
V2 1 fr ” 
= rs G) 
Ja Un Apte ere (V4) pi 
ph nl à 


Prenons un exemple d’une plaque discrétisée en 25 éléments rectangulaires identiques : 
24 


n=25; a=b=1lm; a =b,= 


Fig. A.2. Exemple d’une plaque discrétisée en 25 éléments rectangulaires identiques 


Si nous prenons le cas du déplacement (v,19) au centre de l’élément 7 causé par la charge 


appliquée au centre de l’élément 19 de coordonnées vie c4; b2; b4). Toutes ces 


coordonnées sont calculées de la manière suivante : 
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(x, =&)=-a+(2x px )=-1+(2x2-x(02)--04 


(y =m)=b- f(x) 1)x (2) 1-((2x2-1)x(0,2))=0,4; 
(xo=É,)=-a+(2xp- fs) 1+(2x4-1)x(0,2)=0,4; 


(y =72)=b- fexi-n) 1)x (2) ((2x4-1)x(0,2))=-0,4; 


9 


b. b. a. a, 
bA=y,+—; b2=y,-—;cA4=x,+—; c2=x, —-—; 
Yo 2 Yo 2 7 19 

Finalement, on trouve : 


x,=—0,4 y,=0,4 x,=0,4 y,=-0,4 b4=-0,2; b2=-0,6, c4=0,6, c2=0,2; 
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Annexe : 03 : 


La matrice [A] contenant les éléments de sa première et de ses quatre dernières colonnes, ainsi que ses deux premières et ses quatre dernières lignes lorsque 
la plaque étudiée est discrétisée en n=9 éléments est la suivant : 


e dred 2 2 
Yi +a(F, -MO (EM + EM 2 + EM + EM a + Es + V6 + Ein + Bis + EYYo) Yo +al(r,s MO (FYY, HE 3 + BV 3 + Fa + Es + F6 + Fo Vi + EF oŸ Vs + FsYYo) 7 à 
3 3 
2 0 
PP +al8, -MO (EM, + 2 + 3 + Ed Es + F6 + Er + Rs + EYVo) Yo +a(F Mo (FYY,; HE 5 + Es + Fa + Es + F6 + F7 + Fois + Ln2 2) = 
3 
: as 3 2 2 
a+ a(Rs MO (Eh + Eh Ed + Eos + Fos + F6 + For + Rd Rs + Ro) Po +a(Rs MO (Fi + Eh + Bols + Eh + Eos + Fos + Eh + Fit E bo) 7 . 
3 3 
1 0 
2F, 2F 2F, 2F, 2F, 2F ,.(Û2F 2F, 2F, 2F, 2F, 2F, at 
Mac Mn, SIOE +0F. +08, ga Ag on Ma 19, 2729. 2739 +0, +085 +080 — 7.9 Ag  Æfo9 + mo 
3 3 3 3 3 3 3 3 3 3 
2F, 2F, 2F. 2F, 2F, 2 2F; 2F, 2F, 2F, 2E, 2F, 0 1 2 
Mo] = 275 OF, 31 4,1 OF, 6.1 Ce OF, + 9,1 ol ES à 0.F,9 + 39 ao OF, sg APT, OF+ 9,9 a mo 
3 3 3 3 3 3 3 3 3 3 
2(F 4 - 2 LR à LP à 0 0 
1-Maw (F, : Fi +E, : Ln +F; : Fa +E,; +k; +F,) 1-Maw (Es +E,s ; Fo L Fo +Fs u F9 +, F9 ; Es) 
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